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ABSTRACT 


A computer  model  is  developed  for  simulating  reverse  bias 
thermal  second  breakdown  (TSB)  transients  in  thin  film  diodes.  The 
model  performs  a one-dimensional  electrical  and  a two-dimensional 
thermal  simulation.  Simulations  are  performed  up  to  the  onset  of 
the  TSB  transition  to  a high  conductance  state.  This  condition  is 
defined  as  a maximum  junction  temperature  of  700  °K.  The  model  is 
driven  by  a constant  current  source  and  features  temperature  and 
electric  field  dependent  avalanche  ionization  coefficients,  temper- 
ature, electric  field  and  doping  level  dependent  mobilities,  and 
depletion  region  space  charge  effects.  Simulations  are  defined 
through  271  parameters  which  specify  diode  design,  thermal  conduct- 
ivity perturbations,  and  control  of  the  simulation.  The  program 
generates  graphic  output  and  requires  short  run  times. 


ACKNOWLEDGEMENTS 


The  author  wishes  to  acknowledge  the  technical  guidance 


provided  by  Dr.  David  Mathews  of  the  U.S.  Army  Missile  Command 


the  assistance  by  Larry  Koos,  Lynn  Bilbo,  and  Mark  Halfacre  in 


performing  this  task,  the  efforts  of  the  electrical  engineering  office 


wife  Jennye  for  enduring  yet  another  of  my  many  endeavors 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES 


FIGURE 

1.  Silicon-On-Sapphire  (SOS)  Diode  Structure  . 

2.  Flowchart  for  Simulation  of  the  TSB 

Transient  

3.  Electrical  Analog  for  the  Diode  Thermal 

Model  

4.  Electrical  Analog  for  the  Quasi-Two- 

Dimensional  Header  Thermal  Model  

5.  Electrical  Analog  for  the  Two-Dimensional 

Header  Thermal  Models  

6.  Electrical  Analog  for  the  Combined  Diode 

and  Header  Thermal  Models  

7.  Diode  Doping  Profile  

8.  Flowchart  for  the  Diode  Electrical  Model 

9.  SPARM  and  DPARM  Parameter  Sets  for  a Thermal 
Second  Breakdown  Demonstration  Simulation  . 

10.  Comparison  of  TSB  Delay  Times  for  the  Quasi 
Two-Dimensional  Substrate  Thermal  Model  for 
Zero  and  Twelve  Substrate  Node  Points  . . . 

11.  Comparison  of  TSB  Delay  Times  for  Constant 

and  Blocking  Boundary  Conditions  on  the 
Diode  Temperature  Profile  

12.  Comparison  of  TSB  Delay  Times  for  Zero 

and  None  Zero  Bulk  Region  Electric  Field 
Values  

13.  Comparison  of  TSB  Delay  Times  for  Quasi 

and  Full  Two-Dimensional  Substrate  Thermal 
Models  

14.  Maximum  Junction  Temperature  Versus  Time 

for  Several  Different  Driving  Current 
Values  and  a Full  Two-Dimensional 
Substrate  Thermal  Model  . 


vi 


LIST  OF  FIGURES  (Continued) 


FIGURE  PAGE 

15.  Delay  Time  Versus  Epitaxial  and  Upper 
Region  Doping  Concentration  for  Two 

Different  Driving  Current  Values  48 

16.  Diode  Maximum  Temperature  at  Time  Step 
Number  Ninety  Versus  Maximum  Avalanche 

Breakdown  Integral  Error  (AERMAX)  50 

17.  TSB  Delay  Time  Versus  Avalanche  Ionization 
Coefficient  Thermal  Dependence  Factor 

(ALFATD) 51 

18.  Maximum  Junction  Temperature  Versus  Time 
for  Several  Different  Values  of  Avalanche 
Ionization  Coefficient  Thermal  Dependence 

Factor  (ALFATD)  53 

19.  Maximum  Junction  Temperature  Versus 
Time  for  Several  Different  Values  of 
Avalanche  Ionization  Coefficient  Thermal 

Dependence  Factor  (ALFATD)  , . 54 

20.  Convergence  Characteristics  for  Header 

Thermal  Model  DHT2D  55 

21.  Convergence  Characteristics  for  Header 

Thermal  Model  DHT2D  56 

22.  Convergence  Characteristics  for  Header 

Thermal  Model  DHT2D1  57 

23.  Convergence  Characteristics  for  Header 

Thermal  Model  DHT2D1  58 

24.  A Comparison  of  Convergence  Characteristics 
for  Header  Thermal  Models  DHT2D  and 

DHT2D1 59 

25.  Convergence  Characteristics  for  header 

Thermal  Model  DHT2D2  60 

26.  Convergence  Characteristics  for  Header 

Thermal  Model  DHT2D2  61 


vii 


LIST  OF  TABLES 


Table 

Page 

1 

SIMULATION  CONTROL  PARAMETER  SET 

(SPARM) 

. . . 30 

2 

DIODE  DESIGN  AND  PHYSICAL  PROPERTY 

PARAMETER  SET  (DPARM)  

. . . 33 

3 

SPECIFICATIONS  FOR  TSB  SIMULATION 

SERIES  

. . . 37 

viii 


INTRODUCTION 


A diode  model  is  ueveloped  for  simulating  TSB  (thermal 
second  breakdown)  in  SOS  (silicon-on-sapphire)  type  diode  structures. 
The  model  is  capable  of  evaluating  the  dependence  of  TSB  on  semi- 
conductor parameters  and  diode  design.  Diode  and  simulation 
specifications  for  simulations  are  easily  changed  to  facilitate 
investigation  of  different  diode  designs.  Run  times  of  approximately 
two  minutes  make  investigation  and  characterization  of  different 
designs  economically  attractive.  These  features  make  the  model 
particularly  useful  for  preliminary  investigations  of  different  diode 
designs.  Once  the  general  characteristics  of  these  devices  have  been 
determined,  more  comprehensive  models  [1,2]  that  require  substantially 
longer  run  times  can  be  used  to  further  resolve  the  observed  TSB 
behavior. 

The  SOS  type  diode  structure  [1,6]  consist  of  a thin  layer  of 
silicon  deposited  on  a sapphire  substrate  or  header  (for  the  purposes 
of  this  report  the  terms  substrate  and  header  are  used  interchangeably 
to  identify  the  material  which  supports  the  deposited  semiconductor 
thin  film).  Fig.  1 shows  a typical  SOS  diode  configuration.  The 
diode  model  for  this  structure  is  broken  down  into  three  distinct 
but  closely  coup]ed  models:  diode  electrical,  diode  thermal  and 
header  thermal  models. 

The  significant  electrical  effects  leading  to  the  onset  of  the 
TSB  transition  are  assumed  to  be  one-dimensional  and  along  the  diode 
axis.  The  TSB  transition  time  from  a high  voltage  state  to  a post-TSB 


3 


low  voltage  state  is  assumed  to  be  short  in  comparison  with  the  delay 
time  required  to  achieve  the  critical  TSB  transition  temperature. 

Under  these  conditions,  the  transient  time  required  for  the  diode 
to  obtain  the  critical  temperature  represents  the  TSB  delay  time. 

A value  or  700  °K  has  been  chosen  for  the  critical  temperature 
for  SOS  diode  structures.  This  value  was  predicted  by  simulations 
with  a more  comprehensive  diode  model  [1],  The  electrical  model 
features  temperature,  impurity  and  electric  field  dependent  mobilities, 
temperature  and  electric  field  dependent  avalanche  ionization  coef- 
ficients, and  temperature  dependent  bulk  region  electric  fields. 

The  diode  and  header  thermal  models  combine  to  account  for 
thermal  conduction  along  the  diode  axis  and  perpendicular  to  this 
axis  into  the  substrate.  The  two  model  concept  reduces  model 
complexity  and  maintains  compatibility  with  the  more  comprehensive 
diode  model  developed  previously  [1],  As  a consequence  of  the  thin 
film  structure  of  the  semiconductor,  the  diode  thermal  model  is 
quasi-two-dimensional . This  model  accounts  for  heat  generation  and 
thermal  conduction  along  the  diode  axis  and  thermal  conduction  into 
the  substrate.  Four  different  header  thermal  models  were  developed 
for  thermal  conduction  through  the  header.  One  of  these  models  is 
quasi-two-dimensional  and  was  developed  previously  [5].  The  three 
new  models  feature  two-dimensional  thermal  conduction  and  differ 
only  in  the  numerical  implementation  of  the  energy  continuity 
equation. 

A flowchart  for  the  diode  model  is  shown  in  Figure  2. 


DIODE  THERI-IAL  MODEL 


Quasi-two-dimensional  thermal  conduction  for  the  diode  is  modeled  by 
the  one- dimensional  energy  continuity  equation  which  may  be  written 
as : 

CDCD-!f  ■ ICD0+  iJEl  - ‘ (1> 

where : 

PD  - diode  density 
Cp  - diode  specific  heat 
T - diode  Temperature 
t - t ime 

Kp  - diode  thermal  conductivity 
x - position 
J - current  density 
E - electric  field  intensity 

$ - heat  loss  by  mechanisms  other  than  conduction  along 

the  diode  axis. 

Heat  loss  into  the  diode  header  through  thermal  conduction  is  accounted 
for  through  the  $ term  to  yield  quasi-two-dimensional  thermal  conduction. 


<t>  is  defined  as: 
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where: 

Ky  - header  thermal  conductivity 
- diode  film  thickness 

Ay  - header  node  spacing  perpendicular  to  diode  axis 
T(N)  - diode  temperature  profile 
T(N,1)  - adjacent  header  temperature  profile 


An  electrical  analog  of  this  thermal  model  is  shown  in  Figure  3. 
Thermal  conduction  from  the  diode  into  the  header  is  a function 
of  the  temperature  difference  between  the  respective  diode  and  header 
node  points.  The  numerical  algorithm  for  the  diode  thermal  model 
was  developed  by  combining  equations  (1)  and  (2),  and  applying  finite 
difference  techniques  in  a fully  implicit  formulation  [6,7],  This 
procedure  yields  the  following  system  of  linear  equations: 


At  A T(N  - 1)S+1  + [AtAx2  B - (Ax2  +2  At  A)]  T(N)S+1 


+ At  A T(N  + 1)' 


Ax2  [T(N)S  + At  C(N) ] + B T(N,1)  Ax2At 


where : 


PDCDAyXDT 


C(N)  = 


present  iteration  number 


- next  iteration  number 


diode  node  number,  1 - N - NN 


total  number  of  diode  nodes 


- time  step  size 

distance  between  adjacent  nodes  along  diode  axis 


f 


! I 


or, 

A1  T(N-1)S+1  + A2  T(N)S+1  + A3  T(1W-1)8+1  - A4  (4) 

where : 

A^  = AtA 

A2  E [At  Ax2  B - (Ax2  + 2At  A)] 

A^  = AtA 

A4  E-Ax2[t(N)s  + At  C(N) ] . 

The  second  order  energy  continuity  equation  requires  two 
boundary  conditions  on  temperature.  The  model  assumes  constant 
ambient  temperature  boundary  conditions  by  default.  Blocking  or 
insulating  boundary  conditions  may  be  specified  through  the  appropriate 
simulation  parameter. 

The  resulting  system  of  NN-2  linear  equations  must  be  solved 
for  each  iteration  to  generate  an  improved  approximation  for  the 
diode  temperature  profile  at  the  next  point  in  time.  Iterations 
are  performed  until  a specified  maximum  rmS  (root-mean-square)  change 
in  the  diode  temperature  profile  is  achieved  between  successive 
iterations,  or  a specified  maximum  number  of  iterations  are  performed. 

With  the  completion  of  each  iteration  sequence,  the  diode  temperature 
profile  is  advanced  by  one  time  step.  The  iterative  sequence  between 
time  steps  is  required  only  if  the  energy  continuity  equation  coefficients 
are  temperature  dependent.  Although  the  computer  program  has  been 
formulated  with  the  full  iterative  capability,  the  energy  continuity 
equation  coefficients  are  presently  assumed  constant.  Hence,  sim- 
ulations presently  require  single  iterations  in  the  diode  thermal  model. 
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Perturbations  in  the  diode  and  header  thermal  conductivities  may  be 
specified  through  the  diode  simulation  parameters.  Two  perturbation 
factors  are  assigned  to  each  of  the  diode  node  points.  One  is 
associated  with  thermal  conduction  between  diode  node  points.  The 
other  is  associated  with  thermal  conduction  between  diode  node  points 
and  the  respective  header  node  points.  In  both  cases,  the  perturbation 
factors  are  multiplied  by  the  respective  thermal  conductivity  values 
to  yield  position  dependent  thermal  conductivity  values  for  the  diode 
thermal  model.  The  perturbation  factors  are  assigned  default  values 
of  one.  Any  of  the  perturbation  factors  can  be  redefined  through  the 
simulation  parameters.  The  demonstration  execution  of  the  diode 
model  presented  in  appendix  A has  a diode  thermal  conductivity  per- 
turbation specified. 
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HEADER  THERMAL  MODELS 


The  computer  program  offers  four  different  header  thermal  models. 

The  first  model  was  developed  earlier  [5]  and  is  a quasi-two-dimensiional 
thermal  model  identified  as  DHTEMP.  This  model  features  one-dimen- 
sional thermal  conduction  from  each  of  the  diode  temperature  node 
points  through  the  header  to  the  ambient  temperature  heat  sink.  No 
thermal  conduction  is  allowed  between  these  one-dimensional  thermal 
paths  through  the  header.  An  electrical  analog  of  this  header 
thermal  model  is  shown  in  Figure  4. 

Implementation  of  this  model  is  very  similar  to  the  diode 
thermal  model.  However,  since  only  thermal  conduction  with  heat 
storage  occurs  through  the  header,  the  energy  continuity  equation, 
equation  (1),  reduces  to: 


o c 
mH  H 


3 T 
3 t 


“h 


3 T 
3x2 


(5) 


header  density 
header  heat  capacity 
header  thermal  conductivity 

The  boundary  temperatures  for  the  model  are  specified  as  the  cor- 
responding diode  temperature  and  ambient  temperature,  respectively. 

The  numerical  algorithm  for  solving  equation  (5)  is  developed 
using  the  same  techniques  as  used  previously  for  the  diode  thermal 
model.  The  resulting  system  of  equations  can  be  written  as: 


where : 


“H  " 


1 

i 
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A1  T(N,M-1)S+1  + A2  T(N,M)S+1  + A3  T(N,M+1)S+1  - A 4 (6) 


where : 


T(N,M) 


header  temperature  profile  beginning  at  diode 
temperature  node  N,  (0<  N<  12). 


- header  node  spacing. 


and, 


A = At  A 

1 

A2  = (Ax*  + 2 At  A) 

A3  - At  A 

a4  * - Ax*  T(N,M)S 

After  the  diode  temperature  profile  has  been  updated,  the  above 
algorithm  is  applied  to  each  header  temperature  profile  to  generate 
the  new  header  temperature  profile  for  the  next  point  in  time. 

The  remaining  three  header  thermal  models  represent  two-dimen- 
sional thermal  conduction  with  in  the  header.  All  three  of  these 
models  are  quite  similar  except  for  the  respective  numerical  algorithms 
employed.  In  all  three  cases  the  two-dimensional  continuity  equation 
Is  written  as: 


T 
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header  density 
header  neat  capacity 

longitudinal  header  thermal  conductivity 
transverse  header  thermal  conductivity 

Figure  5 shows  an  electrical  analog  for  these  header  thermal  models. 

Note  that  as  with  the  previous  thermal  model  'DHTEMP' , there  is  no 
power  dissipation  within  the  header  yielding  an  energy  continuity 
equation  of  reduced  complexity.  This  equation  requires  two  boundary 
conditions  on  temperature  for  each  of  the  two  axes  of  thermal  con- 
duction. The  boundary  conditions  on  thermal  conduction  perpendicular 
to  the  diode  axis  are  the  same  as  for  DHTEMP,  i.e.,  the  respective 
diode  temperatures  and  the  ambient  temperature.  The  temperature 
boundary  values  for  thermal  conduction  parallel  to  the  diode  axis  are 
always  the  same  as  those  specified  for  thermal  conduction  along  the 
diode  axis,  i.e.,  ambient  temperature  by  default  and  blocking  boundary 
conditions  when  specified  through  the  appropriate  simulation  control 
parameter. 

All  three  models  incorporate  an  iterative  solution  procedure 
which  is  controlled  by  a specified  maximum  RMS  change  in  temper  iture 
between  successive  iterations  on  a maximum  number  of  iterations.  If 
the  maximum  number  of  iterations  is  achieved  without  obtaining  the 
specified  convergence,  a convergence  failure  message  is  generated. 

The  first  two-dimensional  header  thermal  model  formulated  is 
'DHT2D'.  This  model  is  implemented  using  a vertical  line  procedure 


where : 

PH  “ 
CH 

’Six  - 

V - 


and  incorporates  the  latest  temperature  values  on  an  iteration  basis. 


T(N-1,M)  A T(N,M)  A T (N+1,M) 
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The  vertical  line  approach  requires  that  the  partial  derivative  of 
temperature  along  tu«t  vertical  axis  (axis  perpendicular  to  the  diode 
axis)  be  formulated  by  an  implicit  finite  difference  scheme  and  that 
the  horizontal  partial  derivative  (axis  parallel  to  the  diode  axis) 
be  formulated  by  an  explicit  finite  difference  scheme.  The  restriction 
placed  on  the  usage  of  the  latest  temperature  values  relates  to  the 
implicit  formulation  and  means  that  although  some  new  values  of 
temperature  do  become  available  during  each  iteration  through  the 
temperature  field,  these  values  are  not  to  be  used  in  evaluating  the 
remaining  temperatures  for  the  same  iteration.  Rather,  the  temperature 
values  evaluated  during  a given  iteration  are  used  by  the  succeeding 
iteration  only.  Application  of  this  algorithm  to  the  two-dimensional 
continuity  equation  yields  the  following  system  of  equations: 

A T^J  ( N , M+l ) + B T^J  (N,M)  + B T^J  (N.M-l)  = (8) 


where : 


time  step  number 


iteration  number 


node  number  along  the  diode  axis 
- node  number  perpendicular  to  the  diode  axis 
header  temperature  field 
node  spacing  along  x-axis 


node  spacing  along  y-axis 
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and, 

A = -Kny  At  Ax2 

B = 2 ^ At  Ax?+  Ay2Ax2 

C = At  Ax2 

°N,M  = Ax2Ay2TJ(N»M>  + ^ Ay2At  [Tg+1(N-1,M)  + Tg+1(N-1,M) 

- 2 Tg+1(N,M) ] 

Note  that  solving  the  above  system  of  equations  yields  a single 
temperature  profile  through  the  header.  Accordingly,  one  iteration 
of  this  algorithm  requires  that  the  above  system  of  equations  be 
solved  for  each  node  point  along  the  diode  axis  except  for  the 
respective  end  nodes  which  are  accommodated  through  the  boundary 
conditions  applied  at  these  points.  Thermal  coupling  between  these 
linear  profiles  is  provided  through  successive  iterations. 

The  second  two-dimensional  header  thermal  model  is  developed 
to  accelerate  convergence  and  improve  stability  of  the  above  algorithm. 
The  only  change  is  to  use  the  latest  temperature  values  as  they 
become  available  during  an  iteration  rather  than  delaying  until  the 
next  iteration.  Only  the  'D'  coefficient  of  the  previous  formulation 
is  changed  by  this  innovation  and  it  becomes: 

°N,M  = Ax2  Ay2rJ(N.M>  + ^ Ay2  t [Tg+1(N+1,M) 
t Tg^(N-l.M)  - 2 Tg+1(N,M)  ] 


(9) 
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Note  that  the  Tg+^(N-1,M)  temperature  represents  the  latest  approxi- 
mation for  temperature  at  that  node  rather  than  the  temperature 
value  evaluated  during  the  previous  iteration.  It  is  shown  in  the 
simulation  results  section  that  this  change  does  increase  both  stability 
and  the  rate  of  convergence. 

The  third  two-dimensional  header  thermal  model  is  formulated  like 
the  previous  model  except  that  a horizontal  line  technique  is  employed 
rather  than  the  previously  used  vertical  line  approach.  In  this  case, 
tne  horizontal  partial  derivative  is  formulated  implicitly;  whereas, 
the  vertical  partial  derivative  is  formulated  explicitly.  Applications 
of  the  line  technique  in  analysis  of  two-dimensional  semiconductor 
impurity  diffusions  [8]  suggested  that  the  horizontal  line  approach 
would  be  more  effective  for  the  header  thermal  model  than  the  vertical 
line  technique.  The  convergence  test  presented  later  supports  this 
conclusion.  As  with  the  previous  model,  DAT2D1,  this  model,  named 
DH2D2,  incorporates  the  new  temperature  values  as  they  become  available 
rather  than  on  an  iteration  basis.  The  subsequent  implementation  of 
DHT2D2  yields  the  following  system  of  equations: 

A T(N+1,M)^J  + B T(H,M)£j;  + B T(N-l,M)j£J  = D (10) 

where: 

A ■ - Sx  *c  iy2 

B = [2  Kjjx  At  Ay2  + pH  cH  Ax2Ay2  ] 

c « - ^ At  Ay2 

°N,M  “ PHCH  Ax2  Ay2T<N’M>J  + % At  ax2[T(N,M+1)^+1 
+ TCN.M-lfg1  - 2 T(N,M)g+1] 
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Whereas  the  previous  two  two-dimensional  header  thermal  models  resulted 
in  a system  of  linear  equations  for  temperature  profiles  perpendicular 
to  the  diode  axis,  this  model  evaluates  temperature  profiles  parallel 
to  the  diode  axis.  That  is,  the  above  system  of  equations  must  be 
solved  once  for  each  row  of  nodes  in  the  header,  except  for  the 
bottom  row  which  is  maintained  at  ambient  temperature.  This  procedure 
constitutes  one  iteration.  Thermal  coupling  is  maintained  between 
the  row-wise  temperature  profiles  through  the  iterative  solution 
procedure. 

Different  header  thermal  conductivities  were  defined  for  the 
header  vertical  and  horizontal  axes  such  that  the  new  two-dimensional 
header  thermal  models  could  conveniently  simulate  the  earlier  quasi- 
two-dimensional  model  by  defining  = 0.0.  This  provided  a means 
of  testing  the  two-dimensional  models  during  their  early  development. 

In  all,  four  header  thermal  models  have  been  developed,  and  are 
named:  DHTEMP,  DHT2D,  DHT2D1 , and  DHT2D2 . DHTEMP  is  a quasi-two- 
dimensional  model  and  DHT20,  DUT2D1  and  DHT2D2  are  two-dimensional 
models.  The  essential  difference  between  the  three  two-dimensional 
models  involves  the  numerical  techniques  used.  A comparison  of 
convergence  characterist ics  and  TSB  transient  simulations  using  the 
above  header  thermal  models  is  made  in  the  simulation  results  section. 
Figure  6 shows  an  electrical  analog  for  the  total  diode  thermal  model 
which  includes  both  the  diode  and  two-dimensional  header  thermal 
models  connected  together. 


T(N-1,M)/1  T(N,M)  A T(N+1,M) 


Fig.  6 . Electrical  Analog  for  the  Two-Dimensional  Diode  Thermal  Model.  (T  - temperature,  C-specific 
heat,  R-  thermal  resistance,  H-  header,  D- diode,  Y-v  axis,  X-x  axis,  6 - heat  generation, 
N-  x axis  node  number,  M-y  axis  node  number). 
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I 

DIODE  ELECTRICAL  MODEL 

The  one-dimensional  diode  electrical  model  is  time  dependent 
and  is  designed  to  maintain  the  condition  of  avalanche  breakdown 
throughout  the  TSB  transient  simulations.  As  a consequence  of  the 
temperature  dependence  of  the  avalanche  ionization  coefficients,  hole 
and  electron  mobilities,  and  intrinsic  carrier  concentrations,  the 
diode  electrical  model  is  closely  coupled  to  the  thermal  transient 
which  accompanies  TSB. 

The  maximum  depletion  region  electric  field  is  assumed  to  occur 
at  the  diode  metallurgical  junction  and  the  depletion  region  electric 
field  profile  is  evaluated  analytically  from  this  value.  The  diode 
impurity  profile  is  assumed  to  be  an  epitaxial  configuration  as  shown 
in  Figure  7.  This  configuration  readily  reduces  to  the  classical 
two-sided  abrupt  junction  by  assigning  equal  values  to  the  epitaxial 
and  upper  region  (background)  doping  levels.  Space  charge  effects 
in  the  depletion  region  are  taken  into  consideration  by  using  effective 
space  charge  densities  in  evaluating  the  depletion  region  electric 
field  profile.  The  effective  space  charge  densities  are  defined  as: 


PM  = 

J 

V 

(ID 

sat 

PL  = 

q NL  - PM 

(12) 

pEp 

q NEp  - PM 

(13) 

PU  = 

q NU  " PM 

(U) 

.J 
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The  diode  is  assumed  to  be  in  avalanche  breakdown  when  the  avalancne 

breakdown  integral  [9]  is  equal  to  one,  plus  or  minus  a specified  error. 

W 

J * dx  =*  1 ± y (15) 

o 

where: 

W - depletion  region  width 

* - avalanche  ionization  coefficient  corresponding  to  the 

lower  region  impurity  type  (n  or  p) . 

Y - specified  error  for  avalanche  breakdown  integral 

The  avalanche  ionization  coefficient  in  the  above  integral  is  chosen 
to  correspond  to  the  lower  impurity  region  type,  assuming  that  this 
region  represents  the  high  doped  side  of  the  junction.  Under  these 
conditions,  the  carriers  injected  from  the  epitaxial  or  background 
regions  will  dominate  the  avalanche  breakdown  mechanism  and  will  cor- 
respond to  the  majority  carrier  type  for  the  lower  diode  region.  The 
avalanche  ionization  coefficients  are  evaluated  by  the  following 
relations: 

-b  /E 

“ - a II  - 0 (T-To) ] e P (16) 

-b  /E 

«n  - an  [1  - 6(T-To)]  e n (17) 

where : 

a ” 3.8  x 10^  c 

P ' 

b - 1.74  x 106  V c -1 

p m 


24 


a = 2.25  x 10  c -1 

n m 


b = 3.26  x 106  v c _1 

n m 


The  above  description  of  the  ionization  coefficients  was  derived 
from  Sze's  [10]  characterization  of  these  coefficients  and  represents 
an  empirical  approximation  of  their  thermal  dependence.  The  avalanche 
ionization  thermal  coefficient  6 has  units  of  °R  ^ and  specifies  a 

-3 

linear  dependence  on  temperature.  B is  assigned  a value  of  2.5  x 10 
and  may  be  changed  slightly  to  enhance  the  thermal  dependence  of  the 
avalanche  ionization  coefficients.  B was  named  ALFATD  for  the  computer 
program  version  of  the  diode  model.  It  should  be  emphasized  that 
although  this  formulation  for  the  thermal  dependence  of  the  ionization 
coefficients  is  not  completely  accurate,  it  does  exhibit  a reasonably 
good  description  of  the  thermal  behavior  of  the  coefficients. 

Iterations  on  the  maximum  electric  field  are  controlled  by  an 
interval  halving  routine.  Lower  and  upper  limits  for  the  maximum 
electric  field  are  established  through  simulation  parameters  and  are 
reinitialized  at  the  beginning  of  each  series  of  iterations.  The 
average  of  these  two  limits  is  used  as  a trial  value  for  the  maximum 
electric  field.  Next,  the  electric  field  profile  and  the  avalanche 
breakdown  integral  are  evaluated.  If  the  integral  yields  a value 
greater  than  one,  the  trial  value  for  the  maximum  electric  field 
is  too  large  and  the  trial  value  replaces  the  upper  limit.  If  the 
integral  is  less  than  one,  the  trial  value  is  too  small  and  the  lower 
limit  is  replaced  by  the  trial  value.  This  procedure  is  repeated 


until  the  avalanche  breakdown  integral  takes  on  the  specified  value 


T 
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or  a maximum  number  of  iterations  is  performed.  If  the  latter  case 
occurs  a convergence  failure  message  is  generated. 

Once  the  depletion  region  electric  field  has  been  determined, 
the  bulk  region  electric  fields  are  evaluated.  This  computation  is 
performed  assuming  that  the  currents  with  in  these  regions  consist 
entirely  of  drift  components  and  that  the  carrier  concentrations 
correspond  to  the  thermal  equilibrium  values. 


J = q (yn  n + Up  P ) E 


(18) 


J - diode  current  density 
q - unit  charge 
Up  - hole  mobility 
U - electron  mobility 

n - thermal  equilibrium  electron  concentration 
p - thermal  equilibrium  hole  concentration 
E - electric  field 


Since  the  carrier  mobilities  are  stongly  dependent  on  impurity 
concentration,  electric  field,  and  temperature,  it  is  necessary  that 
the  mobility  coefficients  be  formulated  in  terms  of  these  quantities. 
Accordingly,  the  mobilities  are  evaluated  using  the  following 
equation  [1,  11,  12]: 


u 


1 + 


N 


+ 


05) 
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r 


where  for  silicon: 


y 

a 

b 

c 

d 

e 

Holes 

480 

2.5 

4xl016 

81 

6. lxlO3 

1.6 

2.5x10* 

Electrons 

1400 

2.5 

3xl016 

350 

3.5xl03 

8.8 

7.4xl03 

u - mobility 
T - temperature 
N - impurity  concentration 
E - electric  field 

y,  a,  b,  c,  d,  e - parameters 

The  mobility  dependence  on  electric  field  prevents  an  analytic 
solution  of  equation  (16)  for  the  electric  field  as  a function  of  cur- 
rent density,  temperature,  and  impurity  concentration.  Rather  an 
algorithm  based  on  the  Newton -Raphson  technique  is  used  to  iteratively 
evaluate  the  electric  field  at  each  node  point  along  the  diode  axis. 

Once  the  depletion  region  and  bulk  region  electric  field  profiles 
have  been  evaluated,  they  are  combined  to  yield  the  diode  electric 
field  profile. 

A flow  chart  for  the  diode  electrical  model  is  shown  in  Figure 
8.  The  relationship  between  the  electrical  model  and <.he  rest  of 
the  program  is  best  demonstrated  in  the  system  flow  chart  presented 
in  Figure  2.  In  this  figure,  the  electrical  model  is  represented  by 
the  'Evaluate  the  Electric  Field'  block. 
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PROGRAM  DESCRIPTION 

The  numerical  diode  model  developed  in  the  previous  section  is 
programmed  in  the  Univac  FORTRAN  V computer  language.  This  program  is 
designed  to  simulate  TSB  transients  up  to  the  onset  of  the  TSB  transition. 
The  program  is  very  versatile  and  allows  convient  specification  of  diode 
design  and  simulation  conditions.  This  section  defines  the  diode  and 
simulation  parameters  and  provides  both  instructions  and  suggestions  for 
executing  the  diode  model. 

The  diode  model  consists  of  a main  program  and  several  subroutines. 
Listings  for  these  programs  are  presented  in  Appendix  B.  Both  the  main 
program  and  the  subroutines  are  well  documented  with  comment  statements. 

The  TSB  transient  simulations  performed  by  the  diode  model  are 
defined  through  two  groups  of  parameters.  The  first  group  controls  the 
simulation  and  is  named  SPARM.  The  second  group  specifies  the  diode 
design  and  is  named  DPARM.  Default  values  and  definitions  for  these  two 
parameter  sets  are  presented  in  Tables  1 and  2. 

As  many  parameters  as  desirable  may  be  specified  for  a simulation 
with  all  unspecified  parameters  taking  on  their  respective  default  values. 
Parameter  specifications  must  be  made  on  a group  basis  with  appropriate 
starting  and  terminating  sentinels  for  each  group.  An  example  simulation 
parameter  set,  with  each  line  representing  an  input  line  is  shown  in 
Figure  9.  The  resulting  simulation  printout  for  this  example  is  pre- 
sented in  Appendix  A.  The  parameter  groups  must  appear  in  the  order 
and  format  shown.  As  many  simulations  as  required  may  be  specified  in  a 
single  run  by  adding  additional  sets  of  parameter  specifications.  When 


« 
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$SPARM 
CUR  = 7.5E3 

First  Parameter  NNH  = 12 


Set  Group  Name 


First  Parameter 
Set  Terminator 


Second  Parameter 
Set  Group  Name 

Second  Parameter 


DTIME  = 5.0E-9 
LTRNtN  = 1 
LTEMP  = 1 
LDHTEM  = 1 
PLOT  = 1 
PVOLTX  = 500 
IDHTMO  = 4 
SMAX  = 20 
PLTIME (1)  = 100.0 
$END 
$DPARM 

THKDPF(31)  = 0.5 
$END 


Set  Terminator 


S SPARM  parameter  set 


DPARM  parameter  set 


Firgure  9.  SPARM  and  DPARM  Parameter  Sets  for  a Thermal  Second 

Breakdown  Demonstration  Simulation  (simulation  output 
listing  presented  in  Appendix  A). 
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TABLE  1 

SIMULATION  CONTROL  PARAMETER  SET  (SPARM) 


Parameter  Default  Value 

1.  CUR  = 1.25  x 104  amps/cn2 

2.  NND  = 101 

3.  NNH  = 0 

4.  IBND  = 0 

5.  TIMEMX  = 1.0  sec 

6.  TMAX  = 700. 0°K 

7.  TPMAX  = 800. 0°K 

8.  ITSNMX  = 200 

9.  DTIME  = 1.0  x 10-9  sec 

10.  LTRNIN  = 2 

11.  LTEMP  = 0 

12.  LDHTEM  = 0 

13.  PLOT  = 1 

14.  PVOLTX  = 250.0  volts 

15.  IGRID  = 0 


Parameter  Definition 


diode  current  density 

number  of  diode  node  points  along 
diode  axis,  < 101 

number  of  node  points  in  diode 
header,  maximum  of  twelve,  and 
zero  assumed  for  NNH  < 4 

temperature  boundary  condition 
sentinel,  = 0 - constant  temperature 
boundary  conditions,  > 1 - blocking 
boundary  conditions 

maximum  simulation  time 

maximum  simulation  temperature 
(TSB  temperature) 

maximum  value  for  temperature  plots 

maximum  number  of  time  steps 

time  step  increment 

transient  data  print  increment 

= 1 - diode  temperature  pro- 
files are  printed  at  times  specified 
in  parameter  array  PLTIME 

= 1 - temperature 

listings  are  printed  at  times 

specified  in  parameter  array  PLTIME 

= 1 - temperature,  impurity,  and 
electric  field  profiles  and  voltage 
transients  are  plotted  at  the 

times  specified  in  parameter  array 
PLTIME 

maximum  value  for  voltage  plots 

grid  sentinel  for  plots,  = 1 - grid 
generated 
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TABLE  1 (continued) 

SIMULATION  CONTROL  PARAMETER  SET  (SPARM) 


Parameter  Default  Value 

16.  IRMSMI  = 1.0  x 103 

17.  ITERMX  = 10 

18.  ITLST  = 10000 

19.  IDHIMO  = 1 

20.  SMAX  = 10 

21.  DTHMAX  = 1.0  x 10_A 

22.  ITPRH  = 0 

23.  PLTIME(l)  = 100.0  sec 

• » * 

• • • 

• • • 

32.  PLTIME(IO)  = 100.0  sec 

33.  EMAXL  = 1.0  x 105  volts/cm 

34.  EMAXU  = 1.0  x 106  volts/cm 


Parameter  Definition 

maximum  RMS  temperature  change 
allowed  between  successive  iter- 
ations for  diode  axis  temperature 
profile 

maximum  number  of  iterations 
between  time  steps  for  diode  axis 
temperature  profile 

iteration  diagnostic  print  interval 
for  diode  axis  temperature  profile 

header  thermal  model  sentinel, 

1 - DHTEMP,  quasi-two-dimensional 

2 - DHT2D,  two-dimensional  transverse 

3 - DHT2D1,  two-dimensional  transverse 

4 - DHT2D2,  two-dimensional  logitu- 

dinal 

maximum  number  of  iterations  between 
time  steps  for  header  thermal  model 

maximum  RMS  temperature  change 
allowed  between  successive  iterations 
for  header  temperature  model 

iteration  diagnostic  print  interval 
for  the  header  thermal  model 

diode  temperature  profile  plot  and/or 
list  times  (diode  temperature  profiles 
are  available  at  simulation  termi- 
nation by  default),  timers  must  be 
stored  in  array  PLTIME  in  Chronolo- 
gical order 


lower  bound  for  electric  field  at 
metallurgical  junction  (maximum 
junction  electric  field) 

upper  bound  for  electric  field 
at  metallurgical  junction  (maximum 
junction  electric  field) 


J 


TABLE  1 (continued) 
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SIMULATION  CONTROL  PARAMETER  SET  (SPARM) 


Parameter  Default  Value 
35.  EMAXB  = 1.0  x 10^  volts/cm 


36.  EINT  = 100.0  volts/< 


37.  UPE  = 1 


38.  AERMAX  = 0 1 


39.  EERMAX  = 1.0  x 10' 


40.  ITCMAX  = 10 


41.  IDB0 

42.  IDB1 


43.  IDB2 


44.  IDB3 


45.  IDB4 


46.  PDTP 


47.  PDIP  = 0 


48.  PDEP  = 0 


49.  PDVP  = 0 


Parameter  Definition 

maximum  electric  field  value  for 
bulk  regions 

starting  value  of  electric  field 
for  determining  bulk  region  electric 
field  profiles 

time  dependent  electric  field  sentinel 

0 - time  independent  electric  field 

1 - time  dependent  electric  field 

maximum  error  allowed  in  avalanche 
breakdown  integral 

maximum  RMS  change  in  bulk  region 
electric  field  allowed  between 
successive  iterations 

maximum  number  of  iterations  for 
header  temperature 

diagnostic  print  sentinel  for  subroutines: 
Not  Used 

INAFLA 

BKDEPL 

EPROF 

DOPLG 

diode  temperature  profile  plot 
sentinel,  1-profiles  at  times 
specified  in  parameter  array 

PLTIME 

impurity  profile  plot  sentinel. 

1 - plot  generated 

electric  field  plot  sentinel, 

1 - plot  initial  and  final  electric 
field  profiles 

diode  transient  voltage  plot  sentinel, 

1 - plot  final  diode  voltage  pro- 
files. (A  - lower  bulk  region 
voltage,  B - depletion  region  voltage, 

C - upper  bulk  region  voltage, 

D - total  diode  voltage) 
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TABLE  2 

DIODE  DESIGN  AND  PHYSICAL  PROPERTY 
PARAMETER  SET  (DPARM) 


1. 

Parameter  Default  Value 

THKD  = 1.0  watts/cm-°k 

2. 

THKHX 

= 

0.46  watts/cm-°k 

3. 

THKHY 

= 

0.46  watts/cm-°k 

4. 

DDEN 

= 

2.3  gm/cm^ 

5. 

HDEN 

= 

2 

4.0  gm/cm 

6. 

DSPEC 

= 

0.7  J/gm-°k 

7. 

HSPEC 

= 

0.79  J/gm-°k 

8. 

XDT 

= 

1.0  x 10  ^ cm 

9. 

XDH 

= 

2.5  x 10  ^ cm 

10. 

XLDEP 

= 

0.0  cm 

11. 

XMET 

= 

2.0  x 10” 3 cm 

12. 

XEPDEP 

= 

4.0  x 10  3 cm 

13. 

XUDEP 

= 

0.0  cm 

14. 

XL 

4.0  x 10  ^ cm 

15. 

DOPL 

1.0  x 1017  cm-3 

16. 

DOPEP 

» 

1.0  x 1016  cm-3 

17. 

DOPU 

1.0  x 10'6  cm'3 

18. 

VEL 

= 

1.0  x 107  cm/sec 

19. 

ALFARD 

- 

2.4  x 10~3  °k-1 

20. 

NP1 

s 

1 

Parameter  Definition 

diode  thermal  conductivity 
(silicon) 

header  longitudinal  thermal 
conductivity  (sapphire) 

header  transverse  thermal 
conductivity  (sapphire) 

diode  density  (silicon) 

header  density  (sapphire) 

diode  heat  capacity  (silicon) 

header  heat  capacity  (sapphire) 

diode  thickness 

header  thickness 

lower  depletion  region  boundary 
(not  used) 

metallurgical  junction  location 

epitaxial  layer  boundary 

upper  depletion  region  boundary 
(not  used) 

diode  length 

lower  bulk  region  doping  concentration 

epitaxial  region  doping  concentration 

upper  bulk  region  doping  concentration 

mobile  carrier  saturation  velocity 
(silicon) 

avalanche  ionization  coefficient 
thermal  dependence  factor 

diode  orientation,  0 - PN,  1 - NP 


MMH*1 


TABLE  2 (continued) 
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DIODE  DESIGN  AND  PHYSICAL  PROPERTY 
PARAMETER  SET  (DPARM) 


Parameter  Default  Value 


Parameter  Definition 


21.  THKDEP(l) 


diode  thermal  conductivity 
perturbation  factor 


121.  THKDPF(lOl)  = 1.0 


122.  THKHPF(l) 


diode-to-header  thermal  conductivity 
perturbation  factor 


222.  THKHPF(lOl)  = 1.0 


A 


f 
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multiple  simulations  are  requested,  only  parameter  changes  between 
successive  simulations  should  be  specified.  The  default  values  are  used 
for  the  first  simulation  only.  For  succeeding  simulations,  the  previous 
simulation  parameter  values  are  treated  as  default  values. 

The  program  requires  36 K words  of  memory  for  execution.  Typical 
run  times  on  a Univac  1107  computer  are  from  one-to-two  minutes. 

Several  guide  lines  are  listed  to  facilitate  program  execution: 

1.  If  the  number  of  node  points  in  the  header  (NNH)  is  less 
than  four,  quasi-two-dimensional  thermal  conduction  without 
intermediate  header  node  points  is  assumed. 

2.  A simulation  is  terminated  when  one  of  the  following  con- 
ditions is  achieved: 

A.  Simulation  time  > TIMEMX 

B.  Time  step  count  > ITSNMX 

C.  Maximum  diode  temperature  > TMAX 

3.  The  transient  data  print  increment  (LTRNIN)  must  be  chosen 
such  that  no  more  than  500  transient  data  lines  are  printed. 

4.  Two-dimensional  temperature  listings  for  the  diode-header 
combination  is  generated  at  the  times  specified  in  the  PLTIME 
array  and  at  program  termination  by  default,  if  LDHTEM  = 1. 

5.  The  diode  axis  temperature  profile  is  stored  at  the  times 
specified  in  the  PLTIME  array  and  at  program  termination  by 
default.  A listing  of  these  profiles  is  generated  if 
LTEMP  = 1 and  a plot  is  generated  if  IPLOT  = 1. 

6.  The  time  step  increment  should  be  specified  to  yield  a 
one  hundred  to  two  hundred  time  step  simulation. 

7.  For  most  simulations,  the  diode  maximum  temperature  transient 
plot  and  listing  are  adequate.  Both  of  these  outputs  are 
generated  by  default. 

8.  The  thermal  conductivity  perturbation  factors  are  listed 
only  when  one  or  more  of  these  factors  are  assigned  a value 
other  than  one. 


36 


SIMULATION  RESULTS 

The  results  presented  in  this  section  are  intended  to  characterize 
the  diode  model.  Several  different  type  simulations  are  presented 
to  demonstrate  the  capabilities  and  flexibility  of  the  model.  The 
convergence  characteristics  of  the  four  header  thermal  models  and 
the  sensitivity  of  the  diode  to  various  model  parameters  is  also 
presented. 

Ten  different  series  of  simulations  involving  in  excess  of 
one  hundred  simulations  were  performed  to  accumulate  the  data 
presented  in  this  section.  A simulation  series  is  considered  to  be 
a series  of  closely  related  simulations  which  involve  systematic 
changes  in  a parameter  or  parameters.  A summary  of  these  simulation 
series  is  presented  in  Table  3.  This  table  shows  the  changes  made 
in  the  default  parameter  values  in  order  to  perform  the  respective 
simulations.  The  parameters  that  were  varied  in  a particular 
simulation  series  are  indicated  by  PAR.  The  different  series  are 
numbered  in  the  table  and  are  identified  on  the  various  graphs  pre- 
sented in  this  section  by  the  respective  simulation  number  preceded 
by  SS,  e.g.,  SS-5.  There  are  several  common  parameter  values  among 
the  simulations  presented.  These  are  summarized  here  for  convenience: 


NND  = 

101 

number  of  diode  node  points 

NNH  = 

0 or  12 

number  of  header  node  points 

TMAX  = 

700°K 

TSB  temperature 

THRU  = 

1.0  watts/cm  - °K 

diode  thermal  conductivity 
(silicon) 

THKHX  = 

0.46  watts/cm  - °K 

longitudinal  header  thermal 
conductivity  (sapphire) 
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THKHY 

0.46  watts/cm  = 1 

DDEN 

= 

2.3  gm/cm^ 

HDEN 

= 

3 

4.0  gm/cm 

DSPEC 

s 

0.7  J/gm  - °K 

HSPEC 

= 

0.79  J/gm  - °K 

XDT 

= 

1.0  pm 

XDH 

= 

250  or  4.8  pm 

XMET 

= 

20  pm 

XL 

DOPL 

= 

40  pm 

17  -3 

1.0  x 10  cm 

DOPEP 

= 

DOPU 

DOPU 

= 

. _15  ,.17 

>_  10  , _<  10  cm 

VEL 

= 

1.0  x 10^  cm/sec 

NP1 

= 

1 

transverse  header  thermal 
conductivity  (sapphire) 

diode  density  (silicon) 

header  density  (sapphire) 

diode  heat  capacity  (silicon) 

header  heat  capacity  (sapphire) 

semiconductor  film  thickness 

header  thickness 

location  of  metallurgical 
junction 

diode  length 

lower  region  doping  concentration 

epitaxial  region  doping 
concentration 

upper  region  doping  concentration 
carrier  saturation  velocity 
np  diode  orientation 


Simulation  SS-1  was  a demonstration  simulation  and  the  output 
listing  for  this  simulation  is  presented  in  appendix  - A. 

A comparison  between  TSB  delay  times  for  the  quasi-two-dimen- 
sional  header  thermal  model  with  zero  and  twelve  header  node  points 
is  shown  in  Fig.  10.  For  short  delay  times,  these  models  yield 
similar  results;  whereas,  for  long  delay  times  the  results  differ 
by  two  orders-of-magnitude.  This  behavior  is  expected  since  increasing 
the  number  of  node  points  in  the  header  tends  to  increase  the  rate 
of  heat  transfer  out  of  the  diode.  Of  course,  this  effect  saturates 
quickly  as  the  accuracy  of  the  transverse  header  temperature  derivative 
increases.  The  unexpected  result  is  the  inflection  point  exhibited 
by  the  more  comprehensive  quasi-two-dimensional  model  (NNH=12).  This 
anomaly  is  suspected  to  be  associated  with  the  quasi-two-dimensional 
thermal  model  concept  [13]. 
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Fig.  11,  shows  a comparison  of  TSB  delay  times  between  constant 
and  blocking  temperature  boundary  conditions.  Both  cases  are  for 
comprehensive  quasi-two-dimensional  thermal  conduction  (NNH=12)  and 
both  curves  exhibit  the  previously  noted  inflection  point.  As 
expected,  the  boundary  conditions  on  temperature  at  the  diode  contacts 
have  virtually  no  effect  on  the  delay  time  for  high  currents.  For 
low  currents  the  constant  temperature  boundary  conditions  exhibit 
slightly  longer  delay  times  corresponding  to  a greater  heat  loss  at 
the  contacts. 

Fig.  12,  presents  a comparison  of  TSB  delay  times  for  zero  and 
non  zero  bulk  region  electric  field  values.  Both  curves  represent 
comprehensive  quasi- two-dimensional  thermal  conduction  (NNH=12)  and 
exhibit  the  previously  observed  inflextion  points.  The  bulk  region 
electric  fields  are  seen  to  make  a significant  contribution  to  the 
TSB  delay  times.  The  difference  is  essentially  one  order-of-magnitude 
at  low  current  values. 

A comparison  of  TSB  delay  times  for  quasi  and  full  two-dimensional 
header  thermal  models  is  shown  in  Fig.  1J.  The  two  models  agree  well 
for  high  currents.  However,  they  differ  greatly  for  low  currents 
as  expected.  Further,  the  two-dimensional  model  does  not  exhibit  the 
inflection  point  observed  for  all  the  comprehensive  quasi-two-dimensional 
simulations.  It  appears  that  the  inflection  point  is  associated 
with  the  quasi-two-dimensional  thermal  model  concept  as  predicted 
earlier.  The  full  two-dimensional  model  exhibits  a very  strong 
delay  time  dependence  on  current  for  low  current  values. 

Fig.  14,  presents  maximum  junction  temperature  transients  for 
several  different  driving  currents  and  the  two-dimensional  header 
thermal  model.  These  curves  emphasize  the  effects  of  the  avalanche 
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ionization  coefficient  thermal  dependence.  For  low  current  values, 
the  temperature  transient  almost  saturates  before  reaching  the  TSB 
temperature.  However,  temperatures  sufficiently  high  to  cause  an 
appreciable  decrease  in  the  avalanche  ionization  coefficients  are 
finally  achieved  after  long  delay  times.  Once  the  avalanche  ion- 
ization coefficients  begin  to  decrease,  the  junction  electric  field 
is  forced  to  increase  rapidly  to  maintain  avalanche  breakdown.  This 
increase  in  electric  field  causes  a corresponding  increase  in  power 
dissipation  with  in  the  junction.  A thermal  runaway  condition  occurs. 
These  effects  are  manifest  in  the  diode  temperature  transient  as 
the  TSB  temperature  is  approached.  These  effects  also  occur  at  high 
currents,  but  are  masked  out  by  the  rapid  increase  in  diode  temper- 
ature associated  with  the  high  current  level. 

Delay  time  versus  doping  concentration  characteristics  for  the 
reduced  quasi-two-dimensional  header  thermal  model  (NNH=0)  are  shown 
in  Fig.  15.  The  delay  time  is  shown  to  increase  as  the  low  doping 
concentration  is  increased.  This  corresponds  to  a decreasing  depletion 
region  width  and  a decreasing  resistivity  for  the  low  doping  con- 
centration, both  of  which  tend  to  increase  the  delay  time.  For 
comparison,  one  point  for  each  of  the  two  driving  current  values  is 
also  shown  for  the  full  two-dimensional  thermal  model.  Note  that  the 
two-dimensional  model  increases  the  delay  time  only  slightly  at  the 
high  current  value;  whereas,  the  delay  time  becomes  infinite  for  the 
low  current  value.  The  infinite  delay  time  implies  that  the  diode 
achieves  a thermal  steady  state  without  reaching  the  TSB  temperature. 
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Delay  Time  Versus  F.pitaxial  and  Upper  Region  Doping 
Concentration  for  two  Different  Driving  Current  Values. 
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Fig.  16  characterizes  the  dependence  of  TSB  transients  on 

the  specified  avalanche  breakdown  integral  error,  AERMAX.  The 

maximum  junction  temperature  at  1.8  ps  into  the  simulation  is 

plotted  as  a function  of  this  error.  The  simulations  were  made 

with  the  reduced  quasi-two-dimensional  model  and  with  a driving 
4 2 

current  of  10  amps/cm  . The  results  indicate  that  an  avalanche 
breakdown  integral  error  less  than  or  equal  to  two  tenths  is  most 
desirable.  The  trade  off  for  further  decreases  is  a corresponding 
increase  in  execution  time. 

The  effects  of  the  avalanche  ionization  temperature  coefficient 

(6  or  ALFATD)  on  the  TS1  transients  are  shown  in  Figs.  17  - 19. 

Fig.  17,  shows  the  TSB  delay  time  as  a function  of  6 for  the  reduced 

3 

quasi-two-dimensional  model  and  a driving  current  of  7.0  x 10 
2 

amps/cm  . As  expected,  increases  in  6 result  in  a decreased  delay 
time.  Further  refinement  of  the  value  for  6 may  be  warranted. 

Figs.  18  and  19,  show  maximum  junction  temperature  transients  for 
several  of  the  data  points  in  Fig.  17.  The  critical  temperature 
for  the  avalanche  ionization  coefficients  is  set  equal  to  700  °K 
for  the  diode  model.  This  is  the  temperature  at  which  the  avalanche 
ionization  coefficients  take  on  a value  of  zero  as  a consequence  of 
their  thermal  dependence.  8,  the  avalanche  ionization  temperature 
coefficient,  simply  determines  the  sensitivity  of  the  ionization 
coefficients  to  changes  in  temperature  as  the  critical  temperature 
is  approached.  The  formulation  for  the  ionization  coefficients  is 
shown  in  Equation  fb . Although  the  TSB  temperature  and  the  avalanche 
ionization  critical  temperature  have  both  been  assigned  a value  of 
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700  °K  for  the  simulations  presented,  further  refinement  may  show 
that  one  or  both  of  these  values  require  adjustment.  Although  the 
avalanche  ionization  coefficient  critical  temperature  breakpoints 
are  not  obvious  for  the  logorithmic  time  scale  in  Fig.  18,  they  can 
be  clearly  seen  on  the  linear  time  scale  of  Fig.  19. 

Figs.  20  _ 26j  show  the  convergence  characteristics  for  the 
three  two-dimensional  header  thermal  models.  All  the  curves  shown 

represent  iterations  for  the  first  time  step  of  a simulation  with 

3 2 

a driving  current  of  8.0  x 10  amps/cm  . Fig.  20,  shows  the  con- 
vergence behavior  for  DHT2D  for  eight  different  time  step  sizes. 

This  model  uses  the  vertical  line  technique  and  fails  to  converge 
for  time  steps  greater  than  3 nsec.  Further,  the  convergence  rate 
is  very  sensitive  to  the  time  step  size.  Fig.  21,  compares  the 
convergence  behavior  of  DHT2D  for  two  different  header  thicknesses. 

The  graph  shows  that  the  model  convergences  more  rapidly  for  the 
larger  transverse  header  node  spacing.  On  the  other  hand,  the 
stability  behavior  is  very  similar  for  the  two  cases. 

Figs.  22  and  show  the  convergence  behavior  for  DHT2D1  which 
features  the  vertical  line  formulation  and  incorporates  the  very 
latest  values  of  temperature  during  the  iteration  procedure. 

A comparison  between  the  two  vertical  line  models,  DHT2D  and 
DHT2D1,  for  two  different  header  thicknesses  is  shown  in  Fig.  24 . 

In  both  cases  the  formulation  incorporating  the  most  recent  temper- 
ature values  exhibits  the  greater  stability. 

Figs.  25  and  26  show  the  convergence  behavior  for  DHT2D2  which 


employs  the  horizontal  line  formulation  and  uses  the  most  recent 
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temperatures  values.  These  curves  show  a further  improvement 
in  both  rate  of  convergence  and  stability.  Of  the  three  two-dimensional 
header  thermal  models  developed,  DHT2D2  which  uses  the  horizontal 
line  technique,  exhibits  the  best  convergence  and  stability  performance. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

A computer  model  has  been  developed  for  simulating  TSB  in 
thin  film  diodes.  The  model  performs  a one-dimensional  electrical 
and  either  quasi  or  full  two-dimensional  thermal  simulation.  The 
onset  of  the  TSB  transition  from  a low  conductance  to  a high  con- 
ductance state  is  assumed  equivalent  to  the  diode  achieving  a 
maximum  temperature  of  700  °R.  Further,  the  TSB  delay  time  is 
assumed  to  be  approximately  equal  to  the  time  required  for  the  diode 
to  reach  this  critical  temperature  in  response  to  a constant  current 
reverse  bias  overstress.  Simulations  are  defined  through  271  para- 
meters which  specify  the  diode  design,  thermal  conductivity  perturbations, 
and  simulation  control  parameters.  Execution  times  are  usually  less 
than  two  minutes  on  a Univac  1107  system  and  the  program  has  a 
graphic  output  capability  to  facilitate  data  analysis. 

The  diode  electrical  model  maintains  avalanche  breakdown  by 
varying  the  electric  field  maximum  value  to  satisfy  the  avalanche 
b~eakdown  integral.  The  junction  electric  field  profile  is  evaluated 
from  the  electric  field  maximum  value  and  the  net  space  charge 
density  in  the  respective  doping  regions  of  the  diode.  The  net 

I 1 

space  charge  density  accounts  for  space  charge  effects  in  the 
junction.  Bulk  region  electric  fields  are  evaluated  by  assuming 
that  the  currents  in  these  regions  are  primarily  drift  currents 
aiid  determining  the  electric  field  required  to  support  the  diode 
current.  Coupling  between  the  electrical  model  and  thermal  models 
is  maintained  through  the  thermal  dependence  of  the  avalanche 

i li 
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ionization  coefficients,  mobilities,  and  intrinsic  carrier  con- 
centration. 

Four  different  header  thermal  models  were  formulated  and 
characterized  with  respect  to  stability  and  convergence  behavior. 
The  full  two-dimensional  header  thermal  model  with  horizontal  line 
formulation  and  incorporation  of  the  most  recent  temperature 
values  available  demonstrated  the  best  performance.  All  of  the 
two-dimensional  models  eliminated  the  inflection  point  previously 
observed  in  TSB  delay  time  data  obtained  with  the  comprehensive 
quasi-two-dimensional  header  thermal  model.  The  quasi  and  full 
two-dimensional  thermal  models  agreed  well  for  high  currents,  but 
differ  appreciably  for  low  currents,  as  expected. 

A limited  number  of  simulation  results  were  presented.  These 
include  delay  time  as  a function  of  current  density  for  different 
header  thermal  models,  diode  temperature  boundary  conditions, 
avalanche  ionization  temperature  coefficients,  and  with  and  without 
bulk  region  electric  fields.  In  all  cases,  the  model  demonstrated 
the  expected  qualitative  results. 

As  for  recommendations,  a comparison  between  the  simulation 
results  and  experimental  data  for  SOS  diodes  is  required  for  a 
qualitative  evaluation  of  the  model.  This  comparison  would  also 
allow  further  refinement  of  several  of  the  model  parameters,  e.g., 
TSB  temperature,  avalanche  ionization  thermal  coefficient,  etc. 
There  are  also  several  types  of  simulation  which  have  not  yet  been 
performed.  These  include  various  epitaxial  configurations,  thermal 
conductivity  perturbations,  and  changes  in  the  diode  dimensions. 


Several  improvements  or  additions  to  the  diode  model  may  be 
desirable.  First,  evaluation  of  the  reverse  saturation  current 
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would  eliminate  the  necessity  for  a TSB  temperature.  Rather,  the 
diode  voltage  could  be  monitored  for  the  onset  of  the  TSB  transition. 
This  refinement  may  be  required  to  obtain  agreement  between  the 
model  and  laboratory  results.  A more  accurate  numerical  algorithm 
for  evaluating  the  avalanche  breakdown  integral  would  further  enhance 
the  convergence  of  the  diode  electrical  model. 

A primary  application  of  the  model  developed  here  is  the  pre- 
liminary investigation  of  TSB  behavior  in  diodes  for  the  purpose 
of  defining  simulations  for  more  comprehensive  diode  models  which 
require  a corresponding  increase  in  computation  time. 
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APPENDIX  A 

DEMONSTRATION  SIMULATION  OUTPUT  LISTING 
(SS-1) 


DESCRIPTION  Page 

Simulation  Parameters  A1 

Thermal  Conductivity  Perturbation  Factors  A2 

Transient  Data A4 

Diode  and  Header  Temperatures  at  0.330  usee A5 

Transient  Data A7 

Diode  and  Header  Temperatures  at  0.935  usee A8 

Transient  Data A10 

Diode  and  Header  Temperatures  at  2.25  usee All 

Transient  Data A13 

Diode  and  Header  Temperatures  at  5.55  usee A14 

Plot  of  Maximum  Temperature  Versus  Time A16 

Diode  Temperature  Profiles  A17 

Plot  of  Diode  Temperature  Profiles  A19 

Plot  of  Impurity  Profile A20 

Plot  of  Initial  and  Final  Electric  Field  Profiles  ....  A21 

Plot  of  Diode  Voltage  Transient  A22 
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APPENDIX  B 

COMPUTER  LISTINGS  FOR  DIODE  MODEL 

ALGORITHM  Page 

Subroutine  BANDA6  B1 

Subroutine  BKDEPL  B3 

Subroutine  DEPV B5 

Subroutine  DHTEMP  B6 

Subroutine  DHT2D  B8 

Subroutine  DHT2D1  B13 

Subroutine  DHT2D2  B18 

Subroutine  DOPLG  B24 

Subroutine  EFIELD  B25 

Subroutine  EMOBS  B28 

Subroutine  EPROF  B29 

Subroutine  HMOBS  B31 

Subroutine  INALFA  B32 

Subroutine  IONCOF  B34 

Subroutine  LPLOT  B35 


Program  TEMPTRI0 


B38 


SIMULATION  PARAMFTLIJS 


THERMAL  CONUUCTIVITY  PERTURBATION  FACTORS 


£ C cc ccccccccccccccccccccccccccccccccccccccc 


2 ^ 


2 frK‘iTfrK*K-|rr(*-|<-K<i*-|r^^rK-ir'rK,rK*lrprrK-K-K-|' 


- 1 !£  ^ -Sli  * i i * 55  5 * * * i!£55!5  * * * * **  5 * * «c  cc  <r  cc  <r  cc  <t  « cc  <t  a:  r*  r»  r*  r-  r*  r-  f-  r*r*  r*  r*  r*  r*>  r* i*.  r*  r-  r*>  t*  r-  ^ r<  rv 

O C **cOCCOOOCOOOOOOOOOOCOOOOCOCCOOrcOorrnoArn^n^^^i^^m***  + ♦♦*?♦♦♦ 


f ISo822eccccSeSSeccsec2SSe2SeSccSeeoeSS88Sc?SeScSSgS2SS6c2?cS2??S?  vV(\4<V<VCV<Vj<V(\ 

A?^4miiUiA^iAA4AAAAAAAiiiAiLiinininiinninninmnnniiii? 


(\  <V<X  PV.OI  <V(M  (\i  CwOJ  o,  <MCSrf(\j  CVJ  r\«C  <\  Cg(V<V  (\  CV  O.  Og 


CMr  J * »ror»  ®0’  o^<m*o  tJO'  o^ono^m^r-  <00*  o~.<\i*ojr  tONor-coo^  o^iojin^  modern*  in  c0>OHiMm  ♦ . 

KM  «.M.»<M.NIB<OIIMft<R»WR>»lrtK»K>K>->»»»J»,3a»»,r.7inSS!ES57.5!?i,CJ*“2’Z*?<5S?SPf)£; 


THERMAL  CONDUCTIVITY  PERTURBAT U»N  FACTORS  (Continued) 


c c cc cccccooccccooccococccc 

c CCCCC.CC.CCCCCCCCOCCCCCCCCC 

cccccccccccccccccccccccccc 
c c cc  ccccrcccccccccccccccccc 


w C~  - CCCCCOCCCCOOCCOOOCCOCOO 
C CCOCCCCCCOCOOOCCCCCCCCCCO 
C CCCCCCCCCCCCCCCCCCCCCCCCC 


-f  .T  -TX  U »/■ /■  T\T\T  \f  iT  IT  tT  i"  if  if  iTiT  if  Iflf  .f 

•-  c ccccccccccccccc  ccccccccc 

3 J ^3  3 ^3*  * 3333  33  33  3 cf  ^ a Cf  33 

- C S TCdCCCCCCCCCOC^CCCCCCCC 

j-  jurirjrininknuiirkDinnir.jOininLOu^ 

o'  <?€'&<?(?&  o*  !T  cr  c a a c cr^cyo'  aa  <y  cr 


r- r*.  r*- r-' h- r*  r- r- r*  r*  r~  r- 

♦ * ♦♦♦♦*■«■♦*■♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

c-  c cccoocoocoooococcoc  = oooo 

- r cccccccccccccccccccccccc 

- CCw-CwCCCCCCCSCCOCCCCCiCCC 

- - cccccccsccccccccccccoc 


'\i<\JCV.(V\j<\j(\i<\jCV  CVCHOxOaCV^CMCSi 

- =--ccccccccccecceccccc 

• t ••  t I I I I I I I I I I I I • I I I I I I I I 

c:cccccocc:ccccooooocc 
:::c3cco^cc:ccococ5c5c 
-4  ; •»  x >*4 ~ c x<\4.  c 

' t f 3SU'  c 3 Cf^cCCXO  C'C 

• - •^•O>0^OO*0  0*0K»0«0*0l0>0»0f0>0<0f0-')tf 


TRANSIENT  DATA 


s x -ir  ctr 
t :ch4c 
CT  --(TU- 


<\  (\  cv<\<\<\ 
c c ccc  c 

I I I I I I 


1 


c cr  crccccacr 
- c aaaaa 


tr  * ccc— x<\. 

S f-  <\N  K<\ 
> ^ 


X iT  \T  iT  X 
X CCCOft 
X X XXXr* 


- C C CCC 


<V  <\<?<\*iT 
< J*  X CC  C C 

•o 


1“  : roc:: 


- 

a «\x  of-*r- 

c-  rr-r-xx 


f*  r*  r**  ^ hf' 

*.  :.  ccrsc 

i i I i i i i 

►-«  r.  ccccc 
*-  •rj'xirjvr 
- / ■;  ///./ 


• i i i i i 

- - /ic  «ic 
" — C 'J>0 

r.  - —<\j  rs/o 


'Ni*.*  t-v. 


* 


'•ft'i  ,AiU'*u3  . (.ViO«tiS  .^oo*f\3  ,^.in*n''  .3u»*>403  .^ftfi*n3  .3nn*n3  ,3no«n3  .3f»n*ri^  ,^nn*n3 


, *-  r rtr  irr  rriTfririr  r r irtrir  r rtr  ir  r irfricrfrirrtrirr  irfr  irrirririr  r icrr  f-  rtr  r irr  r tr.r  r r **  r r r ir  r r rr  r *■  r r*  **- 

1 ccccccccccccccccc  c ccc  cc  c ccc  ccccccc  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

I » MMH444HH44M44444444444*44HH4444<*M44H44M444««M^M4HMHiHii 

C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCT1.ee  - 

X C CCC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ccccccc ccc 

<\  r r ir  nrjr  r r r ir  r r »*-  t*-  r*-  r*-  r t*-  r*- re  k*  k- r*- r*“  k re  r*-  r**  rere  re  r t*“  r*-  k k*  r*“  r*- 1*-  k t*- tf  r*- k*  k- re  r r*- r r-  r*-  rrrtrr.rmrrrr.rrrrrn'  * nr  r 


re  re  to  re  r*-  re  to  re  t*- 1*-  to  re  to  re  to  to  re  re  re  to  to  re » o to  re  re  »o  rote  §r.  tr.tr  re  to  to  to  re  re  to  to  »e  to  to  re  to  re  re  to  re  re  to  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  **.  r * to  re  to  t* 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccce 

444444444444*44444444444444444-r4444444444444444444444444444444444444*t4*- 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccr 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccrcc 
K K-Kf-K  rrrrr.r  n^nr  r fr  ter  e-t*-t*-  mr  r.rrrr*  r rKrrri'rr'rr-efrr  k k t*-  kk-  k Kf  r*-  r*-  kkk  k KKr*'r*KKKKt*rr*  r 


t»*  t*“»r  nr^iricmr  K*K*K*K*ir  r*-t*-|rr*‘  ir.  nr  *r  nr  tr.K.K  r nr  tr  ir.nr  nr  nr*  r.  nr*",  nr  t*-tr  nr  T tore  ft*-  *r  retOr^t**  t**  t*-tr  •*“  t^t*"  f tr  •*-**- 1*-  f 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

44444444444444444444444444444444444444444444444444444444444444444444*444  « 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccr 
t»*  re  re  re  retire  retire  »r  ieie  r.  re  re  re  K*rr*  re  re  re  re  re  re  re  terete  K re  ter  re  re  re  r.  re  re  r re  re  r^ffie  re  rie  re  re  »e  re  re  re  re  re  re  r’rrr'r*"  re  re  re  r re  Krr  • 


re  re  re  re  »e  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  »e  re  re  t e to  k 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
4444444444444444444444444444444444444444444444444444444444444444444444444 
ccccccccccccccccc  ccc ccc ccc ccccccc ccccccccccccccccc  ccccccccccccccccccccccc 
cccccccccccccccccccc ccccccc ccccccccccccccccccccccc ccccccccccccccccc cc ccc c 
to«OfOt*vototororereteiorotorOfOtoroioroiorororefOiofOrororom^rojOforoiontotoieiorotororOfOtoiororoioioterorororotoroto^iotem'Oio:eio:Oroto./- 


re  re  re  re  k-  re  re  re  re  r*-  re  re  re  re  re  nr  re  re  tr  re  re  re  re  re  re  rerere  re  re  re  re  re  re  terete  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  r*-r*-  nr.tr  f 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
44444444444444444444+444444++44^444444444444444444444444444  4444444*4**44  4 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c c c ccccccc ccccccc ccc ccc ccc ccc ccc c ccc cc ccc cc ccc c c cc cc cc cc ccccc cccc cc  cc ccc  c 

re  10  re  to  10 1*- re  re  re  re  re  re  re  re  k*  re  to  re  .ore  to  re  re  re  re  to  ro  tore  tore  to  re  to  to  re  re  to  to  re  to  to  to  re  to  ret*- to  re  re  ret*- to  re  re  re  re  to  to  re  K-r*- ter*- re  re  re  tore  ret*- t*r*- 


k- tote  tore  re  tete  tore  re  re  re  re  re  re  re  re  tote  re  re  re  re  re  re  re  tore  re  re  re  re  re  to  re  re  re  rete  rote  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  re  terete  re  re  re  K-r^ 
cccccccccc ccccccccccccccccc ccccccccccccccccccccccccccccc ccccccccccccccccc 
4444444444444444444444444444444444444444444444444444444444444444444444444 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccc ccc ccccccccccccccccc ccccccccccccccccccccccccccccc cccccccccc ccc cccc  cccc 
mere  ie»et»*  ter  terete  r remoter  remote  tore  terete  r rero  tor.  re  re  r- re  retire  reterere  re  re  tore  retire  re  re  re  re  rerere  re  re  re  terete  tore  rerere  r'-K-met*'  r 


re  to  re  to  re  re  to  to  to  re  re  re  to  to  ro»e  to  to  to  re  rote  to  to  to  to  to  tore  to  re  to  to  re  to  to  re  »e  re  rote  to  re  to  to  to  to  re  to  re  re  re  re  tore  re  tore  re  re  re  re  re  to  re  re  re  re  re  rorr^ere 
ccc ccccccc cccccccccc cc ccccc cccccccccccccccccccccccccccccccccccccccccccccc 

44444444444444444444  444444444444444444444444444444444444444  44444444444*  4 4 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
ccc cccccccccccccccccccc ccc cccccccccc ccccccc ccccccc ccccccccccccccccccccccc 
re  re  to  re  re  re  re  re  rr- re  re  re  re  re  re  rr  re»ete  terete  K-  r**  r»-  rf  re  re  re  re  re  re  re  re  re  re  re  re  re  re  rr- re  re  re  re  tore  tee*  re  re  re  re  re  re  re  re  re  re  re  re  re  tere»e  ter.  re  re  re  e-  **- 


to  re  re  »o to  rote  k>#o  re  re  ro  to  to  ro  re  re  re  to  ro  ro  10  to  to  to  to  tore  to  ro  ro  ro  ro  ro  e re  ro  rororO  ro  ro  ro  to  ro  ro  to  ro  ro  ro  ro  ro  re  to  to  ro  ro  ro  re  ro  ro  10  to  ro  to  to  ro  k r'rrr, 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccec 

444*444444*4444444444444444444444444444444444444444444444  44  4*4*44  **4*444  4 

c oc ccccc CCCC CCCCC  CCCCCCC CCCCCCCCCC ccccccccocccccc c ccc ccc cccccc  ccccc ccc ccc 
cccccc:cccc:ccc:ccc:ccccccccccccccccc:ccc:cccccc:cccccccccc:c:cc::c cc  c :: c 
re  tore  torero  teK’r.ro  re  torero  ro»e  re  re  roro  re  re  re  re  re  re  rerere  rote  re  re  re  re  rcreret*"  tore  re  re  rerere  k- re  re  re  re  re  re  re  rerere  re  re  re  re  re  rr.K-.rerereK  rrer 


to  to  rO  to  rO  ro  to  te  te  to  to  to  rO  to  re  re  to  re  re  to  ro  to  to  to  to  re  te  to  to  to  to  m to  m m to  re  to  re  to  rO  re  to  rO  to  re  re  to  to  re  e*  re  to  re  to  re  re  re  to  ro  to  re  to  rO  re  re  to  re  ro  m te  re  to 
cccccccccccc  ccccc ccccc ccc cccccccccc ccc ccc ccc cccccc  ccccccccccccccccccccccc 

44444444444444444444444444444444444444444444444444444444444444444*4444444 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccceccccc 
ccccccccccccccccc  ccccccccc  cccccc  ccc  cccccccccocccccccccccccccccccccccccccc 
IT  * IT  tr  If  rr  tr  r rtr  r rr  nr  r r tr  rr  IT  r r r K >r  r r nr  tr.ir  ft-  r ir  tr  nr  tr  K ir  IT  r IT  r nr  rr  r r tr  r *r  tr  tr  K nr  IT  tr  r r ft.  r ft  IT  If  r r IT  r tr  * tr  r 


e.rOtOtOrO*OtOrOtOrerO  to  rorororerororotoroioroiotoiorotororororetororororororororoe'rororore  terete  re  rotoretotorororetotororo  rote  to  tO'OrerorO^r'  — 
ccc  ccccccccc  cccc cccc ccc  ccc  cccccc  ccc cccccc cc ccc  cc  cc ccccccccccccccccccccccc 

*44444444  *444444444444*444444444444444*44444444444444444*  44  44444444*4*4*4 

c cc  ccc cccc ccccocc ccc ccc ccc ccccccc ccccccccccccccccc  cccccccccccc  ccccc  ccccc  c 
ccc  ccccccc cccccccccccccccccccc cccccccccccccccccccc cccccc ccc : cccc ccccc  :cc  : 

e*  re  r»- ^re  re  tore  tee' re  re  terete  rote  e*  re  e' re  re  re  tore  re  rerererorere  ter- rer*  tore  r- retire  re  re  re  re  re  re  re  re  re  re  re  re  re  e*rere  re  re  re  rre*  tore  re  re  e*  re  *^*00*  *• 

r tore  to  tore  reroute  re  re,  tore  re  re  terete  re  tore  re  re  tore  re  re  re  to  tore  re  re  re  re  re.ro  terete  rerere  tore  torero  *e«e*e«»‘rOierO  tore  »ero  tore  «o»o*ero*e  re  «o»eie»o^ 
cccccccccccc  ccccc cccccc  ccccccccc  cc ccccccccc  ccccc  c ccccccccccccccccccccccc c 

44444444444*44444444  4444444444444444444444444*44444444444  44  444444444*444* 

CCC  CCCCCCCCC CCCC G CCC COCO GO ccc cccc ccccccccccccccccc ccc ccc cccc cccc ccc ccccc c 
ccc  ccccccccc  ccc ccccc ccccc  c cccccccccccccccccccc  cc  cc cccccccccccecTccrcr  t :c c 

tr  tr^K-K-^^tr-rtr  K*  I*- rr*K  tTK  tr- re  t»- tr-t*"  e- re  t^KK- ter*- e- re  K-r  re  #TK  re  rote  Ktr- re  err  trrrtr-r  r»“  re  ret*- rr-rrtr  ter  K- ret*- trtr  re  tt-K-K- re  re  re- rt- re  4- t>- 


rO  *0  *0  *00  *OTO  *0  «0  rOrO  rO  rO  tO  rO  tO  *0  to  *0  rO  rO  fO  *0  tO  *OrO  tO  fOrO  *0  rO  *0  rO  rO  rO  to  rOrO  rO  tO*">  tO  t0r0r0r0-0^r0r0r0r0t0*0r0'0r0r0  tO  0r0*0t0r0r")-0*e 

cere ccc cccc ccccccccc cccc ccc ccccccccc cccccc ccccc ccc  cccccccccccccccccccccccc 

► 4 4 ♦ 4 4 4 4 4 4 4 4 4 4 ♦ 4 4 4 ♦ 4 ♦ 4 


T CC  CCCCCCCCCCCCCCCCC  CCCCCCCCCCCCCC:CCCCCCCCCCCC-"""H**-HH— -CCCCCCCCC  c . c 

-cwcccccrccc  cc c cc  cc c ccc  c c c ccc ccc c c c c cc ccccc cc c c c c c ccccocc::ccccccc:cc::: x 

ee,i*‘ror*-«*-^-rOt*-re-->'OtOe:rOrereto**'rererOrerOr*-,trit*'tOrerOtOrOtOt*-,r^rO'Or*-toi*'rOr*'te'OtOe'K'tOrere^reteret*'te*OreK'rer**rOK'K'*^r*-retee  r *.m  .t* 


e*iooe*vieoeeeioeoo>0'oee'no*vneoe>otO'oeeto  0>0t0t0,0re'0*0t0*0r0*0*0'0r0'0'0t0'0f0o,0,0or0t0>0t0«0t0’0t0f0'0*0i0«0!0'0’0 

cccc:ccc:ccccccccccc:ccccccccccc:cccccccccccccccccccccccccc::cccc:cic:cct 
4 44444  444  44444*444444444444444444444444444444444444444444  44  444444444*444  * 
• - * c— —————— <v*o»e  tea  3 c— re* ccc—  me  crier*  — * — r^recec^ir  —a  ir  w <vr»r.  excite©*  evo-x  re  j i\ 

-cccccccccccccccccccc  cc  c ccccc  — ————— <\<v<v<\.ro tea  air  irxr*r*ccc-  cccc  c<?  CO  crxr-r-xx  xirv  if  if>  a a .?  3 

•-toe  retOt^  r*-e'r^*e*ire  tore  rotor*- rorrro  torero  rorororoterororororerororetotororotorerororeiororo  tore*  ^ a a a torero  re  re  re  rerOt*-  toe-  r*~  r*-fomet*-  * 


^«<tOtere»e*ot*'te^5teio*e*ot*-t*‘re'ete'0'e'Ore  '>iOj<vor'N<f>r‘M<vorfvr.\jorrs(Oi.M<MOt<V'VO<  vjrs,  w\.v.\ 

cccccc  cccccccccccc  ccccccccccccccccc  cc  - cccccccccccc  ccccccccccccccxcccc  ccc c 

I I I I I I I I I I I I I I I I • I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I • 

ccc  cc  c cccc  cc  ©cccc  :c:ccccc»i<vfl  ©a*o„c  c»tn  rj«(vx  ea<r«\x,  ca  x rvx  ca  x <\x  rac  *\x  rax  <\ 

oeorxea  r fvc ojr. o*x< oj cc <-4-0  c.-rxorx 00c— — ojcsjcNjtetear a c 000— <\i<\.<veieaaair.if  oxx. r-r  ccoctcc* 
xx  — — orcjoiroe'a a a txO'O1  — — — — — ———————— ———  — a, 


:-fMe 3X  or- j c — <voa».o  oo- wOJ'C—  rvrt»^a*n  oe- ta c—  orte a ooo> c — <voa«r xr- lO'C- csjto 40  jmC'c  — <v<o 4 
— -t'j*vj<M >j*j*o»oro>o-oo‘Oo*ooaaaa  aaa  t aax»  r\s)r>s\x\S)S)  ox  oxx  j 0 ox 


O OK  K OO  OO  K.OKK  00  0.0  00  OOOOOOKOO 

c cccccccccccccccccccccccccc 

4 4-*4*4444444444444444444444 

C CCCCCCCCCCCCCCCCCCCCCCCCCC 
c cccccccccccccccccccccccccc 
K*  OOOOOOOOOOOO.  OOO  OO  KKlT  OO  O OO  O 


o,  ooo.oooooo 


o o oo  oo  o oo  o 
: ccccccc  c c 


O Or-K-jnrK-K-OK'OOOOOO-'OOOOOOOOOO 

c cccccccccccccccccccccccccc 

4 ♦«4«-*^-»4444444444444444444 

c cccccccccccccccccccccccccc 
c cccccccccccccccccccccccccc 
K’  ifrirn<-nrK,ririrrrrrrrir oo ooo oo o 


o o.o  oooooooooooooooooooooooo 
c cc  cccccccccccccccccccccccc 

4 *4  44444  + 444444444444444444 

c cccccccccccccccccccccccccc 
c cccccccccccccccccccccccccc 
O OOIOOO.OOOOOOO.OOOOOOOOOOOOVOO 


o o.  o o o o o o o o o o o o o o o o o o o o o o o o o 
c cccccccccccccccccccccccccc 
♦ *44*44444*4444444444444<4<44 
c cccccccccccccccccccccccccc 
c ccc  c ccccccccccc ccccccc cc cc 
0 00000000000000000000 000000 


10  o o ir  o o o o it  o *0  o o 10  o ooo  o o »r  00  o tom  o 
c cccccccccccccccccccccccccc 

4 *+  444444444444444444444444 

c cccccccccccccccccccccccccc 
o ccc c ccccccc ccc ccccccccoc cc 
<«rnrrm»-irr»rirrmr  00  »*-  00  rr  o.  00  o 


o o. o o o o.  o>  o o o o o o o o o 00  o »oo  00  o 00  o 
ccc  cccccccccccccccccccccccc 

4 ««t444444444  44444444444444 

C CC  CCCCCCCCCCCCCCCCCCCCCCCC 
c CCCCCCCCCCCCCCCCCCCCCCCCCC 

r o.  ir  o.  rr  r.nnrirrr  o.o.  o oo  oo  o oo  o 


o o.o  o*o  k". Kitf.K) o,o; o, o, r> n o, o>o. o, rr, o>»o o< rf, o. 
c cccccccccccccccccccccccccc 

♦♦♦♦  + 44  + 4444  + + 4444444*44+44 

C ccccccccccccccc cccoc coco  c c 
ccc  CCCCCCCCCCCCCCCCCCCCCCCC 

o o.o o o. o o o. o, rr, rrnrrrrr irr oo 


O O-frOO.OOOOOOOOOOO.fOOOJOOOOOOOO 
c CCCCCCCCCCCCCCCCCCCCCCCCCC 
♦ ♦♦  + 4 + + +444444  + 44+4444+4+4  + 
c c ; C CCCCCCCCCCCCCCCCCCCCC  cc 
c cccccccccccccccccccccccccc 
"■  O O O OO  O O KKK  r/(TK  O O O KK'KKK  ^ ITK  O 


o*.  0,00  ooooo  0000000*00000  oooo»no 
c ccc  ccccccccccccccccc  cccccc 


c cccccccccccccccccccccccccc 
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U.  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

- — - - — --jjcca:a;«tccaa:4i(\<\r^r'rs-r-r' 
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* — -4~-  — -4—  — — ~ 

- c CCCCCCCCCCCCCCCCCCCCCCC 
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- - — — — <\<\  <v<v^fTFnj  a a a ax*  siT-r 
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X C CCCCCCCCCCCCCCCCCCCCCCC 
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k »r  k“  k"  k"  --'K'.  fc»c  fOK-.w.  fT  tr  K-.  Kiioir  k*.k:  ir»r»r  rr  k- 
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— c CCCCCCCCCCCCCCCCCCCCCCC 
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♦♦♦♦♦♦♦♦444444444444444444444444444444444444444444444444444444444444444«4-« 

c --cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

If  Cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

<\  K"  K K"  K K"  K"  K‘t^K‘K"r'"K'“K‘K  K'K'K'.  K'l^K'l^'K'K'K'K'K'K-  V r K K K-  lrlflrlf  fTK  K’K'K'KK’  KK'K’KKK 
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c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

^ cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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££ccc cccc cccc c cccc ccccccccccccccccccc ccc cccccccccccccccccccccc cccc ccc cccc 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

♦♦♦44444444^4444444444-4444^4444444444444444444444444444A4444444444A4444-*4-* 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
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^r'  tr4fT.tr. K!  »T  K*.K*,KiK-.K'.K1K;K:ir)Kl»r»0»rirjfr  ^*^^»r»OK;»C»OK;K)ir.iO»r}»0«OK'.K*.tOK:KHOK'.K>rK>CKJ  KMTKJK}K*,»0KM'',*',>»0K>K>.K‘.K-.w\»0»C 
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O ■0^)O‘0‘0K>O*0>0«0-0«0<0K)'0<)K)K)K)K1K)’0K)'n>0'0<0«0K>*0K>«0i0f0Oa0r0K)K>«0-0'0K>K)K)<0K)K)K>f0O->0'0K>^)  fO0*0*0KV0  K) iO K) O-n K) rO K> 

i*  , - --cccc:ccccccccccccccccccccccccccccccccccccccccccccccccccccc:cccccccccc:c 
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c cccccccccccccccccccccccccc 
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c cccccccccccccccccccccccccc 
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«r  cu"  r^.o-  rrcc,f\j^>£r*-x>ia  ircr  o-  c 

<-■  s :*.*<r<\<v— c c<r  j:  f^tr^rrtvco-ac^^rr  — c 
kT  ^ ^ a 3 3 jjtjj  rrnrrr  rr 


(\  m ‘wOeVAr  .\vMf\rsrrkj(v(\fvrg(\ 

- - c ^cccccccccccccccccccscc c 

• ••♦•♦•  I I II  II  I I I I I I I I I I I I I I 

-»C'\iXc.»«<V£c»cc<\.XC4r<r'vxc*x<vt.  c 
0 w • c~-n<v\t<\rrrir.r;*:*jn»r,XX>or*-f*'flO«<£0  0 o 

<v  ^ rr  rr  -r  rT  rr  rr  rT  rO  rr  IT  rr  tr  fTyr  rT  rT  fT  rO  rT  rT  rT  rr  trrr  ^ 


^ 3r^ -oc  c— or* xa  c — 

OJ3XJOO.O  0-0  V)  M'0'3>0'0‘3'3’3'0|CO 


C C cccccc----— --ccccccccccc—  — — — -cccccccc 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

♦ ♦ ♦ ♦ ♦ + + 

^ Kc«c^^ffCeCM^'K■<\c^^c^(\^«^c^'£ccal^r^or^aaoc^ 
Uh^4(r(N^^CrCCCCCC«h^iri^<\(\-- CCCCCff«M'3lfif(\ 

c acoaecra^^— — -cff{racraoaoec-^-'--cTaoao  cao 


i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i • i i i i i i i 1 1 i i i i i i i i i i i i i T 

t (cumccictffacraaatyaffcracccc'off  ---^^-K-KrrKK-rvinrir^vC^ttffctvafv 


W c c ccccccccccc cccccc cccccc ccc ccc cccc ccc cccccc ccc ccc ccc ccc ccc 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I • I I I I I I I I I I I I I I I I 

a cocaaaaaaaao'aao'aaacaa accraacaaco'acoaaoaaacaoaaaaa-c caaaaaao-aa 


•o  »0.^y5.^lOrn^»0irr0^f0.,0K)f0^!O.,Of0t0JOiO^)tOr0rOi0f0K)K)f0»0»0»0i040K)i^i0r0K)K)»0!C.'Of0.',0r0^l0r0K)K;i0.'0^.,Of0rC 
cccrcccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
♦ ♦ ♦ ♦ ♦ ♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦■*•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  ♦♦♦♦♦♦♦ 
— ------- xcctrccc  affaaaacrffa  craaa-c<\.<\c\»rK-»>-r»rK-  rrrirc  — »■  ---xxx  xxitk-.gccc  -£x  v£k- 

^ :t.» ** ;*:» :» ^^■*xxxxxxxxxxxxxx»cx'c'£xx£xx«£xxr^r''r^cc0  <\.j 

<\  <\ c\c\ f\c\(\<\c\f\py,c\c\(\ fs(c\(\f\c\c\<\(\(\(\<\<\f\rs.f\c\(<\f\(\c\ c\<\ (\f\<\(\(\(\(\<\c\c\(\c\c\c\ (\(\<\r\(\.<\<\<\»Ofr 


<\  (NJ<Vf\.C\  (\'»OKVO»r'fO  K'fOK'  KVOKIrOK'fOfrPOmtr  tOfOKIrrK'fTfr  K>OfOK^K'fOfO»OK'P^K^fOK1K'fnfOfOK'K1K'  I0»0»0»0*ri0 

c ccccccocccccc coecc cccoccoccococccccccccccccccocccec cccccccc 

♦ ♦ ♦ ♦ ♦ 44  4*t>44  44^««^4«444**«4444444444444444«*4f«4*444  ♦♦♦♦♦♦♦ 

r-  r*  r-r-  r^r-r^r-  c\K-KK-»^i»'K-K»*'K-KKK»rKKK-Kr^ccccccccccccccc<cta<rcc<i.'£v£«£N£^3j»rK-cjc 
x x xxx xxx c cccccccccccccccccc  — ------------------ -f\rv<\c\c\ k.k  3 3ir«ic<\ 

j <r  c aaaaa  cr----- 


c cccccccc-cccccccccccceccccccccccccccccccccc  cccccccccccccccc 
♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  ♦♦♦♦♦♦♦ 
x i^ecc^(\K'^u-v£(r  C“c—*o;»xxf*‘i''«rCT  cc**<\<\»o^x  rr^cce  c— 3-j~r^cr a c— cvK-^^r^cc  ccv^xcr 
c c cr^f^f^r^r^r^r^f^f^cccrxccxcc'C‘XxxcTTcra?aa'a-cya'aJ-  :ccccc  :ccc c— — — — — — — — c\ <\ <\ <\ <\ -* ?r 

f\  C\  C\r\  C\.C\C\,C\<\f\<\.(\C\.C\  <\(\(\<\0»C\(\C\C\r\(\f\(\C\(\C\C\fV(\f\^r-K'»r.ir.K1K)K-.K:K-  ^r-  K'K-  K'  r-K-,K' K'  f^K-K'r^  f^ 


— — — -*cx.xx.  xx  x.xxxxxxxxxxx<\<rcaa-(MTcr(7<Ta,<ya'a40»OfO»OK'40»o(Ecccccccc  <\<\.xxc^cac 
A.c«itc.c«occcoococcocaaa<Taocyo<ycccccccccccccc<\<\c\<\<\c\<\KKKK-rctrtf.xxc<Ti/r‘ 

"r  if'r«K,r<ri<,K'KK,lrK‘rK,i<,|f«<'rK,KK,34  JJa3J4  13iJS433333a4  1133  4 333344/(f 


K-<^K-rr<^rK-<Kr  *r»r »r  k-.k- K-xinr ^^K*»r^x,po*o»OK}i^»o*oKiiOFOfrrOfOK',K'K-ioK'.»of*->i040K>K-K-.»OK'.*oio»o»o 

c cccccccc cccc  cc  oc.cc  c coccooccccccoccccccccccc  ccccccc  cccc  occc 


— c ccccccc— coc cccccc c ccccccc oococc ccccccc— cc cccc—  ■ 
ccrccccccccc: cc  cc  ccc  cccc  — ccc  ccccccc  cc  ccccccccccc 
♦ ♦♦♦♦♦♦♦♦♦♦♦♦■♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
jr:^-4  — ajc  cr-cf^cr  r.rxo  3 3 30  o a xa  x xa  * <\  rv  — .a  a k-  x x x a 

t c»o*»nx^ir.  cr^x<\<cc®~  ^ *—<r  x rr -c  r-'X  — c<\x  -e-  rocr-  3 <\  cxr  ct'-a:  — x.  occ* 
— <0  c or^r'-xx®  -•3'<o cr-- r->c  xxxx :t  ♦^^»o»oxa®^xxx)ir>«r>^»^-t'r'«^oNsxxx>— • 


♦ ♦44 
X 3 X - 

x<vr*  c 
-^rgrvjx 


r*»  r>*  r*r*r*fs-f*r^rs-r*r*r»  r,‘r^r^r>-r*‘r*r“r^f^r**r^f^r»r‘r'r'r-r-r>f^r^r*‘r^r,*r‘^r^r*r»r'r«‘r‘r^r^r'r‘r‘r>h-r*r‘r*r‘r,'r‘r* 
: ::c cccccc  c ccc  cccc  c cccccccc ccc  cccc  ccc ccc cccc  cccocroccccccccc 

' 1 1 1 1 1 1 l 1 1 1 1 1 1 1 1 l 1 1 ••  1 l 1 l 1 l 1 1 1 1 1 1 l 1 1 1 1 1 1 1 l 1 1 1 1 1 1 1 1 1 1 l 1 1 1 1 1 1 1 

. :uccctcccccccccc:cccccccccc:cc:ccccccccccccccccc:c::ccccccc 

X X XX  XXXXXXXXXX)  XX»X  XX  XXX  XXXXX  XXX  X XXX  XlXXXXXXXXX  XXXXXXXXXX  XX XX  XX 
* X XX.  xx.xx.x.xxxx.x  XX  X.  XX  X X X X.  X X X X X XXX  X XX.  XX.  X.  X X.  X XXX  X X X.  X X X X X X.  X X X X XX  .-•/. 


X X X X X X XX  X X X X XX  XX  X XX  X X X X XX  X X X X.  X X X XX  X X X X X X X X.X X X X X X X X X X X X X X X X X X 
-ccc  cccc cccccc ccc ccc CO cccc ccc ccccccc cccc ccc cccccccc ccc cccccc 

I I I I I I ••  I 1 1 II  I I I I I I I I I I I I I I I I I I I I I I I I I I I I I • I I I I I I I I I I • I I I I I I 

- X C VA-tOSJj/CX-MVSrrO'  .»  C X — 0 AT-  KVJC  X.  O®  f\l®  *OX  3 CX-®  AC'-  *0»©  3&  XC  ® AlCO  400-  ^ OX. 

*■.  : j sr.  z ; c S'J'C  c—  — <vo*o  * 3 xx®  or-r»  c J'J'c  c — — \<\»o  * trx.x  cr-r**  oocc  — — \ AnO*o*xx 
■ g c\  aimajMaj  \j<\i:\4<VNAm  •o*nm»o»r)*n*o»o,o«0'n*o»n«n->»n,o»  33  44.434443444  4444/1^^  XX  xx>xx  XX 


eg  40  *x®f»ac/0-  o-*a*o«x  ®r~<00'O«cgr)»M>®r'>«;CT'  o«f\i»n^x>or^  eocr  o^rgm» x of- ocxj>  aj*o  »x  cr^®cr  0-1 
j * j -t  j 33J  -'xxx.  x.x.  x.x.xx  r-r^h-r*  c«  c®a.  x,  a.*cx»®oo  ao'C'C-xo  C'O  00 


DIODF  AMD  HEADER  TEMPERATURES  AT  5.55  *SEC 


— X X c 


ir  K-ri'-K-K-rK-irK-^rK-Kr.r^^^^irirrK-KirK-^^K-irKirrri'-irtrK-irr^rc-mrK-irK-rKK-rK  rrKri«-KK-rKKi<-irK,irK-KK‘rK- 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
♦ 4 44  4H4M44  4 44<4*444HH44*4M4HH4  4 44H4<H44M4  4444444M4  444  4 44  4 4 «*4^i 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
k K-.tnr  K-K-K^r-^m^nrK'K-K'irirK  rrrr  rrr  * tr.r rrr  trr.r  rrrrrrrrp-rrKiK-rrrKrri<-rri<'irK-KrKK-Krifrrrr^ 

r.  r.  r:  tr  ic  tr  in  tr  ir  in  tcr.tririr  ini':  r:icic,io  in  ritrtnrnic  ir  trtntrtnr.r.tr  »c  »r  m it  »o  k:  k*.  k;  k:  ir  k;  »r  k:  it,  r.  r.  ir  k;  m »r  k.  inline  ir  k,  »r  k-.  it  »r  r.  k:  k- 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
4 «4^4^44  4^44  44^  44  4 4 4 4 44  4 4 4 4 4 4 4 4 4 4 4 4 4 44  44  4 44  44  4 44  4 4 4 44  4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 

c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc  cccc  ccccccccc 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
KKirK'rrrrrK'r^rrfirrKiri'ri'irrrKrrr'.irrrrKrrrrKrrrKrrKirirrrKrfrrrKKrrrrrKrrrrrrrrrK 

»r  k*K‘K‘K*jt  trir  tr-trtnrtr^^trtrinr  ir^  tr^trir.ir  tr  tf  irtrjr.ir  ir.ir  tr  it tr.tr  ^trtrtri^.li'.ir.  ir^-^r  tr_trjr.tr r iC r r W ir.  r r *r  mr  r 
cceccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
4444444444444444444444444444444444444444444444444444444444444444444444^444 
c c cc ccc ccccc ccc cc cccc c ccccccccccccccc ccccccccc ccccccccc cccc cc cccc ccccccccc 
c cccc cccc ccc ccccc cccc cccc ccc ccc cccc ccccc ccc ccccccccc ccc  ccccccccccccccccccc 
k:  icirirKir  innmnrinriririririririririririririririmiririririrmr  r.iriciririri^irirrririririr  r.icK-ir  rrrrKrrrrKrmrnrrK’r 

»r  »r  »r»r»c»oio»rir  kj ir ir  ir  ir ir »r  k.  ir  ir  f. k-.  k; m ir  k:  k-.  ir  k:  ir  r. k.  k. it  »r  k: k*. »r  iot  k:  ir  k:  ir  ir  r.  ir  k. irrr.i'.mr #r  k. k.  r.  k.  k. ir  r. ir  k.  k*.  r.  r.  r.  r.  r.r.r  r k 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccc ccccccccccccccccccc 

4 4444444444444444444444444444444444444444444444444444444444444444444444*4* 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccc ccccccccccccccccccc 
rO  .•OiOK5i^i^.'riOrO»miOr^iOrOK)fOn»OrO.'OK5»OfO»0»0»OlO»niOi<‘.»OfnrOlO»OiOiriOm<V»OnKVO»OK)K)iO»Oir)lOrOK)*0»0»OK)rOK>rr  KJlOfOfOlT 

|T  ir  |T  ir  fTir  |T  ir^fTiTfr  |T  IT  itk:  ir  ir  |T.  ir  (T.irir;  ir  IT  IT  IT  IT  IT  iTfCVT  IT  IT  IT  r.  r.  fT  K^.lTfr.irfr.  ir  IT  fT  IT  friTfr^ir^  IT  «r  IT  IT  IT  IT  r IT  r.  IT  IT  IT  IT  r IT  rr«rr- 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
4444444444*4444444444444444444444444444444444444444444444444444444444*4444 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccc ccccccccccccccccccc 
c ccc ccc cc  ccc ccccc cccc cccc c c ccccccccc cccc ccccccccccccccc ccccccccccccccccccc 

*r  »riririririr»o»ni*,‘»o»nio»otrK;io»o»oio»rioioirKMo»oio»riricio»rio»r»ofOio»r»r»rK'.»o»OK;»rir»oir»n»CK:»rio»r»o»rtCK;irir  lonoiorinoiriooirfrir 

c c ccccc ccc ccccccccccccccc ccccccccc ccc ccc ccccc ccccc ccccc  ccccc c cccc cccc ccccc 

44444444444444444444444444444444444444444444444444444444444444444444444444 

cccccccccccccccccccccccccecccccccccccccccccccccccccccccccccccccccccccccccc 
c ccc ccc ccccc ccccc cccc cococ c cccc ccoc cccc cccc cccc cccc ccoc coccccccccccccccccc 
*■  tnir.tr  nr.  rinr  (ri'-iriririri'-r-iriOKir  K'K!K*r.»rfr  r r r ir.r.r.r  rirnrir  irirr  r.inr  mnr.irir  »r  *■*•»*■  r^r  rrrrrK-  rrr  r k k r ^ r k 

n po»r»o*cirio»rtr  1010^5  f^miOfntOfr»OfO»o»0(ntr»o»o»ofOfOto»o»oio»o»n#o»o»oto»c»o»o»r»oio»o»OfOiotoK'.io»oio*rK)fOK)io»oir  mioioicioioioicroKvioirir 
c c ccccc  cccccccccc ccc ccccc  ccc ccccc  cccc  cc  cccccccccccc cccc ccccccccccccccccccc 
44444444444444444444444444444444444444444444444444444444444444444444444444 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
r.  0-tnnr.ir fnnrinriririrfrfrirfrfririrfriniririrfririrfriririr *:irir»r  irirr-.ir#riririCK-.»ri'-.ir.iriririr*rir ir»rimir»r inr »rr  rirK-  r o-rrK-^ 

r,  r>irM>iw>r.r  iritr  t^friir.r^irr, *:*,*,*, ,r.r,r,r,ir.r, trir  r.r.ir  r.riirjr.r.r.r.r.ir. 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
44444444444444444444444444444444444444444444444444444444444444444444444444 
c ccc ccc  ccc ccccccccc coooccccoccc cococccocccccc ccocccecoc  c ccccc cc  ccc ccc  ccccc 
cccc  ccc cccccccccc ccccccccccccccc ccccccccccccccccccc cccccccccccc  ccccccccccc 

k-  lOincir^iririr  K*iniri»*<»rK*.iriOK>»»'iirfO»ririr»Oin»r  irioirminmfrmpr^fTioro^.^fririr  mmriririr  »nr»nr  n'.nrr- 

*0  ki  »o  K^»TfnKiKi^^KiK:fOK:f^KiK:KifnrOKifO^»nKiK)K'rOfrro»nfnfn  ft  rn'OfOfOfn  »orn  »oK'.fo»OfOfO»OrofOrr»OK>n  mrofOtntnmfrKirr  ro»OK',»ofom»nK',»n»r 
c c ccccc  cccccccccc  ccc ccccc cccccccccc  cccc cccc ccccccccc ccc  ccccccccccccccccccc 
44444444444444444444444444444444444444444444444444444444444444444444444444 
c c ccccc  cc  cccc cccc  cccc cccc cccc cccc cccc cccccccccccccccccc ccccccccc cccccccccc 
c ccc ccc cccccccccc  ccc ccccc ccccccccc ccc cccccccccccccccccc  ccccccccccccccccccc 
K K KK  K.K  * iTKKKlTKKK  KKK  IT  KK.K  »T  KKKK.KK  KK.  K K K K KK  K.KK  K KK  K KKlTK  K KK.KK  KK  K.K  KK  ITK  K Kr  K K r K k K K K 

»0  *0K)K)»0»0»0»0^t0*0KllC»0K^K>K1»0*0l0K)K)K)K1»0l0»0»n»0f0K)KV0KlK1t0*0f0O»0O»0i0»0*0f0K>f0K:K)K1i0^t0»0l0»r*nK,)»n*n»0ir*C»0»0i0»0l0*0»r*n»nr 

c c cc ccc  cc  cccc cccc  cccccccccc  cccc cc  cccc  cccccccccc ccc ccccc cccc cccc ccc ccc ccccc 
4 *444444*4444444444444444444444444444444444444444444444444444  44444444444** 
c c cc occ c c cc c ccccc cccc cccc cc ccccoccocc CO cccccococcc ccccc ccccccccccccccccccc 
cc ccccc cccccccccc cccccccccc cccc ccc ccc  cccccccccccccccccc  ccccccccccccccccccc 
rrrrK'K'irr  K!K'ir*^»nricir»r»r»riririo»r»OfOK:irirK;»niriOK*K)ir»rrrto»o»r»rK;irirfrir»r»r.K*.K;K-.»rK;»r  rrirfCfrirfr »r*r »nr*r»r»rir»rtr»,‘.ir **“ 

•r  *r*rir K)*r»c»o»r*o**“,»c  *o»o»r‘.io*oio»o*oic»o*r)*o*o*o»o»o»o»Cfo»Ofo»o*o»rK',io»o»OK)»OK;K)»o»o»r.io*o*o»o*o»o*r *o»o»oio*;*o»c»o»oir »r»o»opo*r 
c ccc ccc  cccc cccccccccc cccc ccc ccc ccc ccccc  ccccccccccc ccccc  cccccccccc  ccccccccc 
44444444444444444444444444444444444444444444444444444444444444444444444444 
C CCC  CCC  CCCCCLC— H— C\.<\  <\<\C\«<\.C\C\.(\  C\ 

c c CO CCC CCCCCCCCCC  cc c ccccc  cccccccccccc cc  cccc cccc ccc ccccc  ccccccccccccccccccc 

rtrrrrK'rr  k-kik-k-  rrrr  rrrirr.ir  ^rrm^irr,  inr  r »n**ic  nr  k nrrr.ir  K-rrmr^rK'/  K,rrrKrmrrr(»,rKtrK-(fi<,ir 

•o  KlfOfO^OKJO’OOfOfOOfO^^OOOrOKiOiOOOK^OKIKIOlOrOlOOfOfOK^rOKVOOOntOrOOKIfOlOrOiOrOOlOrOKJK^OfOOtOOfOOnfOPOOfOfOfOrOinf'' 

c ccc  ccc  cccccccccc  cccccccccccccccccccccc  ccc ccc ccccc ccccc  ccccccccccccccccccc 
44444444444444444444444444444444444444444444444444444444444444444444444*44 
C C J cJkTlT  XCO  C C — — CVifOK'  ^VDinX  t'-t"  cc  0 r»0C«O  C C— — <\.<\  CVfT'OfOfOK)  K1fO(VC\C\— ■ — — -CCC 

c c cc  ccc  cccccccccc  ccc—  •- 4. 4. 4.  4. — — 4.,—  —.  — 4.tv<\i(W<\iCv<\jo.fVi<vfN.rv;c\((Mir»o»rirfCK;K)»r^.»r»rio*^ir  or  k.»*\  K'mricr.rA 

*r  *^10*0*00 K".K‘,K'K'r  *r*n»rtr*o»^»n»n*r»OfO»o»r:»r»r»o*oiOK'»ofO»r  »OK'4i»oro»r>K'.*rioK:tnir  »mnio*r*r  K)»r  'ok'^k-  k. 

•O  •0-0»0»0*0 »0*0<>*0»'VO'0 *0«0«0»0  -0*0  fO  0*0  <0 *0*0 010 *0*0 *0*0 *00  OTOO  *0*0  *0*0 OO 0»0  *V0*0  K)*0»0'0»0»0»0  OiOOKHOO  oo  o»o <0*0*0 «o«o»o  OiO  -o 

c c cc  ccc  ccccccccccc«cccccccccccccccccccccccccccccccccccccccctccc:ccc:c:cc:: 

4 4444444444444444444444444444444444444444444444444444444444444444444444**4 

c <Mrr*-  c<\«nr*  t(Wcc-3N  3r.>£C4N-»t  ok-.  o»c  — r-.cr^a  — oecr^r^  rxo  c<w.o-^  cr^ir  ^j-r^  ccr-c «£  x * ever  — ^ xc.  — 3 r-  — *r  cc\x  r 
c c c c— — — <v<vcs/cgK-»o»rd  j tiDirc*  X'Xr'-r'cccffC'c— — rsr c -moti r^o: c x «ccr  c— — ccr  ae: r^x  Xi r 3 a »r rv c 

r mmrminini^ini^inirfrmirinirfri^ininimninirininirmss  33  3333  33  jirt*r^ir>krx'inx>x:x»'X x>x:r^r'r^-r^r'f^r^xx:x'vcx:xx  xx- xx  xx  x 

COCCCC  CC  CCC  CCCCC  CCC  C CC  CCCC  CCCC  3C  C CCCCC  C CCC  CCCCCCCCCCCC  CCCCCCCCCCCwCCeCvJCw 

I I I I I I I I • I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

O - CCC  CCC  CC  CCCCC  CCC  CCCC  CCCC  C^C^ViX  C««  fvx  C»«<vx  CJ-5  <VX  C 3CC  P40  C9  X(VX  Otf  C<\X  CJ  X <\X  C^  X <V  X C^<C  (\ 

co<mx  otec^x  o^x  (Nix  cxr  (>X)04(  cvx  000— — rvifgry  ».ioxx»orr««cecC'(r  ooo«4i(M<Mwr,rj  33  xiu.  xxx^f'-exxo 

^ (X3 —— fw  c*nt^)3 3 3Snr>-&-0  -or-  MOCSIMJ1  —— — — — — (\j<vc\jcvk\>  <voj(\,,'\jr\.(\j(v<vrs,rw<M<voj<\ 


- j xi  cr*cca'  c—<vo^xiXf^«5'C  — <vo^x)xr^  XO'C— <vn^»n  x»^  O'c  — 'NrOtf  x»xr^x(rc—  ir  xr^xcr  c— xr~ 

———-i-*—-*—  _^rM<Nj<x'Mrsjr\if'irvj'\jrg»»>o-Ofn>0'0-o>nK>o^  33333  333  3ns) rtf) runsirt >n>n 00  o j ooco 


DIODE  AND  HEADER  TEMPERATURES  AT  5. IS  nSEC  (Continued) 


T 


£ ir  K“.K-KfrK“K,K*riri*“K:irK,K‘»rir 

< c c c cc cc cccc c c c ccc c ccc cccccc 

44444*44444444444444444444 

r cccccccccccccccccccccccccc 
c cc cccccccccccccccccccccccc 
tr  y K-K(fK^rKrirrr^i<-i<-irnrK,rr,rifrrr 


•0  f»"  K"  tnr  tr  ff  tr  tr  Y- Y'.  KW  **  K- lOtO  K“  »0  IT  K*  rr  K*0 

c cccccccccccccccccccccccccc 

4 4444*444444444444444444444 

c cccccccccccccccccccccccccc 
c cccccccccccccccccccccccccc 
k*  r.K  »<■(<“ K-ioirino lOK'inno imo trio KtioiCfncK'. 


c cccccccccccccccccccccccccc 

4 44444444444444444444444444 

c cccccccccccccccccccccccccc 
c cccccccccccccccccccccccccc 
ro  !OlOK)tO»OlOrriOlOK)lO»OlOlO»OiO»OK>lO»QK>»OtO»OlOK) 


•r  ir»r r r. »r k*. io »or r. it k; iririr k; io »o »r k; »r k; »r ir 
c cccccccccccccccccccccccccc 


c cccccccccccccccccccccccccc 
c cccccccccccccccccccccccccc 
X-  •TTlOK,»0»0»OlOfO*OlO»OlOK3K:»OlOlOlO»0»OtOK>!r)fO»0 


*0  fOf0l0lO»r»0»0l0l0t0»0»0»0l0»0K5»0»0»0»0f0»0KJ»0t0»0 

c cccccccccccccccccccccccccc 

4 4444*444444444444444444444 

c cccccccccccccccccccccccccc 
c ccccccecececcoccccccooceec 
r tr  r Y-  nr  |T  r.  r,  »r  »r.  tr.  tr  tr  tr rr  lO  10 IO 10 IO  KJ  k*  »o  tr  tT 


m »r»r  »OfO»ofO»oiOK,,fOto»o»o»o»oio»o»o»o»oto»OK;K)»o»o 

c cccccccccccccccccccccccccc 

4 44444*44444444444444444444 

c cccccccccccccccccccccccccc 
c cccccccccccccccccccccccccc 
*r  »o»r  inrioioiririotr  io»rir  loinoioirioir  lotoinr^io 


c cccccccccccccccccccccccccc 

4 » 4 ♦ 4 4 444444**4444444444444 
c c c cccccccocccccccooccccccc 
c cccccccccccccccccccccccccc 
k r tr  »or,.ior:io»^»oiOK*.»oir»o»rK‘.K,»o»o»o*o»oir»r.»o»o 


O to  lOlOfOlOlOfOiOK)  K>»0  »0  »0»0»0  »OiO»OfO  K)  fO  »OfO  fOfO 
c c cccccccccccccccccccccccc 

4444*44444444444444444444  4 

c cccccccccccc  ccc  cccccccccc 
c c cccc cccc cccccc cccc cccccc 

K K K-K-KKK'KKK.KK'K’rKKKK'KfrK.K-K-K.KK 


•O  lO  fO  fO^POPOlOlO'O'OlOlOfOiOlOfO  K)K)fOfOlO»0*OfOlOlO 

c cccccccccccccccccccccccccc 

♦ *4*4444444*44444*444444444 

c ccccccccocccccccccccccoccc 
c cccccccccccccccccccccccccc 

»0  K,»r*Oi*'fr»OK*iOlO»0»0»^»0»OK*.iOlO»OlO^‘»0»OlO»OlOK‘ 


c *".»o«opo*o»o*o*o»o*o»o»o*OK>*o#o*o»o»OK,)»o*o*o»o*o»r 
c cccccccccccccccccccccccccc 

4 *4*44444444444444444444444 

rv  c\ »\, <\.c\,c\r\. —i— «« occooc 
- :cccccccccc::=cccccccccccc 
k-  r ir  r mr  <r  *■  r r r ir  r ;r r r tr  Y- r.r  ir  fr.r.r 


n -"IfOfOKV'VO  OrOOfOrOKVOK)'OfOfOrO»OfOfOrOK>f^‘0»0 

~ cccccccccccccccccccccccccc 

4 ♦444*444444444*44444444444 

c.:  3 for\  — co  acr^xtr»o<\  — ©ccr^x^'O—c 

<\»  — ~ — — «ccc r ccc 

K"'  «0  »^  K^KV0^*0 ■0*'.»0K)*0  ✓'»0*O*rK'  K*K' <0 


O "0i0f0*<>0’0^)-v0  0«0<0»0<0f00*0'0f0»0,0'0f0i0*0*0 

c cccccccccccccccccccccccccc 

4 4444*4*4444444444444444444 

*r  £ crox.cc  c<\*x  xr^r*  xx io-*cc,rcir,<7'<v*rrk-a:c 
ca:  ocnx  xir  -diorg—ccr  err-  x-^ior^cccr-irio^c 
x x tTiTiTiriOiOir  .»  a ^ 3 33  *o»o»rioK>»o 


■X  ^.■\(\i;>jfMrxcxkvcx.^cx<V-xogv'x.\)rx<\jrxrsjrsjfxt\j 

c c.ccccc csccccccccccccccccc 

» I I I I I I I I I I I I t I I I I I t t t t I I » I 

X c 3 tt<VXCT«<VX>  tn»0  = «3}«V0CtC!VCC 

0 ©csr-4-i<\.rx.>iiOto-»*  3&S)  x»o-o^cxco  ao 
<\  * K'*n»m»o»^»oK>»o»OfO»rn*o*o»o»o»o»orOfO»o*o»o^ 


/ if*  xo'c«- o*»o^x>xr**  «0'C'-<v»o^/)X^<00'c-4 
r*  'tn  t)  0 odd)  c joo*  o-O'O'O'O'O'O'O'j*  o« 


PLOT  OF  MAXIMUM  TEMPERATURE  VERSUS  TIME 


MAY  78  M H CAUSEY  DAAG29-77-G-0190 


UNCLASSIFIED  MSSU-EIRS-EE-78-3  ARO-15304. 1-A-EL  NL 


DIODE  TEMPERATURE  PROFILES 


K)  lOKiinK^^io^oiOKjfO^OfOfOioinioioiOKitOfOKjiooiOfOPOfoioioiOfOiorniOKifOfOiOfoiOfOJOKjiOfOfOfOfOKjiriniOfOfniotOfointoioioiOfntototoiOfoioiotoioK)^ 

c ccccccccccccccccccccccccccccccsccccccccccccccccccccccccccccccccccccccc: ;cc : 

c <\ i-or*-  o<w>f>‘crg»rco-*^r'Oio>oc^r-^vCcino^-Hr>'»oof*‘^^cxor'f^r^coc:(wiO'^or'in^»ni^»rc&r'»o^^^ciro'-<^NDcc^^r‘  «kn®<\j^O'fOvC 
c c c c — ku  ^^iT.ir»r.^xf^f^cccra  c—— tvicx  ^irxr-cc  c-'<\,ir.r.r^cc  cc\  *0x0*  c— — — —.-cc'accr—x.ctr.*  j^Acv-ccr 
»o  io^k;  ioooicionio»r»rroK}K}»on»OK?K:iOKrK>irK>»OK?»o.x*^'3^x*.^x^^tntr>iflir:»nintn«0O'0oo«0r'‘Nr‘*f'-r*-^r*.0vC»0»c.c«0'OO.c>c  ox  sCO«c.r 

*o  oo»o»oo^*oinin»o»oo«o^^^»o»o»o»o*o*OK)io»o*o»oK)io»o»o»oio»o»o»o»o»o»o*o»o»o»o»o»o»o»o»oK)»o*otoK,)io»ofo*o»o»o»oK;fo»o»oK;K5»o»o»o»o^'.*o»o»r»o 
c occcoooccocoocccccccocccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c cNjroir vc<co--<»oji^)<oc<\i^-£co«f01£)0  ^ircc-Hina  »o<cc^ojoc^oNO»o«cof^inir»rirvCflc^^ootoa  ^eo-<K>»oioc\(eco3.»<x©flCtrioer'‘*— r^no  «r 
c c =cccc— «\<v<v<v<x»0'r*o»r^^^^nrirx'Cf^rw*«0‘cc— <\i<v*o^ir»cr*‘XO'^<vfO*n^xao,cccccca>0'0'ro‘*cc*  r*.f*«cxtr.r 
r #r ir^ ,r r.r.r mr.fr.ir .r .*r«nnrr*ir.«nr.r.ir r.inr r.irr.inrrr. 3^34^3^3343 ^trtr.ir ir.tr  ir ir irooxtcxx inr irir  1r1r.tr ir ir tr iTirtr.tr.tr .r 

fO»OOrTfOfOK)iO»000»OK>K)H^O»OrOOlO»OrOK»rOK)»Of>OOK>»00»OlOfOK>»0»OlOK)»OlOfOlOKiK)lOfOK>iOlOtOfOOtOK>fOfOK)rO»00»OrOOOK)fO»OKyOK>fO»0»o 

c :cccccccccccccocccccccccccccccccccccccccccccc;cccccc ccccccccccccccccccccccc 

c <v<virxxirxxh*ccc  cc<v*o^trr-cccc\(^«c«— ir^caroxc^c  30 ir-*f'*x-<cccr^f'-cco «33h^oc  c acr'*tc^»c»*c ^ir.mH«xfrcN3«N 
w ccc  cccccccccc— — — r\*cwcs(<\j<\jtoor»o**trtf,5ir«cxr*'xecac«’-— cv»o*rr^cto*c  cccccccccco  0 j a- a xxaca.r^f-r-  ^ 
< T,o»r*oioK)»^»rto^'Kjtr»r*oioK*»o*r»o»nfO»oK)»o»o»o»o»o*o»o»o»o»nK'»ofr»o»nK,*K'»o^^a^  9 99  99  9j'.*r.T.\r.  tr.ir.tr  ir.trtr.tn^  3333333333  ts 

•O  O^OOK>^^K>-0*0J0«n*0K)»0i0*n'0f0»0»0'0f0»0«0)0KV0O»0'0i0»0*0»0'0K)f0*0»0K>»0»0t0»*T0K>»0»0O‘0f0»0,0»0K)^>(0^-0'0*0»0OK>f0»0*0'0'0«0»0f0»0'-, 
C CCSCCCCCCCCCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCi 

i 3 -00— fNiOjcvoifgfoicr'.^^trtr.c^ccpc— inar»£flC^»oxo'»or^-NX-'r^»ooflC^tr-M^irtr.<vo'tr.c«c<vf*‘ic3'c<N.a'x»rc«c»i»«\rc. 
c cccccoccccccccccccccccccccccc" — — <\<\KMcvfOi033^irxr'r-oco'ccccco,o>o,cocf^r'X'T«Ctrtnir»r 33 333 k- 

tO  fO »0»0 l0»0f0»0i0»0*0»0»0l0»0O»0»0O*0l0l0»0l0l0Ot0K>l0l0»0l0»0l0l0»0t0»0in»0#0»'»0»0»0»0l0»0OK>»0 9 9 999  »OiO»OlO K)»0»0»0»r K>K>*r) lOfOfO'OfO’OfO*’' 


a 3fO'TfOK>ir»OlO'OK)fO»0»0»r)fr)fOfO»OrOK)m*OfOCVt(\(\0(<V<X<V<Ni(N<OJv^«V<XfX<'<<X<V*V''>J<V«\<\01|<V(\JAjfSjCg  (\CSJCgc\J<\<\tC\J<\iv'\CgrtJ(\f^CV.C\iCV  A*  XCVfVCVX  \ 

c-.c;coccccccccccccccccccccccccccccccccccccccccccccccccccccccc  CCCCCCCOCCCwCC 
I t I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I t I I I I t I I I I I I 

3 C r=CSOCCCCCCCCCOCCCCCCCCC3C«M^C3«»WC3CDW-CC3«(V<C3«M4C3(C(M«CC3«^C3«NtC3C(VCC3«'M«  T 
c;  O O3C(VflO3<t(\'0O3<^«£ 03<C(\i^O O ^ ^lT)lTt0'£*0r‘-fv'  OCX <C0  a OCO  ^-^CyjfyjKJK;  *■  ^ cc  CCCO  C 3 

' ^ x—  —csjojrsjfOK)^^  3-^T)>a3xr^r^<o«'tJJ'a'  — — — »-«\(cv<\i<NKV(V<\jcg<\j<vcvPw(\j(v<\<rvj<\i<\jc\joy<<v<\;(\j(vo*K> 


^ — «*<*  <M  <M  <\KM  <V  CJ  <M  <M<\J  (V<\»  CVJ  <VOI  CMAKM  (ViMAKV  <\» 


•«•  ^xr^cca c^o<»o^ir>4Dr*flOO'c— <Mioj3rf)^r.xa>c^<vo^rf)'0r,**a*c^4<NjiO4tn3f,“axrc— <Ni*o^^or‘XO'c»-<MK>x»nsor^eco'c  — c 

aJ'NJOsJTJTJCMOJPJ  AjCVi'O-OrO'O-O’OfOFO^-nx*'  S'*  .3  3-  3-  3 * :3  .Dt/I -HO  tH/l  OfiO  0 00X3. 0 


DIODE  TEMPERATURE  PROFILES  (Continued) 


fj—>|1liHOTin<li— Winn  r ihji 


*o.r*o*o*o*o*o*o-o*o*o*o*o*o*r*o*o*o*o*o*o*o»o*o*o 

cccccccccc:cccccccccccccc 

c *o«occ  ooj»*r.'Cf'f^oin»o-*«krckOO'<\;vr>KcOG 
cr  icr*  «c  -Cir  * *r  cattt'X**rt\.c:*r',f,*o— c 

s 3 s 3 jfrinrtrrttn 


*o  *0*0  *0^5  *o*o  »o»o  *o*o*o  *o*o  *o*o  *o*o  *o*o*o«o*o*o*o 

CCC'CCCCCCCCCCCCCCCCCCCCCC 

c«r.o>ra  *oxa  n*^«^*o»^c>c^n<ain^»o<r»of*'C 
sc<pccf**«c.r**r<\oa,««c**r-<= 
»r  ir*r  *rirxw'-vr*r  a ***^«  *.**»*“  ****r*o*r»or 


O *0*0  *0*0*0  *0*0  *0*0«0*0*0*0  *0(0  *0*0*00*0*0*0*0*0 

c ccccsrcccccccccccccccccc 
•*•■»♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
^ CvC  cvr-*oc  <vr^  ^irac-*^«cccc— *r^xr-co  c 
<£  ‘./if  .* -/rvocv.o.-.cca  ^r-r-vorf-^oo*— sc 
» -fas  s +s  t ^ ^*r*r*o*o*o»o*o*oKvn*o*o 


**  00-">00*00*0*0«0*0000*00*0 *0*0*00*0*0*0 
w cscccccccccccccccccccccc 

L OC  £*f  n.cc/'.K  c 

£ *o  *o  *o  io  r *o  n,  o*  <v  n.  <v<m<v  n< — — • «*»  ~ •*  c.  c c c 
>o  *o*o  *o  *o  io  r»*o  *o*o  *0*0*0  *o*o*o  *o*oo*o*o*o*n*o*o 


c\  x*vo*cwo*«vn*f^n*o»fV<N(nio*OiO*o* 

- cc -c::c-cccco:c ;cccccccl 

• f i i i l I i i I I I « I I I I l I I l l I l • 
£ CO.-C  r »«<V^C*tN«CC*tfVCC^t(ViC 


►*r*p*.0  *3  O 1.0  *:  Ot343  03'3'7'0'J'M'0'>'J>00 


PLOT  OF  DIODE  TEMPERATURE  PROFILES 


PLOT  OF  IMPURITY  PROFILE 


PLOT  OF  INITIAL  AND  FINAL  ELECTRIC  FIELD  PROFILES 


1 


PLOT  OF  DIODE  VOLTAGE  TRANSIENT 


I 


<\ 

Ci 

<v 

<\ 

<v  <\ 

<\  <v 

«_ 

c 

c 

c 

c 

c c 

c c 

c 

♦ 

♦ 

4 

♦ 

♦ ♦ 

4 4 

c 

c 

c 

c c 

c c 

c 

c 

X 

c 

x 

c 

If  c 

IT  C 

c 

c 

• 

• • 

X 

x 

a 

a 

ir 

x 

<v  <\ 

c 

c 

1 

» 

4 - - 

— - 4- - 

- - 4 - - - 

• - 4 •- 

---  4 4 

•>  4 ••  *■ 

4 If 

1 

c 

a 

t 

1 x 

If 

1 

C 

a i 

e x 

I 

c 

<rc 

i 

X 

1 

c 

CK 

1 x 

1 

1 c 

1 

I 

c 

CO. 

» » 

4 C 

1 

c 

err 

1 C 

IT 

I 

c 

C ff 

! x! 

1 

c 

c or 

i 

* 

1 

c 

t a 

1 x 

l C 

1 

1 

c 

C ff 

1 J 

1 

c 

c a 

1 ^ 

1 

o 

U X 

1 » 

1 

c 

C X 

1 

X 

1 

c 

C X 

1 x 

1 C 

c 

1 1 

O' 

4 O 

<c 

1 

o 

U X 

1 « 

» 

•0 

I 

c 

C X 

1 *5 

1 

c 

C X 

1 

X 

1 

c 

1 X 

t 

1 c 

J 

1 

o 

U X 

1 1 

♦ 

4 r. 

I 

c 

C X 

i 

I 

iC 

1 

c 

C X 

i x 

►- 

f 

1 

c 

C X 

i 

S£ 

x 

I 

c 

C X 

I x 

0 

1 

1 C 

a 

c 

O X 

1 1 

L. 

c 

<J  X 

4 © 

> 

• 

1 

c 

O X 

1 • 

if 

1 <v 

U l 

c 

C X 

• 

c 

1 

< 

c 

C X 

1 

o 

O X 

1 9 

o 

c 

<J  X 

i i 

> 

c 

(J  X 

♦ 

• <V 

u 

c 

u x 

• • 

c 

» fc 

c 

c 

C X 

1 

« 

c 

o 

U X 

1 

1 X 

c 

vJ  X 

1 C 

1 

c 

o X 

4 ^ 

• x 

1 

c 

t • 

1 

o 

u X 

1 

x 

1 

1 

c 

u x 

1 

1 X 

1 

£ 

U X 

* c 

♦ 

£ 

U X 

• 

o 

o a 

1 • 

•* 

c 

O X 

» — 

o 

o r 

1 

1 

U X 

1 

» r- 

— 

O C X 

1 c 

J" 

x 

♦ 

1 

2 u a 

l A 

1 X 

O' 

I 

1 

3 Cl 

1 • 

1 x 

Q XU 

1 

1 

1 

2 X 

w 

1 

1 

1 

Q X 

u 

1 

1 _ 

* 

4 —4 

c 

O* 

04 

eg 

CM 

rc  CJ 

<V  <\ 

c 

C 

c 

© 

c c 

o e 

C 

♦ 

♦ 

♦ 

4 4 

© 

© 

© 

© 

© e 

c c* 

>T 

c 

x> 

C 

X>  c 

IT  C 

C 

• 

• 

• 

f. 

S 

•o 

fO 

M <M 

*1  m 

o 

-fiv- 


B 


APPENDIX  B 
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Subroutine  DEPV 
Subroutine  DHTEMP 
Subroutine  DHT2D 
Subroutine  DHT2D1 
Subroutine  DHT2D2 
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SUBROUTINE  BANDA6 ( LRMAX , NCD * NHD » L XROW » LAROW • I ACOL.A. 
*X.NBNOTX»PIVMlN> 

SUB.  »BANDA6»  EVALUATES  THE  SOLUTION  OF  A BANDED 
(OR  DIAGONAL)  SYSTEM  OF  (LRMAX)  LINEAR  SIMULTANFOUS 

HP^iasarww  ?ivsii*w?gs:,iN  tfchn,*if- 

h5S85  : gMo^fiSraiTS?  ?5c»..g{K8WINSSl!SLV 

- c8LfiMNRNUMRtHE^F  CENTRAL  DIAGONAL  IN  *A» 

- COLUMN  NUMBER  FOR  EQUATION  CONSTANTS  IN  »A* 

- COEFFICIENT  AND  CONSTANT  ARRAY 

nbndtx  - n8mbe&°ofA  ^ollTlONS  RESFRVED  FOR  BND  VALUES  AT  TOP 
OF  »X»  ARRAY 

- NUMBER  OF  COEFFICIENT  diagonals  BELOW  central 
(OR  MAIN)  DIAGONAL 

- NUMBER  OF  COEFFICIENT  DIAGONALS  ABOVE  CENTRAL 
(OR  MAIN)  DIAGONAL 

caw : mmm*,  CONSTANT  ARRAY  ROW  DIMENSION 
LACOL  - COEFFICIENT  AND  CONSTANT  ARRAY  COLUMN  DIMFNSION 

DIMENSION  X(LXROW) » A ( LAROW  * LACOL ) 


NCD 
NHD 
A 


NDB 

NDA 


ROWS  ( 2 ) — ( NDB ) 


NDBSNCD-I 
NDA=NHD-NCD-1 
PIVMIN=1.0E35 

UPPER  TRIANGULATTON  OF 
IF (NDB.EQ. 1 ) GO  TO  20 
DO  5 LR=2»  NOB 

bo°10NE?oiL^QF * NDB 

ALFA=A ( LR  # LCO) /A ( LR+LCO-NCD, NCD ) 

DO  15  LC=lfNDA 

A ( LR , LC+LCO ) =A { LR . LC+LCO ) -ALFA* A ( LR+LCO-NCD , LC+NCD ) 
CONTINUE 

A ( LR . NHD ) =A ( LR » NHD ) -ALFA*A ( LR+LCO-NCO » NHD ) 

CONTINUE 

CONTINUE 

CONTINUE 


UPPER  TRIANGULARTZATION  OF  ROWS 
DO  25  LR=NCD. LRMAX 


(NDB+1 )- (LRMAX) 


SEARCH  FOR  MINIMUM  DIAGONAL  PIVOT  ELEMENT 
TA=A (LR-NDB'NCD) 

IF(DABS(TA).LT.DABS(PIVMIN) ) PI VMlNrTA 
DO  30  LC0=1»N0B 

ALFA=A ( LR » LCO ) /A ( LR+LCO-NCD  # NCD ) 

DO  32  LC=1#NDA 

A ( LR . LC+LCO > =A ( LR , LC+LCO ) -ALFA* A ( LR+LCO-NCD  ,1  C+NCD ) 
CONTINUE 

A ( LR . NHO ) =A ( LR  * NHD ) -ALFA*A ( LR+LCO-NCD » NHO ) 

CONTINUE 

CONTINUE 

BACK  SUBSTITUTION  FOR  X ( LRMAX )-X(LRMAX-NDA+l) 
LRSTOP=LRMAX-NDA+l 
~ R=LRMAX*LRSTOP,-1 


.b»NHD) 


DO  45 
TAsA (L 

KKSLR+NBNDTX 
KSTOP=LRMAX-LR 
IF (KSTOP.EQ.O)  GO  TO  40 
00  35  KsltKSTOP 
TAsTA-X  I K+KK ) *A  ( LR  * K4-NCD ) 

CONTINUE 

CONTINUE 

X ( KK ) =TA/A ( LR » NCD ) 

CONTINUE 

BRi^-LfiMAXINDAI0N  F°R  ><,LRHA)(•N0A,•  X<NDB*1) 
DO  55  LRSLRST » 1 # -1 


SUBROUTINE  BANDA 6 (Continued) 


B2 


80 

81 

82 

83 

84 

85 

86 
7 


8 


50 

55 

C 


KKzLR+NBNDTX 

TA=A(LR»NHD) 

DO  50  Ksl.NOA 

TA=TA-X(K+KK)*A(LR»K+NCD) 

CONTINUE 

X ( KK ) =T  A/ A ( LR  r NCD ) 
CONTINUE 

RETURN 

END 
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SUBROUTINE  BKDEPL 


1 


P<1). BKDEPL 


SUBROUT  I NE  BKDEPL  ( NND . XDEPL  * XMFT  • XDEPEP , XDEPl  )»  XI.  OX.  DOPl  , OOPEP . 
ADOPU , ALF  A TD . EMA X , AERMAX . EMAXL » EM A XU . E » T * ALFA  T N , VDFP  , 

AITCMAX. IDB1. IDB2. IDB3.NP1.SDOPL.S0OPEP.SDOPU) 


rEGRAL  I 
[Rations 


ARE  PFRFORMEn. 


SUCH  THAT  THE 
UAL  TO  ONE. 
UNTIL  THF 
UM  NUMBER  OF 


NND 

XDEPL 

XMFT 

XDEPEP 

XL 

Bopl 

DOPEP 

DOPU 

ALFATD 


- NUMBER  OF  NODES  ALONG  AXIS  OF  OIODF 
“ B5ELETION  RE6I0N  LOWER  BOUNOARY 
- POSITION  of  metallurgical  JUNCTION 


position  of  metallurgi 

EPITAXIAL  LAYER  BOUNDARY 
DIODE  LENGTH 

NODE  SPACING  ALONG  DlOOF  AXIS 
LOWER  DOPING  CONC. 

EPITAXIAL  DOPING  CONC. 

UPPER  DOPING  CONC. 


DOPU  - UPPER  DOPING  CONC. 

ALFATD  - AVALANCHE  I0NI7ATI0N  COEFFICIENT  TFMPERATURF  DFPENDENCF 

EMAX  - MA  xim5mfdepl  ETION  REGION  ELECTRIC  FTELD 
AERMAX  - MAXIMUM  AVALANCHE  BREAKDOWND  INTEGRAL  ERROR 
EMAXL  - LOWER  BOUND  FOR  EMAX 
EMAXU  - UPPER  BOUND  FOR  EMAX 

¥ : 

ALFAIN  - AVALANCHE  BREAKDOWN  INTEGRAL 
VDEP  - DEPLETION  REGION  VOLTAGE 
ITCMAX  - MAXIMUM  NUMBER  OF  ITERATIONS  ON  EMAX 

IDB1  - DEBUG  SENTINEL  FOR  SUB  TNALFA 

IDB2  - DEBUG  SENTINEL  FOR  SUB  RKDEPL 

IDB3  - DEBUG  SENTINEL  FOR  SUB  EPROF 

NP1  - JUNCTION  ORIENTATION*  l - NP.  0 - PN 


SDOPL  - LOW 


- DEBUG  SENTINEL  FOR  SUB  RKDEPL 

- DEBUG  SENTINEL  FOR  SUB  FPROF 

- JUNCTION  ORIENTATION*  l - NP. 


SOT : 


ON  ORIENTATION*  l 

depletion  region 

IAL  DEPLFTION  REG 


- NP* 


- PN 


Ton  region  space  charge 

PLETION  REGION  SPACE  CHaRGF 
ION  REGION  SPACE  CHARGE 


DEBUG  NAMELIST 


NAMEL I ST  /DEBUG/NND » XDEPL . XMET . XDEPEP  * XDEPU . XL » DX . DOPL  * DOPEP  * DOPU . 
AALFATD* ALFAIN»EMAX. AERMAX. EMAXL* FMAXU* VDFP. VDEPL.VDEPEP.VOEPU. 


«ITCMAX*NP1 

DIMENSION  E(NND) .T(NND) 

INTEGER  FLAGEB 
INITIALIZATION  COMPUTATIONS 

SET  THE  MAXIMUM  ELECTRIC  FIELD  ITERATION  COUNTER 
ITC=0 

INITIALIZE  ELECTRIC  FIELD  SEARCH  BOUNDARIES 

TEMAXUSEMAXU 

TEMAXL=EMAXL 


200  CONTINUE 


BEGIN  ITERATION  ON  MAXIMUM  ELECTRIC  FIELD 

}tS-?t§  l ITeRATt0N  counter 

CALCULATE  TRIAL  EMAX 
EMAX: ( TEMAXU+TEMAXL ) /2 . n 

GENERATE  NEW  ELECTRIC  FIELD  PROFILE  FOR  EMAX 
CALL  EPROF ( NND . OX . XDEPL . XMET . XDEPEP . XDEPU . XL . SDOPL . 

I SDOPEP » SDOPU » EMAX .E.IDB3* FLAGEB) 

EVALUATE  THE  AVALANCHE  IONIZATION  INTEGRAL  FOR  THE  NEW  ELECTRIC 
FIELD  PROFILE 

CALL  INALF A ( NND . XDEPL  * XMET ♦ XDEPFP. XDEPU, OX  * DOPL . OOPEP . 

IOOPU, ALF ATD . E * T * ALFA IN  * TOB1 . NPl I 


SDOPL. 


AOOPU. 


k TD. E.T. ALFAIN* TDB1. NPl) 


jHECK  FOR  A DEPLETED  BULK  RE6I0N.  SUB  F 
•GUAL  TO  ONE  FOR  A DEPLETED  BULK  REGION. 
tF(FLAGEB.LE.O)  GO  TO  300 


FPROF  SETS  »FLAGFB» 


A BULK  REGION  HAS  BEEN  DEPLETED.  IF  ALFAIN  IS  GRFATER  THAN  ONE 


C 

C 


c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

s 


c 

c 

c 


c 

c 


SUBROUTINE  BKDEPL  (Continued) 

A VALID  SOLUTION  MAY  EXIST.  PERFORM  ANOTHER  ITFRATION J 
IF(ALFAIN.GT.l.O)  GO  TO  300 

BAD  DIODE  DESIGN.  THE  REQUIRED  EMAX  CAN  NOT  BE  ACHIEVED 
WITHOUT  DEPLETING  A BULK  RE^iON.  TERMINATE  XQT. 

WRITF(6» 290 ) FLAGEB. ALFA IN, EMAX 
290  F0RMAT(/,40H  A BULK  REGION  HAS  RFEN  DEPLFTED  WITHOUT. 

A44H  SATISFYING  THE  AVALANCHF  BRFAKDOWN  INTEGRAL./. 

A32H  XQT  IS  TERMINATED  IN  SUB  BKDEPL ./. 

A8H  FLAGEB: . 12. 5X . 7HALFAIN=»E 1 0 . 3. 5X. 5HEMAXS.F10 . 3) 

STOP 

300  CONTINUE 

UPDATE  ELECTRIC  FIELD  SFARCH  BOUNDARIES 
IF(ALFAIN.GT.I,0)  GO  TO  3) 0 

INCREASE  ELECTRIC  FIELD  TRIAL  EMAX 

TEMAXL=EMAX 

GO  TO  320 

310  CONTINUE 

DECREASE  ELECTRIC  FIELD  TRIAL  VALUE 
TEMAXU=EMAX 

320  CONTINUE 

EVALUATE  AVALANCHE  I0NI7ATI0N  INTEGRAL  ERROR 
IF(AlFAlN.GE.l.O)  AERROR=ALFA IN-1 .0 
iF(ALFAIN.LT.l.n)  AFRRORrl . 0/ALFA IN-1 . 0 

DEBUG  LIST  OPTION 

IF ( IDB2.EQ. 1 ) WRITE(fe»505)  ITC*TFMAXL»TFMAXU. 

AEMAX . ALFA I N » AERROR  * XDEPL  # XDEPU 
505  FORMAT (I10.7E10.3) 

HAS  THE  AVALANCHE  IONIZATION  INTEGRAL  BFEN  DFTERMINFD  TO  THF 
AF^o8.Ee?2e^MAX)  GO  TO  600 

THE  REQUIRED  ACCURACY  HaS  NOT  BFFN  ACHIEVED.  IF  THF  MAXIMUM 
NUMBFR  OF  ITERATIONS  HAS  NOT  BEEN  PFRFORMFD.  MAKE  ANOTHER! 
IF(ITC.LT.ITCMAX)  GO  TO  200 

THE  MAXIMUM  NUMBER  OF  ITERATIONS  HAS  BEEN  PFRFORMFD  WITHOUT 

' " iitf  < 


ACHIEVING  THE  SPECIFIED  ACCURACY. 
AND  TERMINATE  XQT. 

WRITFJ6, 51 0 ) EMAX. ALFAIN, AERROR 
FORMAT  f " " - 


WR] 


AN  FRROR  mfssagf 


510  FORMAT (//.45H  THE  MAXIMUM  NUMBER  OF  ITERATIONS  ON  FMAX  HAS. 
A48H  BEEN  PERFORMED  WITHOUT  ACHlFVlNG  THF  SPFCIFIED. 

A27H  ACCURACY  FOR  THE  AVALANCHE./. 

A55H  BREAKDOWN  INTEGRAL.  EITHER  THF  ACCURACY  »A^RMAX«.  THE. 
A41H  NUMBER  OF  ITERATIONS  *ITCMAX».  OR  ONE  OF, 

«19H  THE  ELECTRIC  FIELD./, 

A49H  BOUNDARIES  »EMAXL  OR  EMAXU* , SHOULD  BF  MODIFIFD. »/, 
A/.6H  EMAX=.E10,3.5X.7HAl  FAlNr »E 1 0 . 3.SX , 7HAERR0R=,F1 0 . 3) 

600  CONTINUE 

EVALUATE  DEPLETION  VOLTaGF.  VDE 
CALL  DEPV(SDOPL»SDOPEP.SDOPU.XD 
AVDEPL . VDEPEP . VDEPU , VDEP ) 


FPL . XMET . XDEPFP , XDFPU , 


C1000 

i 


WRITE<6. DEBUG) 


RETURN 

END 


RETURN 

DEBUG  OUTPUT  FROM  St  IB  BKDEPL) 


B5 


SUBROUTINE  DEPV 


*TEMP{ 1 ) .DEPV 
1 
2 

3 

4 

5 

6 
7 


6 

if 


L6 

17 

18 

19 

20 
21 
22 

23 

24 

23 

27 

II 

30 

31 

33 

34 

35 

36 

II 

n 

41 

42 

43 

44 

33 

47 

43 

49 

50 


C 

C 


C 

C 

8 

C 

i 

C 

8 

C 


C 

c 


c 

c 

c 

c 

c 

c 

c 


SUBROUT I NE  DEPV  ( SOOPL » SDOPEP . SDOPIJ . XDEPl. # XMF V # XDEPFP , XDFPIJ  # 
4VDEPL  » VDEPEP# VDFPU# VDEP' 

SUB  DEPV  EVALUATFS  ANALYTICALLY  THE  VOLTAGE  ACROSS  THE  LOWER# 
- - ’TION  F 


EPITAXIAL#  AND  UPPER  DEPLETION  RFGIONS. 
VARIABLE  DEFINITIONS: 

DOPlJ  - UPPER  DEPLETION  REGION  SPACE  CHARGF 


HARGF 


GION  ROUNDARY 


C 

C 


C 

8 

C 


XDEPL  - LOWER  DEPLETION  REGION 

XMET  - METALLURGICAL  JUNCTION 

XDEPFP  - EPITAXIAL  LAYER  BOUNDARY 

XDEPU  - UPPER  DEPLETION  REGION  BOUNDARY 

VDEPEP  - EPlfAX?ALLDE^ETl6NTREGTONTv8LTAGE 

VDEPU  - UPPER  DEPLETION  REGION  VOLTAGE 

VDEP  - TOTAL  DEPLETION  REGION  JUNCTION  VOLTAGF 

CONSTANT  ASSIGNMENTS 
DATA  Q/1.6E-19/ 

DATA  PERM/1. 06E-1 2/ 

PERM  - FARADS/CM 

EVALUATE  LOWER  OFPLETION  REGION  VOLTAGE 
VDEPL=0/2. 0/PERM*SDOPL* ( XMET-XDFPL ) **2 

DOES  UPPER  DEPLETION  REGION  RND  LIE  WITHIN  EPITAXIAL  REG. 
IF(XDEPU.LE.XDEPEP)  GO  TO  100 

UPPER  DEP.  REG.  RND  LIES  OUTSIDF  EPITAXIAL  RFG. 
VDEPEP=Q/PERM* ( SDOPL* ( XMET-XDEPL ) -SDOPEP/2 . 0* ( XDEPEP-XMFT  > ) * 
* ( XDEPEP-XMET ) 

VDEPU=Q/2. 0/PERM*SDOPU* ( XDEPU-XDFPEP) **? 

GO  TO  200 

100  CONTINUE 

DEPLETION  REG.  RND  LIES  WITHIN  EPITAXIAL  REG. 

VDEPFP=Q/2 . O/PERM+SDOPEP* ( XDFPU-XMET ) **2 
VDEPU=0.0 

200  CONTINUE 

EVALUATE  TOTAL  DEPLETION  REGION  VOLTAGE 
VDEP=VDEPL+VDEPEP+VDEPU 

RETURN 

END 
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SUBROUNTINE  DHTEMP 


♦TEMP ( 1 ) 
1 
2 

3 

4 

5 

6 

7 

8 

,8 

U 


.DHTEMP 


It 

36 

37 

38 

39 

40 

i *»1 

4§ 

44 

45 

46 


SUBROUTINE  DHTEMP ( IFLA6 . UNO . NNH . XDH . THKH » HSPFC  * HDFN » DTIMF  # 

ATEMPD#  TEMPH) 

a quasi-two-dimensional  thermal  model  composfd  of  »nnh»  uncouplfd 


ONE  DIMENSIONAL  MODELS  IS  FORMULATED.  THERMAL  CONDUCTION  OCCURS 
PERPENDICULAR  TO  THE  DIODE  AXIS  THROUGH  A PARALLEL 

combination  OF  one  dimensional  conductors. 


VARIABLE 
IFLAG  - 


E DEFINITIONS 

- NEW  SIMULATION  SENTINEL  WHICH  PROVIDES  FOR  SUCCESSIVI 
SIMULATIONS  IN  ONE  RUN  STREAM.  IFLAG  IS  RESET  IN  MA 
PROGRAM  FOR  NEW  SIMULATION  AND  SET  IN  SUR  DHTEMP  DUR 
FIRST  REFERENCE 

NUMBER  OF  NODE  POINTS  Al  ONG  DIODE  AXIS 

number  of  nooe  points  along  y-axis  (Perpfndicular  to 

DIODE  AXIS)  WITH  IN  SUBSTRATE  MATERIAL 
HEADER  MATERIAL  THICKNFSS 
HEADER  MATERIAL  THERMAL  CONDUCTIVITY 
HEADER  MATERIAL  SPECIFIC.  HEAT 
HEADER  MATERIAL  DENSITY 


ROVIDES  FOR  SUCCESSIVE 
IFLAG  IS  RESET  IN  MAIN 
SET  IN  SUR  DHTEMP  DURING 


XDH 

THKH 

HSPEC  - 
HDEN 

□TIME  - 
TEMPO  - 


TEMPH 


VHffi  : riKiAiTfUmi^MPER.TURES  ALONG  DIODF  AXIS.  VALUES 
SUPPLIED  BY  CALLING  PROGRAM  AND  USFO  AS  RND  BY  SUB 

TEMPH  - 9C5edimfnsional  HEADER  TEMPERATURE  ARRAY  FOR  RFTURNING 
HEApER  TEMPERATURES  TO  CALLING  PROGRAM 

NBNDTX  - COLUMN  POINTER  FOR  LOCATING  THF  RESPECTIVE  ONE 

DIMENSIONAL  HEADER  TEMPFRATURE  PROFILES  WITH  IN  THF 
TWO  DIMFNSIONAL  HEADER  TEMPERATURE  ARRAY  FOR  SUB 
BANDA6.  SUB  BANDA6  SOLVES  THE  INDIVIDUAL  HEADER 
TEMPERATURE  PROFILES  AND  STORES  THFM  IN  ARRAY  TEMPH 

DIMENSION  TEMPDt 101) > TEMPH (12» 101 ) .A<12.4) 


THE  MAIN  PROGRAM  SETS  IFLAG  EQUAL  TO  ZERO  AT  THF  RFGINNING  OF  EACH 
NEW  SIMULATION  TO  GUE  INITIALIZATION  CALCULATIONS.  IF  IFLAG  FfiUAL 
ONE  SKIP  INITIALIZATION  CALCULATIONS 
IF ( IFLAG. EG. 1 ) GO  TO  25 


INITIALIZE  SUBSTRATE 
DO  20  Krl.NNH 
DO  15  J=l»NNO 
TEMPH(K» J)=300.0 
15  CONTINUE 
20  CONTINUE 

NNHM1=NNH-1 

NNDM1=NND-1 

DXHSXDH/NNH 

DXDXH=DXH*DXH 

AA=THKH/HSPEC/HDEN 


TEMP  ARRAY 


SET  FLAG  TO  SKIP 
IFLAG=1 


INITIALIZATION  CALCULATIONS 


25  CONTINUE 

INTERMEDIATE  COEFFICIENT  EVALUATION 
AlrDTlME*AA 

A2=- ( DXDXH+2 . 0*DTIME*AA ) 


A3=A1 

A4S-DXDXH 

INITIALAZE 

NBNDTXrO 


TEMPERATURE  COLUMN  POINTER 


evaluation  OF  »NND*  one  dimensional  temperature  PROFILES 
DO  50  J=1»NND 


COEFFICIENT  EVALUATION 

DO  55  Nsl.NNHMl 

A(N»1 )=A1 

A(N.2)=A2 

A(N#3)=A3 

A(N#4)=A4*TEMPH(N. J) 

55  CONTINUE 

INCLUDE  BOUNDARY  CONDITIONS 
A ( 1 . 4 ) =A ( 1 » 4 ) -A  1 * TEMPO ( J ) 


vO  CD  03  CD  CD  CD  CD  CD®  OD  ® 
OvOtt'JO'UI-t  OJfV) *— * O 


T 


SUBROUTINE  DHTEMP  (Continued) 


C 

c 

c 

c 


A ( NNHM1 , 4 ) =A ( NNHM 1*4) -300  3 0* A ( NNHM1 » 3 ) 

EVALUATE  ONE  DIMFNSIONAL  TEMPERATURE  PROFILF 

CALL  BANDA6 ( NNHM1 . 2 . 4 . 1? 12  # 12  * 4 . A . TFMPH . NRNDTX * PI VMIN ) 

INCRFMENT  TEMPERATURE  COLUMN  POINTER 
NBNDTX=NBNDTX+12 
50  CONTINUE 


RETURN 

END 


B8 


SUBROUTINE  DHT2D 


♦TEMP ( 1 ) . DHTPD 


SUBROUT I NE  DHT20 ( IFLAG , NND . NNH , XL » XDH . THKHX , THKHY . HSPEC . 
ftHDEN.  DT I ME . TEMPO . TEMPH . SMAX . DTHMA  X , I TPRH  * IBND  * S ) 

A TWO  DIMENSIONAL  THERMAL  MODEL  FOR  SIMULATING  THERMAL  CONDUCTION 

cON8u|¥lVl?Y?ISpf‘clFlCTHlA^  AND^D^NSI  TY^O^ThI? * SURS^RaVf^MATER I AL 
ARE  ASSUMED  TEMPeRATORE  INDEPENDENT.  THERE  IS  NO  JOULE  DISSIPATION 
WITH  IN  THE  SUBSTRATE.  THE  TEMPeRATURF  PROFILES  FOR  THE  NEXT 
POINT  IN  TIME  (TIME+DTIME)  ARE  CALCULATED  EACH  TIMF  THE 
SUBROUTINE  IS  CALLED.  AN  iterative  SOLUTION  procedure  similar 
to  the  sor  technique  is  used  to  evaluate  thf  new  temperatures 


PROFILES. 

declarative  statemints 

DIMENSION  TEMPD(lOl) ,TEMPH(1?*101 ) t TEMPS ( 12*101) 
DIMENSION  TEMPS l ( 12. 101) .AA( 12.4) 


INTEGER  S.SMAX 
VARIABLE  DEFINITIONS 

IFLAG  - NEW  SIMULATION  SENTINEL  WHICH  PERMITS  SUCCESSIVE 

SIMULATIONS  in  one  run  stream  by  queing  simulation  initialization 

COMPUTATIONS. 

NND  - NUMBER  OF  GRID  OR  NODE  POINTS  ALONG  DIODE  AXIS  (X-AXIS) 

NNH  - NUMBER  OF  GRID  OR  NODE  POINTS  ALONG  Y-AXIS  (PFRPENDICULAR 
TO  DIODE  AXIS)  WITH  IN  HEADER  MATERIAL. 

XL  - DIODE  LENGTH.  (X-AXIS) 

XDH  - HEADER  MATERIAL  ThICKNFSS 

THKHY  - HEADER  MATERIAL  THERMAL  CONDUCTIVITY  ALONG  Y-AXIS 
(ASSUMED  CONSTANT) 

THKHX  - HEADER  MATERIAL  THERMAL  CONDUCTIVITY  ALONG  X-AXIS 
(ASSUMED  CONSTANT).  NORMALLY  THKHXrTHKHY . THKHXrn.fl  FOR 
QUASI-TWO-DIMENSIONAL  HFADER  THFRMAL  MODEL.  FOR  THIS  RFDUCED 
MODEL  ONLY  ONE  ITERATION  IS  REQUIRED.  INORnFR  TO  OVERIDE 
The  CONVERGENCE  CHECK  FOR  THIS  CASE  SET  SMAXrl  AND  DTHMAX-100 . 0 

HSPEC  - HEADER  MATERIAL  SPECIFIC  HEAT  (ASSUMED  CONSTANT) 

(DEN  - HEADER  MATERIAL  DENSITY  (ASSUMED  CONSTANT) 

dtimf  - time  step  increment  FOR  NEXT  time  step 

TEMPO  - LINEAR  ARRAY  FOR  TEMPERATURES  ALONG  DIODE  AXIS.  VALUES 
SUPPLIED  BY  CALLING  PROGRAM  AND  USED  AS  BND  RY  SUB  DHT20 


TEMPH  - TWO  DIMENSIONAL  HEADER  TEMPERATURE  ARRAY  FOR  RETURNING 
HEAOFR  TEMPERATURES  TO  CALLING  PROGRAM.  DURING  ITERATIONS  TEMPH 

represents  Temperatures  at  previous  point  in  time  rut  is 

REDEFINED  AS  NEW  TEMPERATURE  VALUES  BEFORE  RETURNING 
TO  THE  CALLING  PROGRAM 

SMAX  - MAXIMUM  NUMBER  OF  ITERATIONS  FOR  HEADER  TEMPERATURE. 

pTHMAX  - CONVERGENCE  CRITERIA  FOR  HEADER  TEMPERATURES  ALGORITHM. 
MAXIMUM  ACCEPTABLE  NORMALIZED  CHANGE  IN  temperature  between 
iterations  s AND  S+l. 

ITPRh  - header  temperature  iteration  print  sentinil.  WHEN 
ITPRhsI  s ano  dthsmx  arf  printed  for  EACH  ITfRAT! 

I BND  - INDICATES  THE  TYPE  BND  AT  XrO  AND  XrXL.  0 - CONSTANT 
TEMPERATURE.  1 - BLOCKING 

TEMPS  - TWO  DIMENSIONAL  HEADER  TEMPERATURES  FROM  PRFVIOUS 
ITERATION  (S). 

TEMPS1  - TWO  DIMENSIONAI  HEADER  TEMPERATURES  FOR  NEW 
ITERATION  (S+l). 

AA  - COEFFICIENT  ARRAY  FOR  ONE  DIMENSIONAL  TFMPFRATURE  PROFILES. 
NND  - NUMBER  OF  TEMPERATURE  NODES  ALONG  DIODE  AXIS 


ITERATIONS  TEMPH 
E RUT  IS 
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SUBROUTINE  DHT2D  (Continued) 


C 

C 

c 

i 


c 

c 

c 

c 

c 

8 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c. 

c 

c 

c 


ONE  DIMENSIONAL 

perature  array 


NNH  - NUMBER  OF  TEMPERATURE  NODFS  ALONG  Y-AXIS  (PERPENDICULAR  TO 
DIODE  AXIS)  IN  SUBSTRATE.  MAXIMUM  VALUF  OF  NNHMX 

NNHMAX  - NNH  MAXIMUM  VALUE  WHICH  IS  THE  SAME  AS  THE  *ROW  DIMENSION* 
FOR  ARRAYS  TEMPH,  TEMPS.  AND  TEMPS1. 

NNHBND  - COLUMN  POINTER  FOR  LOCATING  THE  RESPECTIVE 

temperature  profiles  with  in  the  two  dimensional  tfmpei 
tempsi  for  subroutine  bandag  which  solves  the  individual  one 

DIMENSIONAL  TEMPERATURE  PROFILES. 

S - ITERATION  COUNTER 

fiWWWfiB  i?*N!?ER»T?8Ns!R  TEHPER,T,me  *T  N00E  ,J'n 

DTHSMX  - MAXIMUM  NORMALIZED  CHANGED  IN  HEADFR  TEMPERATURE 
BETWEEN  THE  S AND  S+l  ITERATIONS. 


JTPRDH  - header  temperature  ITERATION  PRINT  SINTINEL. 
ITPRoHsI  S ANO  DTHSMX  ARE  PRINTED  FOR  EACH  ITERATION 


WHEN 


100  CONTINUE 


C 

C 

C 

C 


SIMULATION  INITIALIZATION  CALCULATIONS 

IFLAG=0  IMPLIES  FIRST  TlMF  STEP  FOR  NEW  SIMULATION  AND  THAT 
SIMULATION  INITIALIZATION  CALCULATIONS  SHOULD  BE  PERFORMED. 

IF ( IFLAG.NE • 0 ) GO  TO  200 

SET  IFLAG  TO  INDICATE  THAT  THE  SIMULATION  HAS  BFEN  INITIALIZED 
IFLAG=1 

INITIALIZE  TEMPERATURE  ARRAYS 

BO  110  1=1.101 
0 110  J=1 . 12 
TEMPH ( J. I ) =300.0 
temps (J. I ) =3oo , o 
TEMPSK  J.  I )=300.0 
110  CONTINUE 


c 

c 


c 

c 

8. 


c 

8 

C 


EVALUATE  GRID  INCREMENTS 
DX=XL/(NND-1) 

THE  NUMBER  OF  HEADER  GRID  POINTS  ALONG  Y-AXIS  IS  (NNH+1 ) 
DYSXDH/NNH 

EVALUATE  MODIFIED  GRID  COUNTS 
NN0M1=NND-1 

NHR3M3 

NNHM?=NNH-2 

INITIALIZE  COLUMN  POINTER  INCREMENT 
NNHMAX— 12 

*20o"cont7nue 

PRE-ITERATION  computations 


partial  coefficient  evaluation  which  will  accommodate  a varilable 
time  step  controlled  BY  the  calling  program 
a=-Thkhy*dtime*dx*6x 

B=-2 . 0*A<fDY*DY*DX*DX*HDEN*HSPEC 

dT=dx*dx*dy*oy*hden*hspec 
op=thkhx«-dy*dy*otime 


C INITIALIZE  ITERATION  COUNTER 


C 

C 


initialize  tempsi  with  temph  from  previous  time  stfp 


88  I!  8 '!{:"$ 

TEMPSK  J.I)=TEMPH(  J.I) 


BIO 


SUBROUTINE  DHT2D  (Continued) 


160 

161 

162 

163 

164 

165 

166 

III 

169 

H? 

172 

m 

}?i 

m 

us 


161 

162 

183 

184 

186 

ill 

189 

190 

m 

i?j 

ill 

197 

198 

H2 

201 

202 

203 


ffl 

211 

115 

gig 

ill 

220 


>1 


223 

225 

228 

231 

li! 

238 

239 


C- 

C- 

C 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

o 

c- 

c 

c 

c 

c 

c 


210  CONTINUE 


300  CONTINUE 


C 

C 


C 

C 

C 

C 

C 

C 


S. 

C- 


iteration  loop 

TRANSFER  TEMPS1  TO  TEMPS  SO  THAT  NEW  TEMPS1  VALUE  MAY  BF  CALCULATFO 
DO  305  1=1. NND 
OO  305  jrl.NNH 
TEMPS( J. I )=TEMPS1 (J.I) 

305  CONTINUE 

INCRFMENT  ITERATION  COUNTER 
S=S+1 

EVALUATE  NND  ONE  DIMENSIONAL  TEMPERATURF  PROFILFS  ALONG  Y-AXIS 
WITH  C0UPLIN6  WITH  ADJACENT  COLUMNS. 

INITIALIZE  COLUMN  POINTfR  FOR  THF  FIRST  ONE  DIMFNSIONAL 
TEMPFRATURE  PROFILE. 

MNHBnD=NNHMAX 

400~CONT INUe” 

~ ~ m — — — jr " i ~r r r t — — — — i m ni  — 1 ■ ai  — i 1- it  M — 1 

IF  BLOCKING  BND  SPECIFIED  FOR  X=n,  EVALUATE  TEMPSl(J.l)  AND 
IRfgtflB.l)  Go  TO  500 

COEFFICIENT  EVALUATION  FOR  TEMPSKJ.2)  WITH  BLOCKING  BND 

DO  410  J=1 .NNHMl 

AA ( J. 1 ) =A 

AA ( J. 2) SB 

AA(J,31=C 

AA( J.4)=D1*TEMPH(J.?)+D?*(TEMPS(J.3)-TEMPS( J.2) ) 

410  CONTINUE 


INCLUDE  BND  FOR  TEMPSl(J.P) 

AA ( 1 , 4 ) =AA  C 1 . 4 ) -A*TEMPD ( 2 ) 

A A ( NNHMl . 4 ) =AA ( NNHMl .4) -C *300.0 


>4. 1212. 12.4. AA» TEMPS 1 .NNHBND.PIVMIN) 


C 

C 


C 

C 


C 

C 

C 

C 


EVALUATE  TEMPSltJ.J 
CALL  BANDA6 ( NNHMl .i 

INCREMENT  COLUMN  POINTER  FOR  NEXT  TEMPERATURF  PROFILF 
NNHBND=NNHBND+NNHMA  X 

EOUATE  TEMPS1  (J.2)  AND  TEMPSUJ.1)  FOR  BLOCKING  BND 

DO  420  J=1 . NNH 

TEMPS1(J.1)=TEMPS1(J.2) 

420  CONTINUE 

500~CONT INUE 

IF  CONSTANT  TEMPERATURE  BND  SPECIFIED  FOR  X=0.  EVALUATE  TEMPSKJ.2) 
IF ( IRND.NE . 0 ) GO  TO  600 

COEFFICIENT  EVALUATION  FOR  TEMPSUJ.2)  AND  CONSTANT  TEMP  BND 
DO  510  J=l. NNHMl 
AA( J. 1)=A 
AA ( J.2 i =B 
AA ( J. 3) =C  _ 

AA ( J.4 ) =D1*TEMPH ( J.2 ) +02* ( 300 ,0+TEMPS( J. 3)— 2 .0*TEMPS ( J.2 ) ) 

510  CONTINUE 


INCLIlOE  BND 
AA(1.4)=AA( 

AA ( NNHMl • 4 ) =AA ( NNHMl .4) 


AA ( 1 • 4 ) =AA (1.4) -A*TEMPD ( 2 ) 

-C*300.0 


EVALUATE 
CALL  BANI 


TEMPS1 (J.2) 

I0A6( NNHMl . 2. 4. 1212. 12. 4. A A, TEMPS 1. NNHBND.PIVMIN) 


INCRFMENT  COLUMN  POINTER  FOR  NEXT  TEMPERATURF  PROFILF 
NNHBND=NNHBND+NNHMAX 


Bll 


SUBROUTINE  DHT2D  (Continued) 


600  CONTINUE 


EVALUATE  TEMPERATURE  PROFILES  FOR  TEMPS 1 ( J# 3)  - TFMPS1 ( J»NNDM2) 
DO  620  I=3»NNDM2 

EVALUATE  COEFFICIENTS  FOR  TEMPSl (J, I) 

DO  610  jrl.NNHMl 
AA< J,l)rA 
AA(J.2)=B 
AA( J,3)rC 

AA<J»4)rDl*TEMPH(J»I)+D2*(TEMPSC J, 1-1 )+TEMPS< J, 1*1  » 
«-2.0*TEMPS<J*I) ) 

610  CONTINUE 

INCLUDE  0ND  FOR  TEMPSl(J.I) 

AA  < 1 . 4 ) =AA (1,4) -A*TEMPD ( I ) 

A A ( NNHM1 , 4)rAA (NNHM1 ,4) -C *300.0 

EVALUATE  TEMPSl(J.I) 

CALL  BAN0A6 ( NNHM 1, 2, 4, 1212, 12, 4, AA, TEMPS1 , NNM8ND , PI VMIN ) 

INCRFMENT  COLUMN  POINTER  FOR  NEXT  TEMPFRATURF  PROFILE 
NNHBNDSNNHBND+NNWMAX 

620  CONTINUE 
700  CONTINUE 

~ t ~ 1 m ■ t i — ■ M M M — — * »>>**■!! 

IF  BLOCKING  BND  SPECIFIFD  EVALUATE  TEMPS1 ( J, NNDM1 ) AND 
II??rn^n£?i!  GO  TO  BOO 

COEFFICIENT  EVALUATION  FOR  TEMPSl (J, NNDMl ) KITH  BLOCKING  BND 
DO  710  Jr 1 , NNHM1 
AA( J,l)rA 
AA(J,2)rB 


AA(J,3)rC 

A A ( J , 4 ) rD 1 *TEMPH ( J > NNDM 1 ) +D2* ( TFMPS ( J , NNDM2 ) -TEMPS ( J , NNDM 1 ) ) 

710  CONTINUE 

INCLUDE  BND 

AA ( 1 , 4)rAA ( 1 » 4) -A*TEMPD (NNDM1 ) 

AA(NNHM1,4)=AA(NNHM1 ,4)-C*300.0 

EVALUATE  TEMPSl (J, NNDMl) 

CALL  BANDA6 ( NNHM1 ,2,4, 1212, 12,4, AA,TEMPS1 ,NNHBND,PIVMIN) 

EQUATE  TEMPS l ( J, NNDM 1)  AND  TEMPS1 ( J.NND)  FOR  BLOCKING  BND 
DO  720  Jr] , NNH 

TEMPS 1(J,NN0)=TFMPS1 (J.NNDMt) 

720  CONTINUE 

’fl00"cONTlNUE*" 

IF  CONSTANT  TEMPERATURE  BND  SPECIFIED  FOR  XsXL,  EVALUATF 
TEMPSl ( J,NNDM1 ) 

IF(IBND.NE.O)  GO  TO  900 

COEFFICIENT  EVALUATION  FOR  TEMPSl (J, NNDM 1 ) AND  CONSTANT  TEMP  BND 

DO  810  Jri.NNHMl 

AA( J,l)rA 

AA(J,2)=B 

AA( J,3)rC 

AA ( J , 4 ) rD 1 *TEMPH ( J , NNDM1 ) *D2*  < TFMPS ( J , NNDM2 ) +300 . 0 
«-2.0*TEMPS< J.NNOM1) > 

810  CONTINUE 

INCLUDE  BND  FOR  TEMPSl ( J , NNDMl ) 

AA ( 1 , 4 ) r A A ( 1 , 4 ) -A*TEMPD ( NNDMl ) 

A A ( NNHM1 , 4 ) rAA ( NNHM1 , 4) -C *300,0 

EVALUATE  TEMPSl (J, NNDMl) 

CALL  BAN0A6(NNHM1 ,2,4, 1212, 12*4, AA, TEMPSl ,NNHRNO,PTVMlN) 
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SUBROUTINE  DHT2D  (Continued) 


C- 

c 

c 


900  CONTINUE 


C 

c 

c 

s 


CONVERGENCE  CHECK 

evaluate  maximum  normalized  change  in  temperature  eor  this 
iteration 

Bthsmx=o.o 

0 9]0  I=2,NNDM1 

8?h1=(  TeM^I? J?  1 > -TEMPS < J» I ) ) /TEMPS 1 1 J # I ) 

IF(DTHS.GT.DTHSMX)  dthsmx=dths 
910  CONTINUE 

PRINT  ITERATION  SUMMARY 
IF(ITPRH.EQ.l)  WRITE (6# 920 ) S»QTHSMX 
920  FORMAT ( 3H  S=. I2.RX.7HDTHSMX=.E8.3> 

IF  DTHSMX  .LE.  OTHMaX  THE  SPECIFIED  CONVERGENCE  HAS  BEEN  ACHIEVED 
IF(OTHSMX.LE.DTHMAX)  GO  TO  1000 

CONVERGENCE  HAS  NOT  BEEN  ACHIEVED.  IF  THE  ITERATION  COUNT 

is  less  than  thf  maximum  specified  start  another  iteration 

IF ( S.LT.SMAX ) GO  TO  300 


930 


C 

C 

8. 

8 

c 


the  maximum  number  of  iteration  has  been  performed  without 

ACHIEVING  CONVERGENCE.  WRITE  ERROR  MESSAGE  AND  TERMINATE  X«T 
WRITE (6.930)  S. DTHSMX 

FORMATI///, 46H  *****  CONVERGENCE  FAILURE  IN  SUBROUTINE  DHT2D, 

A/»27H  *****  EXECUTION  TERMINATED./. 

A9H  *****  S=.I2,5X.7HDTHSMX=»E8.3./. 

A34H  *****  INCREASE  SMAX  AND/OR  DTHMAX ) 

STOP  ************************************************************** 


1000  CONTINUE 


C 

C 


transfer  new  substrate  temperature  valufs  to  array  temph 

DO  1010  1=1. NND 
DO  1010  J= 1 » NNH 
TEMPH(J,I)=TEMPS1 (J.I) 

1010  CONTINUE 

RETURN  TO  CALLING  PROGRAM  WITH  NEW  HEADER  TEMPERATURES 
RETURN 
END 
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SUBROUTINE  DHT2D1 


♦TEMP ( 1 ) .DHT2D1 
1 
2 

3 

4 

5 

6 


3 

56 

57 

II 

II 

63 

II 

67 

68 

69 

70 

71 


r 4 

i 


! 


C 

c 


c 

c 


c 

c 


c 

c 

£ 

c 

c 

s 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  DHT2D1 ( I FLAG • NND.NNH. XL  * XDH*  THKHX » THKHY.HSPEC . 

AHDEN , DT I ME . TEMPO . TEMPH , SMAX . OTHMAX . I TPRH . IBNO  * S ) 

A TWO  DIMENSIONAL  THERMAL  MODEL  FOR  SIMULATING  THERMAL  CONDUCTION 

through  the  diode  substrate  or  hfadfr  Material.  Thf  Thfrmal 

with  in  the  substrate,  the  tempfrature  profiles  for  thf  next 
EP.Ii^IN  time  (TIME+DTIME)  .are  calculated  each  timf  the 
SUBROUTINE  IS  CALLED.  AN  ITERATIVE  SOLUTION  PROCEDURE  SIMILAR 

to  the  sor  technique  is  used  to  fValuaTf  thf  new  Temperatures 


echnique  is  used  to  fvaluatf  the  nfw  temperature* 

0hT20I  DIFFERS  FORM  DHT2D  IN  THAT  THF  VERY  LATEST 
IN  THE  ITFRATIVE  PROCEEDURF. 


PROFILES.  0HT2DI  DlFFE 
VALUES  OF  TEMP  ARF  USED 


BlMENS^O^fEMPDi^oiVI^EMPHt  12#  101 

DIMFNSION  TEMPS 1 (12.101) .AA(12»4) 


) .TEMPS( 12.101) 


INTEGER  S.SMAX 

^^EN^Fi«®2¥?ON  SENTINEL  WHICH  PFRMITS  SUCCFSSIVF 
SIMULATIONS  IN  ONE  RUN  STREAM  BY  QUEING  SIMULATION  INITIALIZATION 
COMPUTATIONS. 

NND  - NUMBER  OF  GRID  OR  NODE  POINTS  ALONG  DIODE  AXIS  (X-AXIS) 

NNH  - NUMBER  OF  GRID  OR  NODE  POINTS  ALONG  Y-AXIS  (PFRPENDICULAR 
TO  DIODE  AXIS)  WITH  IN  HEADER  MATERIAL. 

XL  - OIODF  LENGTH.  (X-AXIS) 

XDH  - HEADER  MATERIAL  THICKNFSS 

THKHY  - HEADER  MATERIAL  THERMAL  CONDUCTIVITY  ALONG  Y-AXIS 
(ASSUMED  CONSTANT) 

THKHX  - HEADER  MATERIAL  THERMAL  CONDUCTIVITY  ALONG  X-AXIS 

''  LLY  THKHXrTHKHY.  THKHXsO.O  FOR 


(ASSUMED  CONSTANT). 
QUASI-TWO-  " 

Mor  ' ■ 

TH! 


NORMAL 


ASI-TwO-DIMENSIONAL  HFADER  THFRMAL  MODEL.  FOR  THIS  RFDl 
DEL  ONLY  ONE  ITERATION  IS  REQUIRED.  INORDFR  TO  OVFRIDE 
E CONVERGENCE  CHECK  FOR  THIS  CASF  SET  SMAXrl  AND  DTHMAX: 


RFDUCFD 
1IDE 

rHMAXSIOO.O 

HSPEC  - HEADER  MATERIAL  SPECIFIC  HEAT  (ASSUMED  CONSTANT) 

HDEN  - HEADER  MATERIAL  DENSITY  (ASSUMED  CONSTANT) 

DTIMF  - time  step  increment  FOR  NEXT  TIME  STFP 

TEMPO  - LINEAR  ARRAY  FOR  TEMPERATURES  ALONG  DIODE  AXIS.  VALUES 
SUPPLIED  RY  CALLING  PROGRAM  AND  USED  AS  BNO  RY  SUR  DHT2D 

TEMPH  - TWO  DIMENSIONAL  HEADER  TEMPFRATURE  ARRAY  FOR  RETURNING 
HEADFR  TEMPERATURES  TO  CALLING  PROGRAM.  DURING  ITERATIONS  TEMPH 
reprfsents  temperatures  AT  PREVIOUS  POINT  IN  TIME  BUT  IS 
REDEFINED  AS  NEW  TEMPERATURE  VALUES  BEFORF  RFTURNING 
TO  THE  CALLING  PROGRAM 

SMAX  - MAXIMUM  NUMBER  OF  ITERATIONS  FOR  HEADFR  TEMPFRATURE. 

DTHMAX  - CONVERGENCE  CRITERIA  FOR  HEADER  TEMPERATURFS  ALGORITHM. 
MAXIMUM  ACCEPTABLE  NORMALIZED  CHANGE  IN  TEMPFRATURE  BETWFEN 
ITERATIONS  S AND  SM. 

ITPRH  - header  temperature  iteration  print  sfntinfl.  when 

ITPRHSl  S AND  DTHSMX  ARF  PRINTED  FOR  EACH  ITFRATION 


IBNO  - INDICATES  THE  TYPE  BND  AT  XrO  AND  X=XL. 
TEMPERATURE.  1 - BLOCKING 


0 - CONSTANT 


TEMPS  - TWO  DIMENSIONAL  HFADFR  TFMPERATURES  FROM  PRFVIOUS 
ITERATION  (S). 

TEMPS1  - TWO  DIMENSIONAL  HEADER  TEMPERATURES  FOR  NFW 
ITERATION  (S-fl). 

AA  - COEFFICIENT  ARRAY  FOR  ONE  DIMENSIONAL  TFMPFRATURE  PROFILES. 
NND  - NUMBER  OF  TEMPERATURE  NODFS  ALONG  DIODF  AXIS 


B14 


SUBROUTINE  DHT2D1  (Continued) 


NNH  - NUMBER  of  TEMPERATURE  NODFS  ALONG  Y-AXTS  < PERPENDICULAR  TO 
DIODF  AXIS)  IN  SUBSTRATF.  MAXIMUM  VALUF  OF  NNHMX 

NNHMAX  - NNH  MAXIMUM  VALUF  WHICH  IS  THE  SAMF  AS  THF  ‘ROW  DIMENSION* 
FOR  ARRAYS  TEMPH,  TEMPS.  AND  TEMPS1. 

NNHBND  - COLUMN  POINTER  FOR  LOCATING  THF 


THF  RESPECTIVE  ONE  DIMENSIONAL 
DIMFNSIONAL  tfmperaturf  ARRAY 
LVFS  THF  INDIVIDUAL  ONF 


tempfRature  PROFILES  WITH  IN  THF  TWO  

TEMPS1  FOR  SUBROUTINE  BANDA6  WHICH  SOLVFS 
DIMENSIONAL  TEMPERATURE  PROFILES. 

S - ITERATION  COUNTER 

DTHS  - NORMALIZED  CHANGF  IN  HEADFR  TEMPFRATURE  AT  NODE  CJ.I) 
BETWFEN  THE  S AND  S+l  ITERATIONS. 

DTHSMX  - MAXIMUM  NORMALIZED  CHANGED  IN  HEADFR  TFMPERATURF 
BETWEEN  THE  S AND  S+l  ITERATIONS. 


ITPRDH  - header  temperature  ITERATION  PRINT  SINTINFL, 
ITPRDHsI  S AND  DTHSMX  ARE  PRINTED  FOR  EACH  ITERATION 


WHEN 


100  CONTINUE 


w - 

cc 


i. 


simulation  initialization  calculations 

iflag=o  implies  first  ttmf  step  for  new  simulation  and  that 
simulation  initialization  calculations  should  bf  pfrformfd. 

IF(IFLAG.NE.O)  GO  TO  200 

set  iflag  to  indicate  that  the  simulation  has  bffn  initialized 

IFLAG=1 

INITIALIZE  TEMPERATURE  ARRAYS 

do  no  in. 101 
DO  lto  J=1 • 12 

TEMPH(J.I)=300.0 
TEMPS(J.I)=300.0 
TEMPS1 (J»I)= 300.0 
110  CONTINUE 

evaluate  grid  increments 

DX=XL/(NNO-l) 

THE  NUMBER  OF  HEADER  GRID  POINTS  ALONG  Y-AXIS  IS  (NNH4-1) 
DY=XOH/NNH 

EVALUATE  MODIFIED  GRID  COUNTS 

NNDM1=NND-1 

NNDM2=NND-2 

NNHM1SNNH-1 

NNHM?=NNH-2 

INITIALIZE  COLUMN  POINTER  INCREMENT 
NNHMAXS12 

"20o"cONTINUE  " 

■ ■»  • «■»  • « • * a • a ■ a « • « • mm  w a « ■■  v mm  ■*  • • • tm  • mm  • mm  • mm  * • « g.  • tm  • mm  w 

pre-iteration  computations 

PARTIAL  COEFFICIENT  EVALUATION  WHICH  WILL  ACCOMMODATE  A VARILABLE 

time  step  CONTROLLED  BY  THE  CALLING  PROGRAM 

a=-thkhy*dtime*dx*dx 

B=-2 . 0*A+DY*OY*OX*DX*HDFN*HSPEC 

C=A 

Dl=DX«DX*DY*OY*HDEN*HSPFC 

D2=THKHX*DY*DY*DTIME 

INITIALIZE  ITERATION  COUNTER 
5=0 

INITIALIZE  TEMPS 1 WITH  TEMPH  FROM  PREVIOUS  TIME  STFP 
DO  210  in.NNO 
DO  210  J=1.NNH 


SUBROUTINE  DHT2D1  (Continued) 
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TEMPSl  ( J#  I ) =TEMPH ( J . I ) 
210  CONTINUE 


300  CONTINUE 


ITERATION  LOOP 

TRANSFER  TEMPS 1 TO  TEMPS  SO  THAT  NEW  TEMPS1  VALUE  MAY  BF  CALCULATFD 
DO  305  I=1»NND 
DO  305  JS1.NNH 


TEMPS(J»I)=TEMPS1 (J.I) 

CONTINUE 

INCRFMENT  ITERATION  COUNTER 
S=S+1 

EVALUATE  NND  ONE  DIMENSIONAL  TEMPERATURF  PROFILFS  ALONG  Y-AXIS 
WITH  COUPLING  WITH  ADJACENT  COLUMNS. 

INITIALIZE  COLUMN  POINTFR  FOR  THF  FIRST  ONE  DIMFNSIONAL 
TEMPERATURE  PROFILE. 

NNHRnD=NNHMAX 


400  CONTINUE 


IF  BLOCKING  BND  SPECIFlFD  FOR  X=n,  EVALUATE  TEMPSl (J,l)  AND 
ii^lBN^'.NE.l)  GO  TO  500 

coefficient  evaluation  FOR  TFMPSKJ.2)  WITH  RLOCKING  BND 

DO  410  J=1 * NNHMl 

AA< J.1)=A 

AA(J.2)=B 

AA( J.3)=C 

A A ( J , 4 ) sDl *TEMPH ( J » 2 ) +D2* I TEMPS ( J * 3 ) -TEMPS ( J • 2 ) ) 

410  CONTINUE 

INCUIOE  BNO  FOR  TEMPS1(J»2> 

AA ( 1 , 4 ) =A A (1,4) -A*TFMP0 ( 2 ) 

A A ( NNHM1 , 4 ) =AA ( NNHM1 ,4) -C *300.0 


EVALUATE  TEMPS1 ( J,2) 
CALL  BANDA6(NNHM1 .2,4, 


1212, 12,4, AA.TEMPS1 ,NNHBND,PIVMlN) 


INCRFMENT  COLUMN  POINTER  FOR  NEXT  TEMPERATURF  PROFILE 
NNHBND=NNHBND+NNHMAX 

EQUATE  TEMPS  l ( J,  2 ) AND  TEMPSl (J.1)  FOR  BLOCKING  BND 


DO  420  J=1,NNH 

TEMPS 1 ( J, 1 )STEMPS1 ( J,2) 

CONTINUE 


500  CONTINUE 


IF  CONSTANT  TEMPERATURE  BND 
IF( IrND.NE.O)  GO  TO  600 


SPECIFIED  FOR  X=(l»  FVALUAtE  TEMPSl tJ»2> 


COEFFICIENT  EVALUATION  FOR  TEMPS1(J,2)  AND  CONSTANT  TEMP  BND 

DO  510  j:1 , NNHMl 

AA( J, 1)SA 

AA(J,2)=B 

AA(J,3  SC 

AA(  J,4»sDl*TEMPH(J,2U02*(300.0+TEMPS(J,3)-2.0*TEMPS<J,2»  ) 


mm 

510  CONTINUE 


INCLUDE  BND 

AA ( 1 , 4 ) =AA (1,4) -AWTEMPD ( 2 ) 

AA ( NNHMl , 4 ) =A A ( NNHM 1,4) -C*300 . 0 

EVALUATE  TEMPSl 1 J,2) 

CALL  BANDA6 ( NNHMl ,2,4, 1212, 12, 4, AA, TEMPSl , NNMBND , PI VMlN ) 

INCRFMENT  COLUMN  POINTER  FOR  NEXT  TEMPERATURE  PROFILE 
NNHBND=NNHRND+NNHMAX 
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SUBROUTINE  DHT2D1  (Continued) 


C 

C- 

c- 

c 

c 

c 

c 


600  CONTINUE 


C 

C 


C 

C 

C 

C 


C 

C- 


C 

8 


C 

C 


c 

c 


c 

c- 

c- 

c 

c 

c 

c 

c 


evaluate  temperature  profiles  for  temps i c j# 3)  - temps k j * nnom2 ) 

DO  620  1 = 3*  NNDM2 

EVALUATE  COEFFICIENTS  FOR  TEMPSK  J.  I) 

DO  610  JZ1.NNHM1 
AA( J* 1 ) =A 
AA(J.2)=B 
AA(J.3)=C 

AA( J.4)=D1*TEMPH( J, I ) +D2* ( TEMPS! < J # I - 1 ) +TEMPS ( J » I + 1 ) 

A-2.0*TEMPS( J.I) ) 

610  CONTINUE 

INCLUDE  BND  FOR  TEMPSK  J.I) 

AA ( 1 . 4 ) =AA ( 1 . 4 ) -A*TEMPD ( I ) 

AA(NNHMl ,4>=AA(NNHM1»4)-C*300.0 

VALUATE  TEMPSl(J.I) 

ALL  BANDA6 ( NNHM1 *2»4*1212»12*4.AA* TEMPS 1 * NNHRND  *PIVMlN) 

INCRFMENT  COLUMN  POINTER  FOR  NFXT  TEMPERATURF  PROFILE 
NNHBND=NNHBND+NNWMAX 

620  CONTINUE 

*700~CONTTNUE  " " 

• v mm  mm  tm  m w mm  mm  • mm  mm  ^ mm  • « ^ mm  ^ mm  mm  ^ mm  ^ mm  ^ ^ ^ ^ **  « ^ • mm  ^ •“  ^ «»  • mm  mm  • a*  ■■  «•  w ••  • ■■ 

IF  Bl  OCKING  BND  SPECIFIFD  EVALUATE  TEMPS1 ( J*NNDM1 ) AND 
TEMPSK  J.NND) 

IF ( IRND.NE . 1 ) GO  TO  800 

COEFFICIENT  EVALUATION  FOR  TEMPS1 ( J.NNDM1 ) WITH  BLOCKING  BND 
DO  710  jri.NNHMl 
AA ( J* 1 ) =A 
AA( J.2)=B 

AA ( J ! 4 1 =8 1 *TEMPH ( J . NNDM 1 ) +D2* ( TFMPS l ( J * NNDM? ) -TFMPS ( J » NNDM1 ) ) 

710  CONTINUE 

INCLUDE  BND 

AA ( 1 . 4 ) =AA (1*4) - A*TFMPD ( NNDM 1 > 

AA ( NNHM1 ,4)=AA(NNHM1 *4)-C*300.0 

EVALUATE  TEMPSK J.NNDM1) 

CALL  BAND A6 ( NNHM 1.2. 4*1212. 12*4. AA* TEMPS1 . NNHBNO . PI VMl N ) 

EQUATE  TEMPS l (J, NNDM 1 ) AND  TEMPS 1 (J.NND)  FOR  BLOCKING  BND 
DO  720  J=1*NNH 

TEMPSK  J.NNO)=TFMPSl  ( J.NNOM1 ) 

720  CONTINUE 

800  CONTINUE  _ — 

IF  CONSTANT  TEMPERATURE  BND  SPECIFIED  FOR  XrXL.  EVALUATF 
TEMPSK  J.NNOM1) 

IF( IRND.NE. 0)  GO  TO  900 

COEFFICIENT  EVALUATION  FOR  TEMPS1 ( J*NNDM1 ) AND  CONSTANT  TEMP  BND 
DO  810  JSl.NNHMl 


C 

C 


DO  810  J= 

AA ( J* 1 )=A 

aaIj’§  -? 

A A ( J * 4 ) =5 1 *TEMPH ( J » NNDM 1 )+D2*(TEMPSl( J*NNDM2)+300.n 
8“2.0*TEMPS( J.NNOMl ) ) 

810  CONTINUE 

INCLUDE  BND  FOR  TEMPSK .J.NNDM1 ) 

AA ( 1 » 4)=AA ( 1 * 4 ) -A*TFMPD ( NNDM1 ) 

A A (NNHM1 » 4 ) =AA ( NNHM1 . 4 ) -C*300 . 0 


C 

C 


EVALUATE  TEMPSK  J.  NNDM  1 ) 

CALL  BAN0A6 ( NNHMt  »2*4»1212»12.4*AA*  TEMPS 1 » NNHBNO * PI VMlN) 
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ill 


C 

C" 

c- 

c 

c 

c 

c 

c 


900  CONTINUE 


CONVERGENCE  CHECK 

EVALUATE  MAXIMUM  NORMALIZED  CHANGE  IN  TEMPERATURE  EOR  THIS 

iteration 
dthsmx=o . 0 

00  910  I=2»NNDM1 


DO  910  J=1,NNHM1 

feAKg??4W4ii7TB¥5iWii«(TFHPS" '■ J 


.1) 


C 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c- 

c 

c 

c 


c 

c 


910  CONTINUE 

PRINT  ITERATION  SUMMARY 
IF(ITPRH.EQ.l)  WRITE(6#920)  S*DTHSMX 
920  FORMAT ( 3H  S=# I2*SX»7H0THSMX=,E8.3> 

IF  DTHSMX  .LE.  OTHMAX  THE  SPECIFIED  CONVERGENCE  HAS  BEEN  ACHIEVFD 
IE(DTHSMX.LE.DTHMAX)  GO  TO  1000 

CONVERGENCE  HAS  NOT  BEEN  ACHIEVED.  IF  THE  ITERATION  COUNT 
IS  LESS  THAN  THE  MAXIMUM  SPECIFIED  START  ANOTHER  ITERATION 
IF(S.LT.SMAX)  GO  TO  300 

I^t^JA8u?0W^.,TE81U0^SSI8ii§.^B«RS{5S?r«T 

WRITE (6*930 ) S, DTHSMX 

930  FORMAT!///, 48H  *****  CONVERGENCE  FAILURE  IN  SUBROUTINE  DHT2D, 

A/»27H  *****  EXECUTION  TERMINATED* /• 

«9H  *****  S=*I2#SX*7hOTHSMX=#E8.3#/» 

«34H  *****  INCREASE  SMAX  AND/OR  DTHMAX) 

STOP  ************************************************************** 
i000“cONTTNUE 

TRANSFER  NEW  SUBSTRATE  TEMPERATURE  VALUFS  TO  ARRAY  TEMPH 
DO  lOlQ  1=1, NND 
DO  1010  Jsl.NNH 
TEMPH ( J . I 1 =TEMPS l ( J , I) 

1010  CONTINUE 

RETURN  TO  CALLING  PROGRAM  WITH  NFW  HEADER  TFMPERATURES 
RETURN 
END 


"I O'  O' CT'C'OO'O' U»U1  Ul  WUKfl CW1 Ul UW?  -P -C  -P-l?  ^-P-C-C -CO*  OlUU  1^0101040*0  l'orol\>f'Jf'or\>r\>«'OfV>»-»~^t->fc««-»v->k~  »-•*-• 
voa)'40'a>-coN>>-ovoao'><a'(ji^o*\jfc-o>oo)'vJCMJvtoM^cxocS'Jcrcn-co4Mk-'Ov£K»'jacn^04rv)k-o>£)CD-viO'  t*.P04f'o*-ovOCD~gaui-P04rv>*-ovr)a>^icr'cn-eo4  ru^- 
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SUBROUTINE  DHT2D2 


♦TEMP ( 1 ) .DHT2D2 


C 

C 

C 

C 

C 

C 

C 

C 

C 

8 

C 

C 

c 

c 


c 

c 

8 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

8 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUT I NE  DHT2D? ( I FL AG . NND . NNH , XL  * XDH . THKHX , THKH Y , HSPEC . 

AHDEN . DT I MF » TEMPI) , TEMPH , SMAX . DTHMaX  , I TPRH . I BND  . S ) 

A TWO  DIMENSIONAL  THERMAL  MODEL  FOR  SIMULATING  THERMAL  CONDUCTION 
THROUGH  THE  DIODF  SUBSTRATE  OR  HFADFR  MATERIAL.  THF  THERMAL 
CONDUCTIVITY * SPECIFIC  HEAT  AND  DENSITY  OF  THE  SURSTRATF  MATERIAL 
ARE  ASSUMED  TEMPERATURE  INDEPENDENT,  THERE  IS  NO  JOULE  DISSIPATION 
WITH  IN  THE  SUBSTRATE.  THE  TEMPfRATURF  PROFILES  FOR  THF  NEXT 


POINT  IN  TIME 

it  I 


__  . _ _ _ . . _ PROFILES  FOR 

...  (TIME+OTIME)  ARE  CALCULATfD  EACH  TIMF  THE 

SUBROUTINE  IS  CALLED.  AN  ITERATIVE  SOLUTION  PROCEDURE  SIMILAR 
TO  THE  SOR  TECHNIQUE  IS  USED  TO  EVALUATF  THE  NEW  TEMPERATURE 
PROFILES.  FOR  0HT202  THE  HEADER  TEMPERATURF  I 
ROW  WISE  (ALONG  X-AXIS)  RATHER  THAN  COLUMN  WTS 
THE  VERY  LATEST  VALUE!-" 

PROCEDURE. 


i)  rather  than  column  wise  asAinAdh?2di . 

OF  TEMP  ARE  USEO  IN  THE  ITERATIVE 


DECLARATIVE  STATEMINTS 
DIMENSION  TEMPD U 01) .TEMPH (12 
DIMENSION  TEMPS1(12.101).AA(1 


12*101 ) • TEMPS ( 12*101) 

01.4) .TEMP(lOl) 

INTEGER  S.SMAX 
VARIABLE  DEFINITIONS 

IFLAG  - NEW  SIMULATION  SFNTINEL  WHICH  PERMITS  SUrCFSSIVF 
SIMULATIONS  IN  ONE  RUN  STREAM  BY  QUEING  SIMULATION  TNITI ALI7ATION 
COMPUTATIONS. 

NND  - NUMBER  OF  GRID  OR  NODE  POINTS  ALONG  DIODE  AXIS  (X-AXIS) 

NNH  - NUMBER  OF  GRID  OR  NODE  POINTS  ALONG  Y-AXIS  (PERPENDICULAR 
TO  DIODE  AXIS)  WITH  IN  HEADER  MATERIAL. 

XL  - DIODE  LENGTH.  (X-AXIS) 

XDH  - HEADER  MATERIAL  THICKNESS 

THKHY  - HEADER  MATERIAL  THERMAL  CONDUCTIVITY  ALONG  Y-AXIS 
(ASSUMED  CONSTANT) 

THKHX  - HEADER  MATERIAL  THERMAL  CONDUCTIVITY  ALONG  X-AXIS 
( ASSUMED _ CONSTANT ) . NORMALLY  THKHXrTHKHY.  ThKHX=0 . 0 FOR 


QUASI-TwO-OIMENSIONAL  HEADER  thermal  model 
MODEL  ONLY  ONE  ITERATION  IS  REQUIRED. 

The  CONVERGENCE  CHECK  FOR  THI  — - 


urn  FOR  THIS  RFDUCFD 
- INORdfR  TO  overide 

CASE  SET  SMAXrl  AND  DTHMAXslOO.O 


HSPEC  - HEADER  MATERIAL  SPECIFIC  HEAT  (ASSUMED  CONSTANT) 

HDEN  - HEADER  MATERIAL  DENSITY  (ASSUMED  CONSTANT) 

DTIMF  - TIME  STEP  INCREMENT  FOR  NEXT  TIME  STFP 

te;mpd  - linear  array  for  temperatures  ALONG  DIOOE  AXIS.  VALUFS 
SUPPLIED  BY  CALLING  PROGRAM  AND  USED  AS  BND  RY  SUR  DHT2D 

TEMPH  - TWO  DIMENSIONAL  HFADFR  TFMPERATURF  ARRAY  FOR  RETURNING 
HEADER  TFMPFRATIJRES  TO  CALLING  PROGRAM.  DURING  ITERATIONS  TEMPH 
REPRESENTS  TEMPERATURES  AT  PREVIOUS  POINT  IN  TIME  BUT  IS 
REDEFINED  AS  NEW  TEMPERATURE  VALUES  BEFORE  RETURNING 
TO  THE  CALLING  PROGRAM 

SMAX  - MAXIMUM  NUMBER  OF  ITERATIONS  FOR  HEADFR  TEMPFRATURE. 

DTHMAX  - CONVERGENCE  CRITERIA  FOR  HEADER  TEMPERATURFS  ALGORITHM. 
MAXIMUM  ACCEPTABLE  NORMALIZED  CHANGF  IN  TEMPFRATURF  BETWEEN 
ITERATIONS  S AND  S+l. 


PRINT  SFNTINFL.  WHFN 
EACH  ITFRATION 


0 - CONSTANT 


ITPRh  - HEADER  TFMPERATURF  ITERATION 
ITPRh=1  S ANO  DTHSMX  ARE  PRINTED  FOR 

I BND  - INDICATES  THF  TYPE  BND  AT  X=0  AND  X=Xl . 

TEMPFRATURE*  l - BLOCKING 

TEMPS  - TWO  DIMENSIONAL  HFAOf R TFMPFRATURES  FROM  PRFVIOUS 
ITERATION  (S). 

TEMPS1  - TWO  DIMFNSIONAL  HEADER  TEMPERATURES  FOR  NFW 
ITERAtlON  (S+l). 

AA  - COEFFICIENT  ARRAY  FOR  ONE  DIMENSIONAL  TFMPFRATURE  PROFILES. 
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NND  - NUMBER  OF  TEMPERATURE  NODES  ALONG  DIODF  AXIS 

NNH  - NUMBER  OF  TEMPERATURE  NODFS  ALONG  Y-AXTS  ( PFRPENDICULAR  TO 
DIODF  AXIS)  IN  SUBSTRATE.  MAXIMUM  VALUF  OF  NNHMX 

NNHMAX  - NNH  MAXIMUM  VALUE  WHICH  IS  THE  SAMF  AS  THE  *ROW  DIMENSION* 
FOR  ARRAYS  TEMPH.  TEMPS.  AND  TEMPS1. 

NNH8ND  - COLUMN  POINTER  FOR  LOCATING  THF  RESPECTIVE  ONE  DIMENSIONAL 
TEMPFRATURE  PROFILES  WITH  IN  THF  TWO  DIMENSIONAL  TEMPERATURE  ARRAY 
TEMPS1  FOR  SUBROUTINE  BANDA6  WHICH  SOLVFS  THF  INDIVIDUAL  ONF 
DIMENSIONAL  TEMPERATURE  PROFILES. 

S - ITERATION  counter 

DTHS  - NORMALIZED  CHANGF  IN  HEADER  TEMPERATURE  AT  NODE  (J.I) 

BETWFEN  THE  S AND  S+l  ITERATIONS. 

DTHSMX  - MAXIMUM  NORMALIZED  CHANGED  IN  HEADFR  TFMPERATURE 
BETWEEN  THE  S AND  S+l  ITERATIONS. 

ITPRDH  - HEADER  TEMPERATURF  ITERATION  PRINT  SINTINEL . WHEN 
ITPRDH=1  S AND  DTHSMX  ARE  PRINTfd  FOR  EACH  ITERATION 


100  CONTINUE 

SIMULATION  INITIALIZATION  CALCULATIONS 

IFLAG=0  IMPLIES  FIRST  TlMF  STEP  FOR  NEW  SIMULATION  AND  THAT 

simulation  initialization  calculations  should  bf  pfrformfd. 

IF < IFLAG.NE • 0 ) GO  TO  200 

set  tflag  to  indicate  that  the  simulation  has  bfen  initialized 

IFLAG=1 

INITIALIZE  TEMPFRATURE  ARRAYS 
DO  110  1=1 .NND 
DO  110  J=1.NNH 
temphi  j.  d=3oo.o 
TEMPS (J* I ) =300.0 
TEMPSl(J»I)=300.n 
110  CONTINUE 

EVALUATE  GRID  INCREMENTS 
DX=XL/ ( NND-1 ) 

THE  NUMBER  OF  HEADER  GRID  POINTS  ALONG  Y-AXIS  IS  <NNH«-1) 
DY=XDH/NNH 

EVALUATE  MODIFIED  GRID  COUNTS 


NNDM1=NND-1 
NNDM?=NND-2 
NNHM1 =NNH-1 
NNHM?=NNH-2 


200  CONTINUE 


pre-iteration  computations 

partial  coefficient  evaluation  which  will  accommodate  a varilable 

TIME  STEP  CONTROLLED  BY  THE  CALLING  PROGRAM 

A=-THKHX*DTIME*DY*DY 

B=-2 . 0*A+DY*DY*DX*OX*HDFN*HSPEC 

C=A 

Dl=DX*DX*DY*DYWHOEN*HSPFC 

D2=THKHY*DX*DX*DTIME 

INITIALIZE  ITERATION  COUNTER 
S=0 

INITIALIZE  TEMPS1  WITH  TEMPH  FROM  PREVIOUS  TTME  STEP 

DO  210  1=1 .NND 

00  210  J=1 . NNH 

TEMPS l (J.I) = TEMPH (J.I) 


SUBROUTINE  DHT2D2  (Continued) 


B20 


160 

161 

162 

163 

164 

165 

166 

1U 

m 

171 

172 

m 

175 

176 

m 

179 

180 
181 
182 

183 

184 

185 

11$ 

188 

189 

190 

191 


210  CONTINUE 


196 

111 

ffl 

201 

ill 
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in 

IW 
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C 

C- 

C- 

C 

C 

C 

C 


C 

C 

C 

C 

C 

c 

c« 

c- 

c 

c 

c 

c 


c 

c 

c 


c 

c 


300  CONTINUE 


ITERATION  ! OOP 

TRANSFEk  1TEMPS1  TO  TEMPS  50  THAT  NEW  TEMPS1  VALUE  MAY  BF  CALCULATED 

DO  305  1=1. NND 

DO  305  J=1#NNH 

TEMPS ( J*  I ) = TEMPS 1 (J.I) 

305  CONTINUE 

INCRFMENT  ITERATION  COUNTFR 
5=5+1 

EVALUATE  NNHMl  ONE  DIMENSIONAL  TFMPFRATURF  PROFILES  ALONG  X-AXIS 
WITH  COUPLING  WITH  ADJACENT  ROWS. 

400  CONTINUE 

IF  BlOCKING  BND  SPECIFIED  FOR  X=0  AND  X=XL»  FVALUATF  TEMPSl(l.I) 

IF ( IRND.NE . 1 ) GO  TO  500 

COEFFICIENT  EVALUATION  FOR  TEMPSltl.I)  WITH  RLOCKIN6  BND 
DO  410  1=2  * NNDM1 
11=1-1 
AA ( I T * 1 ) =A 
AA ( 1 1 * 2 ) =B 
“ l =C 


410 


C 

c- 

c- 

c 

c 

c 

c 

c 


8 

8 

8 


wwlii#  =D1*TEMPH (1*1) +D2* ( TFMPn ( I ) +TEMPS ( 2 . T > -2 . 0*TFMPS < 1 . 1 )) 

INCLUDE  BLOCKING  BND  FOR  TEMPS1 ( 1 • I ) 

AA ( 1 # 2 ) =AA ( 1 .2) +AA (1.1} 

A A ( NNDMP ,P)=AA( NNDMP ,2)+AA ( NNDM2 , 3 ) 

EVALUATE  TEMPS1 < 1 » I ) 

CALL  BANDA6 ( NNDM2 .2. 4*101. 101*4. A A . TEMP #1*PIVMIN) 

EQUATE  1EMP(1)  AND  TEMP<2>.  AND  TEMP ( NND ) AND  TFMP(NNDMl) 

FOR  RLOCKING  BND 
TEMP ( 1 ) =TEMP (2 ) 

TEMP ( NND ) =TEMP ( NNDM1 ) 

TRANSFER  HEADER  ROW  TEMPERATURE  TO  THE  HEADFR  TFMPFRATURF  ARRAY 
DO  420  1=1. NND 
TEMPS1 ( 1 » I ) =TEMP ( I ) 

420  CONTINUE 

500~CONTTNUE 

• * mm  • • • mm  ai  mm  • « w ■ • v mm  am  • • * ■.  • • * ••  «p  w w * • mm  a • ■ • * mm  w ^ a • mm  a mm  a « mm  • flp  w a w m mm  v 

IF  CONSTANT  TEMPERATURE  BND  SPECIFIED  FOR  X=0  AND  X*XL. 

EVALUATE  TEMPSltl.I) 

IF( IRND.NE. 0)  GO  TO  600 

COEFFICIENT  EVALUATION  FOR  TEMPSl(l.I)  AND  CONSTANT  TEMP  BND 
DO  510  I=2.NNDM1 
11=1-1 
AA ( I I . 1 ) =A 
AAt II .2)=E 
AA(fl»3)=( 

AA(I!,4)=C 
510  CONTINUE 


) 1 *TEMPH C 1 . 1 ) +D2* ( TFMPD { I ) +TEMPS ( 2 . 1 ) -2 . 0*TEMPS  1 1 . 1 ) ) 


INCLUDE  CONSTANT  TEMPERATURE  BND  FOR 
AA(1.4)=AA(1.4)-AA( 1. 1) *300.0 
A A ( NN0M2 » 4 ) = A A ( NNDM2 . 4 ) — A A ( NN0M2 . 3 ) * 


TEMPS1 t 1 . I > 
*300.0 


EVALUATE  TEMPSltl.I) 
CALL  BANDA6(NNr 


6 ( NNDM2. 2 .4.101. 101. 4. AA. TEMP. 1 »PI VMlN) 

TRANSFER  HEADER  ROW  TEMPERATURE  TO  THE  HFADFR  TFMPFRATURF  ARRAY 
DO  520  I=2.NNDM1 
TEMPSl (1*1) =TEMP ( I ) 
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520  CONTINUE 


600  CONTTNUF 


* < 
Hi 


g. 


g 


EVALUATE  TEMPERATURE  PROFILES  FOR  TFMPS1 (2*1)  - TFMPS1 ( NNHM? ,11 
WITH  BLOCKING  BND 
DO  650  J=2,NNHM? 

EVALUATE  COEFFICIENTS  FOR  TEMPS! {J, I) 

DO  610  I=2,NNDMl 
11=1-1 
AA(II,1)=A 
AA ! I I ,2) =B 
AA ( 1 1 * 3 ) =C 

AA 1 1 1 ,4 1 =D1  *TEMPH  1 J . 1 1 +02*  ( TEMPS!  1 J-l , 1 1 + TEMPS  ( J+1 « 1 1 
A-2.0*TEMPS!J,I) 1 
610  CONTINUE 

INCLUDE  BLOCKING  BND  FOR  TEMPSl <J, I) 

AA< 1,2)=AA(1,11+AA! 1,2) 

AA ( NNDM? , 2 ) =AA ( NNDM2 , 2 1 +A A < NNDM? , 3 ) 

EVALUATE  TEMP! J, I) 

CALL  BAN0A6 ( NN0M2 ,2.4, 101, 101, 4, AA, TEMP , 1 »PIVMlN) 

EQUATE  TEMP! I ) AND  TEMP! 21 » AND  TEMP ! NNDM 1)  AND  AND  TEMP 1NND ) 
FOR  BLOCKING  BNO 
TEMPI 1)=TEMP!2) 

TEMP  1 NND 1 STEMP ( NNDM1 1 

TRANSFER  header  row  temperature  to  THE  headfr  tfmperaturf  array 

DO  620  1=1, NND 
TEMPSl ( J , 1 1 =TEMp ( I ) 

620  CONTINUE 
630  CONTINUE 

too’contTnue""  ~~ 

EVALUATE  TEMPERATURE  PROFILES  FOR  TEMPSl !?, I)  - TEMPSl ( NNHM2, I ) 
WITH  CONSTANT  TEMPERATURE  BND 
DO  730  J=2 , NNHM2 

EVALUATE  COEFFICIENTS  FOR  TEMPSl (J,I) 

DO  710  I=? » NNDM1 
11=1-1 
AA ! I I , 1 ) =A 
AA ! I I , 2 ) =B 
AA( IT,3)=C 

AA ! 1 1 , 4 1 =0 1*TEMPH ! J, I ) +02* ! TEMPS! ! J- 1 » 1 1 + TEMPS ! J+ 1 . 1 1 
«-2.0*TEMPS! J,I) 1 
710  CONTINUE 


INCLUDE  CONSTANT  TEMPERATURE  BND  FOR  TEMPSl! J, I) 
AA! 1 ,4)=AA ! 1,4) -AA! 1,1 ) *300.0 
AA ! NNDM2 , 4 ) = A A ! NNDM? , 4 1 — AA ! NNDM? ,31*300.0 


EVALUATE  TEMPIJ.I) 

CALL  BAN0A6(NNDM2,2»4, 101,101 ,4. AA , TEMP, 1 ,PI VMIN) 

TRANSFER  HEADER  ROW  TEMPERATURE  TO  THE  HEADFR  TEMPFRATURE  ARRAY 
DO  7?0  I =2, NNDM 1 
TEMPSl ( J , 1 1 =TEMP ( I ) 

720  CONTINUE 
730  CONTINUE 


800  CONTINUE 

■ ■»*••*****■•*  w mt  mt  w • mm  mm  • rnmmt  mm  tm  mm*m  w tm  mm  » mm  •rnmmmm  ^mm  mmmm^mmmm 

IF  BIOCKING  BNO  SPECIFIFD  EVALUATE  TEMPSl 1NNHM1 , 1 ) 
IFURND.NE.l)  GO  TO  900 

COEFFICIENT  EVALUATION  FOR  TpMPSI 1NNHM1 , 1 1 WITH  BLOCKING  BND 
DO  810  I=?,NNDM1 


11=1-1 


SUBROUTINE  DHT2D2  (Continued) 
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AA(II#1)=A 

AA(II#2)=B 

AA(II#3)=C 

AA  ( 1 1 » 4 ) =D 1 *TEMPH  ( NNHMl , 1 ) +D2*  ( TFMPS1  ( NNHM2  , T ) +300 . D 
A-2.0*TEMPS(NNHMl#I>) 
ftlO  CONTINUE 

INCLUDE  BLOCKING  BND  FOR  TEMPS1 ( NNHMl * I ) 
AA(1»2)=AA(1'2)+AA(1'|) 

A A ( NNDM? # 2 ) =AA ( NNDM2 , 2 ) +AA ( NNDM2 , 3 ) 

EVALUATE  TEMPSl (NNHMl# I) 

CALL  BAnDA6(NNuM2»2»4» 101*101 #4, AA'TEMP, 1 #PlVMlN) 

EQUATE  TEMP ( 1 ) AND  TEMP(2>#  AND  TEMP ( NNDM 1)  AND 
TEMP ( NND ) FOR  BLOCKING  RND 


TEMP ( 1 ) = TEMP ( 2 ) 

TEMP (NNn)=TEMP( NNDM  1 ) 


TRANSFER  HEADER  ROW  TEMPERATURE  TO  THE  HEADER  TFMPERATURF  ARRAY 
DO  820  1=1 *NND 
TEMPSl (NNHM1 . I )=TEMP( I ) 

820  CONTINUE 


900  CONTINUE 


IF  CONSTANT  TEMPERATURE  BND  SPECIFIED#  EVALUATE 
TEMPSl ( NNHMl • I ) 

IF(IBND.NE.O)  GO  TO  100n 

COEFFICIENT  EVALUATION  FOR  TFMPS1 (NNHMl . I ) WTTH  CONSTANT 
TEMPFRATURE  BND 
DO  910  I =2' NNDM 1 
11=1-1 
AA( IT# 1)=A 
AA ( I T # 2) =8 
AA ( I I ' 3) =C 

AA ( 1 1 # 4 ) =D1*TEMPH ( NNHMl . I ) +D2* ( TFMPS1 < NNHM2 . T ) +300 . 0 
4-2.0*TEMPS(NNHM1 . I) ) 

910  CONTINUE 

INCLUDE  CONSTANT  TEMPERATURE  BNO  FOR  TEMPSl (NNHMl, I) 

AA ( 1 » 4 ) =AA ( 1 » 4 ) -300 . 0*AA (I'l) 

AA ( NNDM2  * 4 ) = A A ( NNDM? , 4 ) -300 . 0*A A ( NN0M2 # 3 ) 

EVALUATE  TEMPSl (NNHMl# I) 

CALL  BAnDAMNNDM2'2#4#  101,101  #4# AA# TEMP.  l.PIVMlN) 

TRANSFER  HEADER  ROW  TEMPERATURE  TO  THE  HEADFR  tfmpfraturf  ARRAY 
DO  920  I =2 » NNDM 1 
TEMPSl ( NNHMl . I ) =TEMP ( I ) 

920  CONTINUE 


:TEMP(I> 


1000  CONTINUE 


convergence  CHECK 

EVALUATE  MAXIMUM  NORMALTZFO  CHANGE  IN  TEMPERATURE  FOR  THIS 


IterWionMAXIMUM  normaltzfo  change  IN  TEMI 

DTHSMXSO.O 

no  mio  1=2 # nndmi 

d9hI= ( TFMpli? J?I ) -TFMPS ( J' I M /TEMPSl ( J# I) 
IF(DTHS.GT.DTHSMX)  dthsmx=dths 


1010  CONTINUE 


PRINT  ITERATION  SUMMARY 
IF ( ITPRH.EQ, 1 ) WRITE (6,1 0?0 ) S.OTHSMX 
1020  FORMAT (3H  S=, 12, SX, 7HDTHSMX= ,E8. 3) 


IF  DTHSMX  .LE.  DTHMAX  THE  SPEC I 
IF (DTHSMX.LE.DTHMAX)  GO  TO  2000 


IFIED  C0NVER6FNCE  HAS  BEEN  ACHIEVED 


CONVERGENCE  HAS  NOT  BEEN  ACHIEVED.  IF  THE  ITERATION  COUNT 


T 


B2  3 


SUBROUTINE  DHT2D2  (Continued) 


IS  LFSS  than  thf  maximum  sfecififd  start  another  itfration 

IF(S.LT.SMAX)  GO  TO  300 

THE  MAXIMUM  NUMBFR  OF  ITERATION  HAS  BEEN  PERFORMED  WITHOUT 
ACHIFVING  CONVERGENCE.  WRITE  ERROR  MESSAGE  AND  TERMINATE  XGT 
WRITF(6» 1030)  S.DTHSMX 

1030  FORMAT (///»46H  *****  CONVFRGFNCF  FAILURE  IN  SUBROUTINE  DHT2D. 

A/.27H  *****  EXECUTION  TERMINATED./. 

A9H  *****  S=»I2»5X.7HDTHSMX=.EB.3»/# 

A34H  *****  INCREASE  SMAX  AND/OR  DTHMAX ) 

STOP  ************************************************************** 

20 oo ’continue'' 

TRANSFER  NEW  SUBSTRATE  TEMPERATURE  VALUFS  TO  ARRAY  TEMPH 
00  2020  irl.NNO 
DO  2020  Jzl.NNH 
TEMPH ( J , I ) sTEMPS t < J . I) 

2020  CONTINUE 

RETURN  TO  CALLING  PROGRAM  WITH  NEW  HEADER  TFMPERATURES 
RETURN 
END 


C 

C 

c 


c 

c 

8. 

c 

8 


C 

c 


SUBROUTINE  DOPLG 


B24 


*TEMP(1) . DOPLG 


i 

5 

6 

l 

10 

tt 

13 

14 

15 

I! 

19 

IS 

22 

23 

24 

25 

26 

27 

28 
29 

I? 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 


11 


SUBROUT I NE  DOPLG ( NND , DX , XME T , XDFPEP » DOPL . DOPU » DOPEP , DOPL OG » I DB4 ) 
SUB  DOPLG  GENERATES  A LOGARITHMIC  DOPING  PROFILF  IN  ARRAY  DOPLOG 


NND 

DX 

XMET 

XDEPFP 

88pep 

DOPU 

DOPLOG 

IDB4 


- NUMBER  OF  NODES  ALONG  DIODE  AXIS 

- NODE  SPACING  ALONG  DIODE  AXIS 

- METALLURGICAL  JUNCTION  LOCATION 

- EPITAXIAL  LAYER  BOUNDARY 

- L06FR  DOPING  CONCENTRATION 

- EPITAXIAL  LAYER  DOPING  CONCENTRATION 

- UPPER  DOPING  CONCENTRATION 

- LOGARITHMIC  DOPING  PROFILE 

- DEBUG  PRINT  SENTINEL 


DIMENSION  DOPLOG (NND) 

DOPLGL=ALOGl 0 ( DOPL ) 

DOPLGESALOG 1 0 ( DOPEP ) 

DOPLGU=ALOGl 0 ( DOPU ) 

GENERATE  IMPURITY  PROFILE 
DO  10  1:1 » NND 
X:<I-1)*DX 

IF ( X .LE.  XMET ) DOPLOG ( I ) rDOPLGL 

IF ( ( X .GT.  XMET). AND. (X  .LE.  XDFPEP) ) DOPLOG ( I ):DOPLGF 
IF ( X .GT.  XDEPEP) DOPLOG ( I ):DOPLGU 

IF ( IDB4  .NE.  0)WRITE(6»11)I.DOPLOG(I) 

FORMAT ( IX. 7HDOPLOG( . I3» ?H):*F6. 3) 


10  CONTINUE 

?65URN 


SUBROUTINE  EFIELD 


B25 


*TEMP(1  ) 
1 
2 

3 

4 

5 

6 
7 

a 

9 

10 


EFIELO 


& 


C 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


s 

c 

c 

c 

c 


c 

c 


c 

c 


SUBROUT I NE  EF IELD ( TEMP , XL , NND . I TFR » R AC . XLDEP . XUDEP . CUR . DOPL . DOPU . D 
AOPEP . XMET . XOEPEP . NP  1 , VLRR . VRBR , VTOT . ECRIN . EmAXB . I0R4 . EIN t T ) 

SUB  ffield  evaluates  thf  diode  ei  ectric  field  from  contact  TO 
contact  assuming  that  the  diode  current  consist  totally 
OF  A DRIFT  CURRFNT.  THF  RESULTING  F FIELD  IS  TfMPFRaTURE 
and  doping  concentration  depfndfnt.  thermal  fouilirrium 
majority  and  minority  carrier  concentrations  arf  assumed. 


variable 

ECRINC I ) - 
fi)  - 


TEMP 

iter 

&> 

RAC 

XLDEP 

XUDEP 

DOPL 

DOPU 

CUR 

VLBR 

VRBR 

VTOT 

XN 

P 

XNI 

XMET 

XDEPFP 

DOPEP 

NP1 

IDB4 

EINIT 

0 

XK 

F 

PF 


DEFINITIONS 

ELECTRIC  FIELD  INTENSITY  AT  NODF 
TEMPERATURE  AT  NODE  I 
NO.  OF  ITERATIONS  USEO 


IN  CALCULATING  F FIFLD 


“I? 


-TOTAL  LENGTH  OF  DlOOl 

’N0*-d^sirKSa8curacy  in  Computation  of  fcrincd 

-LEFT  DEPLETION  REGION  BOUNDARY 
-RIGHT  DEPLETION  REGION  BOUNDARY 
-DOPING  CONCENTRATION  OF  LEFT  BULK  RFGION 
-DOPING  CONCENTRATION  OF  RIGHT  BULK  REGION 
-DIODE  CURRENT  DENSITY 
-VOLTAGE  ACROSS  LEFT  BULK  REGION 
-VOLTAGE  ACROSS  RIGHT  BULK  REGION 
-TOTAL  VOLTAGE  ACROSS  DIODE 
•CONCENTRATION  OF  ELECTRONS 
-CONCENTRATION  OF  HOLES 
•INTRINSIC  CARRIFR  CONCENTRATION 
-METALLURGICAL  JUNCTION 
-EPITAXIAL  boundary 
•epitaxial  impurity  concentration 

*1  FOR  NP  CONFIG. »0  FOR  PN  CONFIG. 

calculating  flectric  field 

-UNIT  CHARGE  ON  AN  ELECTRON. IN  COULOMBS 
-BOLTZMANN  CONSTANT, ELECTRON  VOLTS/DFGRFE  KFLVIN 
•TOTAL  CURRENT  DENSITY  EOUATION 
-DERIVATIVE  OF  F WITH  RFSPECT  TO 


ELFCTRIC 


FIFLD 


DIMENSION  ECRINCNND) .TEMP(NND) 

DATA  G/1.6E-19/XK/8.62E-5/ 

CALCULATE  DEPLETION  REGION  BOUNDARIES  TO  NEAREST  NODE 
DX=XL/(NNO-l) 


NSXLDEP/DX  +1.! 
MSXUDEP/OX  +1.! 


C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

? 

c 


50 


10 


DEBUG  OUTPUT 

IF ( IDB4  .EG.  1)WRITE(6*S0)DX • N.M.NND 

FORMAT(1X,3HDX=,F10.5» 1X,?HN=»I3. IX. 2HMr , 13. 1 X. 4HNNDS, 13) 


DO  LOOP  FOR  CALCULATING  ELECTRIC 
DO  20  Lsl.NND 


FIFLD  AT  EACH  NODE  L 


11 


CALCULATE  XNV.XNC. EFFECTIVE  DENSITY  OF  STATF 
XNVSl . 02E19* ( TEMP(L) /300.0) **1 .5 
XNC=2. 8E19*( TEMP (LT/300.0) **1.5 

CALCULATE  INTRINSIC  CARRIER  CONCENTRATION  XNI 
XNI=(XNV**.5)*<XNC**.5)*EXP(-<1.1/<2*XK*TEMP(L>  > ) ) 

CALCULATE  HOLE  AND  ELECTRON  CONCENTRATIONS 
LL=L-1 

IF  N-P  DIODE  CONFIG. (NPtsI )GO  TO  18 
IF ( NP1  .EG.  1 )GO  TO  18 


FOP  P-N  DIODE  CONFIG. 

IF (LL*OX  .LE,  XMFTJgO  TO  il 

IF(<LL*DX  ,GT.  XMET).AND.ILL*DX  .LE.  XOFPEP) )GO  TO  12 
IF(LL*DX  .GT.  XDEPEP1GO  TO  13 

FOR  LOWER  (P)^SIDE 

H0L=.5*D0PL+.5*I (TDOPL**2)^4*CXNT**P) )**.5) 
ELE=(XNI**2)/HOL 


r 


icu 


133 

146 


1 


47 


48 


>0 
151 

i 1! 
ip 
It* 
ill 


C 

C 


c 

c 


c 

c 

l 


C 

C 


C 

C 


c 

c 


c 

c 

c 


c 

c 


c 

c 


1 
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SUBROUTINE  EFIELD  (Continued) 


60  TO  16 


FOR  FPITAXIAL  LAYFR ( N-TYPF ) 
12  ELEz.5*D0PE  “ 


P+.5*<  < CD0PEP**2)+4*(XNI**2)  )**.5) 
H0Lz(XNI**2)/ELE 
GO  TO  16 


FOR  UPPER  <N)  SIDE 


13  ELE=.5*D0PU+.5*< ( <D0PU**2)+4*(XNT**2) )**.5) 
HOLz { ; 


;XNI**2)/ELF 
GO  TO  16 


18 


FOR  N-P  DIODE  CONFIG. 


IF ( LL*DX  . LE  , XMFTJGO  TO  14 

IF((IL*DX  -GT.  XMET) . AND. (LL*DX  .LE.  XDFPEPHGO  TO  15 
IF(LL*DX  .GT.  XDFPEP)GO  TO  21 


FOR  LOWER  <N)  SIDE 

14  ELE=.5*D0PL+.5*( ( <D0PL**2)+4#(XNT**2) )**.5) 
HOL= (XNI**2) /ELF 
GO  TO  16 


FOR  FPITAXIAL  LAYER (P-TYPE ) 

15  H0L=.5*D0PEP+.5*(  7 (00PEP)**2)+4*(XNI**2)  >**.51 
ELEs<XNI**2)/H0L 
GO  TO  16 


FOR  UPPER  (P)  SIDE 
.5*<(T 


21  H0LS.5*D0PU+ 

ELEz ( XNI **2 ) /HOL 


D0PU**2)+4*(XNT**2> )**.5) 


INITIALIZE  ELECTRIC 
16  EZEINIT 


FIELD 


DETERMINE  DOPING  FOR  POSITION  ALONG  DIODE  AXTS 
IF(LL*DX  .LT.  XMETJDOPlzDOPL 

!F((LL*DX  .GT.  XMET) .AND. (LL*DX  .LE.  XDEPEP) IDOPISDOPEP 
IF ( LL*DX  .GT.  XDEPEP) DOPlzDOPU 


DEBUG  OUTPUT 

IF(IDB4  .EQ.  1 ) WRITE (6*65) ELE * HOI  »TEMP(L) 

95  FORMAT(/»1X»4HELFS,E10.6,1X*4HHOLS,F10.5.1X*5HTEMPS.F6.1 ) 


C 

C 


CHECK  FOR  CONVERGENCE  OF  NEWTON-RAPHSON  METHOD 
CALL  EMOBS ( DOPl . FMAXB . TFMP ( L ) » EMO . EMOE * FMOT ) 
CALL  HM0BS(D0P1.EMAXB*TfMP(L) »HMO,HMOE»HMOT) 
curprzq*Tele*emo+hol*hmo ) *EMAXB 

IF(CURPR.GE.CUR)GO  TO  40 


c 

c 


SET  ELECTRIC 
ESEMAXB 
GO  TO  19 


FIELD  TO  DEFAULT  VALUE 


C 

C 


DO  LOOP  FOR  CALCULATING  ELECTRIC 
40  DO  17  izl.ITER 


FIELD 


C 

C 

c 


CALL  SUBROUTINES  FOR  FINDING  MOBILITIES 
CALL  EM0BS(D0P1.E*TFMP(L) .EMO. EMOE. EMOT) 
CALL  HMOBSIDOPl.F.TFMP(L) .HMO.HMOE.HMOT) 


C 

C 

C 


TE  F(E>. TOTAL  CURRENT  DENSITY  EQUATION 
F*EM0*E ) ♦ ( Q*HOL*HMO*E ) 


•Z(Q* 

CALCULATE  PF(E) .DFRIVATIVF  OF  TOTAL  CURRENT  DENSITY 
° TO  FLI 


C 

C 


WITH  RESPECT  To  ELECTRIC  FIELD 

PFs ( Q*ELE*EMO ) ♦ ( Q*h5l*HM0 ) ♦ ( Q*ELF*E*FMOF ) ♦ ( Q*HOL*F*HMOE ) 


W 


BUG  OUTPUT 

(IDB4  .EQ.  1)WRITE(6»45)E*F*PF.DE»EMO*HMO.FMOE.HMOE 


45  FORMAT ( 1 X.2HES.F9.4. 1X.2HFz*F9»4* lX»3HPFz.E9,4. 1X.3HDEb.E9. 4. IX. 4H 
AEMOz  »E9. 4 » IX • 4HHMOZ.F9.4 . 1 X . RHEMOEz , F9 . 4 » IX. 5HHM0FZ.E9.4) 


SUBROUTINE  EFIELD  (Continued) 


CALCULATE  NEW  VALUE  OF  F 
E=E+OE 

„ CHECK  FOR  DESIRED  ACCURACY 
17  IF(ARS(DE/E) .LE.  RAC ) 60  TO  19 

DEBUG  OUTPUT 

19  IF ( IDB4  .EQ.  1 > WRITE (6»47 )E 
47  FORMAT ( 1X»7HEFIFLDS»E10.5»/) 

SET  ELECTRIC  FIELD<ECRIN(L>>  EQUAL  TO  FINAL  VALUE  OF  E 

20  ECRIN(L)=E 

CALCULATE  VOLTAGF  ACROSS  LEFT  BULK  REGION 
VLBRrECRIN< 1>*.5*0X 
DO  25  J=2.N 

25  VLBRrVLBR+ECR IN ( J ) ♦DX 

CALCULATE  VOLTAGF  ACROSS  RIGHT  BULK  REGION 
VRBRrECR I N ( NND ) * „ 5*DX 
NL=NND-1 
DO  30  J=M,NL 

30  VRBRsVRBR+ECRlN( J)*DX 


CALCULATE  DEPLETION  REGION  VOLTAGE 

NXLSN+1 

NUL=M-1 

VDEPrO 


DO  35  LL=NXL»NUL 
35  VDEP=VDEP+ECRIN(LL)*DX 

CALCULATE  TOTAL  VOLTAGE 
VTOTrVLBR+VRBR+VDEP 


RETURN 

END 


I 


SUBROUTINE  EMOBS 


B28 


♦TEMP(l) .EMOBS 
1 

! 

4 


S 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 


c 

c 


c 

c 


c 

c 


c 

c 


c 

c 


SURROUTINF  EMOBS ( OOP 1 »Et  »TEMP»EMO»EMOE»FMOT) 


EMOBS  EVALUATES  ELECTRON  MOBILITY » DERIVATIVE  OF 
WITH  RESPSCT  TO  THE  ELECTRIC  FIELD  AND  THE  DFRIV 
MOBILITY  WITH  RFSPECT  TO  temperature. 


electron  mobility 
dfrivative  of  electron 


DOP1 


w 


TEMP 

EMO 

EMOE 


- DOPING  CONCENTRATION 

- ELECTRIC  FIELD 

- TEMPERATURE 


- D^iyATIVE°8JLFLECTRON  MOBILITY 


EMOT 

IFLAG 


__  WITH  RESPECT  TO 

- deriva^veTSfDflectron  mobility  with  respect  to 
temperature 

- sentinel  for  initialization  calculations 


ASSIGN  CONSTANTS  FOR  ELECTRON  MOBILITY  CALCULATIONS 
DATA  EMOO.EN.ES.FA»EF.EB.ALFA#IFl AG/1400 . 0# 3. 0E1 6. 350 .0 . 3.5F3» 
Afl.8f7.4E3*2.5»0/ 


INITIALIZATION  calculations 
IF(IELAG.NE.O)GO  TO  10 


TB=  FF*EA*EA 
TC=EB**(-2) 

TE=  1 .0/<EMOO*3O0.0**ALFA) 


SET  SENTINEL  TO  SKIP 
IFLAGrl 


INITIALIZATION  CALCULATIONS 


10  CONTINUE 


E=ABS(E 

DOP=AB! 


1) 

>(D0P1 ) 


EVALUATE  REMAINING  PARAMETERS 
TA=DOP/ ( DOP/ES+EN ) 

TD=EA*E+TB 


evaluate  intermediate  function 

H=1 . 0+TA+ ( 1 . 0/TD+TC > *E*F 


DERIVATIVE  OF  H WITH  RESPECT  TO  f 
HE= ( ( 2 . 0*TD-EA*E ) ✓ ( TD*TD ) *2 . 0*TC ) *E 


G=TE*TEMP**ALFA 
DERIVATIVE  OF  G 
GT=ALFA*G/TEMP 


WITH  RESPECT  TO  T 


EVALUATE  ELECTRON  MOBILITY 
EMOsi . 0/ ( G*SQRT (HI ) 


evaluate  derivative  OF  ELECTRON  MOBILITY  WITH  RESPECT  TO  E 
EMOEs-HE/ ( 2 . 0*G*H** 1.5) 


evaluate  derivative  OF  FLECTRON  MOBILITY  WITH  RFSPFCT  TO  T 
EMOT=-GT/ ( SORT ( H ) *G*G ) 


RETURN 

END 


vDO>>i<M»toiiN»-o<oo9'*JO‘y>«oiiN>t-o»ooB«-j(ruuco«v)^osoo>^ia>cH"Con\)k-o<oOD'jo'uw?oMw*-o\Ocx>^io'uwro«\>^-o\ooB^O'Ui-co»w*-o\Doc^jO'm.ci>iv)k-  s 
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SUBROUTINE  EPROF 


P<1> .EPROF 

C 

C 

8 

C 

c 

c 

c 

c 

c 

c 

c 

8 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


8 

c 

c 

c 

c 

c 

c 


SUBROUT I NE  EPROF ( NND . DX . XDEPL » XMFT . XDEPEP . XOFPU . XL . SOOPL . 

H SDOPFP . SDOPU . EMA  X . E * IDB3  # FLAGEB l 

SUB  EPROF  ANALYTICALLY  GENERATES  THE  DEPLETION  REGION  ELFCTRIC 

EKttLBWflf  3KB.8H. ’flWxT r^Tklta^rokeiFo. 
^ftSSlRfSi^El^Sb.rTS'V^^c^Sg^S^glgVSI'gguSg^?'? 

VARIABLE  DEFINITIONS 


NND 
DX 

XDEPL 

XHET 

XDEPU 

IkoPL 

SDOPEP 

SDOPU 

EMAX 

E 

IDB3 

FLAGEB 


- NUMBER  OF  NODFS  ALONG  DIODE  AXIS 

- NODE  SPACING  ALONG  DIOOF  AXIS 

- lower  depletion  REGION  boundary 

- METALLURGICAL  JUNCTION  I OCATION 

- UPPER  DEPLETION  REGION  BOUNDARY 

: ESSgB  hEPSYio N region  SPACE  CHARGE 

- EPITAXIAL  DEPLFTION  REGION  SPACE  CH 

- UPPER  OFPLFTION  REGION  SPACE  CHARGF 
““C  FIELD  AT  METALLURGICAL 


- ELECTRIC  _____ 

- ELECTRIC  FIELD  ARRAY 

- SUB  EPROF  DEBUG  PR 

- DEPLETED  BULK  REGI 


CHARGF 

GF 

JUNCTION 


PRINT  SFNTINEL 
‘ ON  SFNTINEL 


NAMEL 1ST  /DEBUG/NNQ » OX  * XDEPL . XMFT  * XOEPEP . XOFPU  # XL » SOOPL . 
ASDOPFP » SDOPU # EMAX »Q » PERM . FEP 

DIMENSION  E(NND) 

INTEGER  FLAGEB 

Bata  PERM/f7o6E-i2/ 

EVALUATE  LOWER  DEPLFTION  REGION  FDGE 

xdepl=xmet-peRm/q*emax/sdopl 

EVALUATE  UPPER  DEPLETION  REGION  WIDTH 

xdepu=xmet+perm/q*emax/sdopep 

DOES  UPPER  DFP  REG  EDGE  FALL  WITHIN  EPIT  REG? 

IF ( XDEPU. LE. XOEPEP)  60  TO  100 


C 

C 

C 

C 

C 


C 

C 

c 

c 


c 

c 


UPPER  DEP  REG  EDGE  OCCURS  WITHIN  BACKGROUND  DOPING 

iEv^WK4^iJSiio5rfkS,AxSE^^Tx5?¥( 

EVALUATE  UPPER  OEP  REG  EDGE  WITHIN  BACKGROUND  OOPING 
XDEPU=XOEPEP+PERM/Q*EEP/SDOPU 

100  CONTINUE 

NEW  XDEPL  AND  XOFPU  HAVE  BEEN  EVALUATED 
ARE  DEP  REG  EDGFS  VALID? 

RESET  FLAGEB 
FLAGFB=0 

IFIXDEPL.GE. 0.0. AND. XDEPU. LE. XL)  GO  TO  200 

A BULK  REG  HAS  BEEN  DEPLETED.  SFT  FLAGEBrl  AND  SET 
THE  ERRONEOUS  DEP  REG  BNO  EQUAL  TO  THE  RESPFCTIVF  BULK 
BND.  IF  THIS  CONDITION  YIELDS  Al  FAIN  .LT.  AFRMAX 
THE  DIODE  OESIGN  IS  BAD  AND  XQT  WILL  BE  TERMINATED 
IN  Si  IB  BKDEPL 
FLAGEBrl 

IFIXDEPL.LT. 0.0)  XI 
IF ( XDEPU. GT . XL)  XI 

200  CONTINUE 


(DEPLrfl.O 

(DEPUrXL 


GENERATE  ELECTRIC  FIELD  PROFILES 
DO  200  Nr 1, NND 


DETERMINE  NODE  POSITION 
Xr(N-l)*DX 

EVALUATE  E FIELD  AT  NODF 
IFfX.GE. XDEPL)  GO  TO  210 


I OR  X 


B30 


SUBROUTINE  HMOBS 


B31 


♦TEMP ( 1 ) 
1 
2 

3 

4 

5 

6 

7 

8 
9 


•HMOBS 


SUBROUTINE  HMOBS (DOPl . El .TEMP.HMO.HMOE.HMOT) 


HMOBS  EVALUATES 
WITH  RESPSCT  TO 
MOBILITY  WITH 


HOLE  MOBILITY#  DERIVATIVE  OF  HOLE  MOBILITY 
_ THE  ELECTRIC  FIELD  AND  ThF  DERIVATIVE  OF  HOLE 

RESPECT  TO  TEMPERATURE. 


DOP1 

El 

TEMP 

HMO 

HMOE 

HMOT 


- DOPING  CONCENTRATION 

- ELECTRIC  FIELD 

- TEMPERATURE 

- HOLE  MOBILITY 

- DERIVATIVE  OF  HOLE  MC 


VE  OF  HOLE  MOBILITY  WITH  RESPECT  TO 
FIELD 

VE  OF  HOLE  MOBILITY  WITH  RESPECT  TO 


IFLAG  - sIntInElT^R  INITIALIZATION  CALCULATIONS 

ASSIGN  CONSTANTS  FOR  HOLE  MOBILITY  CALCULATIONS 
DATA  HMOO# HN#HS#HA»HF # HR# ALFA# IFI  AG/480 .0 #4. OElB#flt .0»6. 1E3# 
A1.6»?.5E4#2.3»0/ 

INITIALIZATION  CALCULATIONS 
IF ( IFLAG.NE. 0)GO  TO  10 
TB=  HF#HA*HA 
TC=HB**<-2) 

TE=  1 .0/(HMOOA300.0**ALFA) 

set  sentinel  to  skip  initialization  CALCULATIONS 
IFLAG=1 

10  CONTINUE 

E=ABS(E1 ) 

DOP=ABS(DOPl) 

evaluate  remaining  parameters 

TArDOP/ ( DOP/HS+HN » 

TD=HA*E+TB 

EVALUATE  INTERMEDIATE  FUNCTION 
H=1 . 0+TA+ < 1 . 0/TD+TC > ♦E^F 

DERIVATIVE  OF  H WITH  RESPECT  TO  F 
HE= ( ( 2 . 0#TD-HA*F ) / ( TO*TO > +2 . 0*TC ) *E 


G=TE^TEMP**ALFA 
DERIVATIVE  OF  G 
GT=ALFA*G/TEMP 


WITH  RESPECT  TO  T 


EVALUATE  HOLE  MORILITY 
HMO=1.0/(G*SQRT(H) ) 

EVALUATE  DERIVATIVE  OF  HOLE  MOBILITY  WITH  RFSPECT  TO  E 
HMOE=-HE/ ( 2 . 0*G*H** 1.5) 

EVALUATE  DERIVATIVE  OF  HOLE  MOBILITY  WITH  RFSPECT  TO  T 
HMOT=-GT/ ( SORT ( H ) *G^G  > 


SUBROUTINE  INALFA 
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♦TEMPI 1> 
1 
2 
3 


INALFA 


SUBROUT I NE  I NALF A ( NND » XnEPL . XMET . XDEPEP . XDEPU . DX , DOPL » DOPEP , 
ADOPU, ALFATD. E.T.ALFAIN.TDB1.NP1) 

SUB  INALFA  EVALUATES  THF  INTEGRAL  OF  ALFA  OR  THE  AVALANCHE 
BREAKDOWN  INTEGRAL  FROM  THE  FIRST  GRID  POINT  PASS  XDEPL  TO 
THE  LAST  GRIO  POINT  BEFORE  XDEPU. 


NND 

XDEPFP 

XDEPU 

DX 

DOPL 

DOPEP 


- NUMBER  OF  NODFS  ALONG  DIODE  AXIS 

- ^os|ti8nPofem|talEurgical  junction 

- epitaxial  layer  boundary 

- upper  oep  region  bnd 

- node  spacing  along  diodf  axis 

- lower  doping  concentration 


DOPL  - LOWER  DOPING  CONCENTRATION 
DOPEP  - EPITAXIAL  DOPING  CONCENTRATION 
OOPU  - UPPER  DOPING  CONCENTRATION 

ALFATD  - AVALANCHE  IONIZATION  COEFFICIENT  TFMPERATURF 


_ DEPENDENCE  PARAMETER  „ 

E - ELECTRIC  FIELD  ARRAY 

T - TEMPERATURE  ARRAY 

ALFA  IN  - SOCIAL  INTEGRAL  OF  ALFA 

IDB1  - i FBUG  PRINT  SENTINEL 

NP1  - JUNCTION  ORIENTATION*  1 - NP»  0 - PN 

DEBUG  NAMELIST 

NAMEL I ST  /DEBUG/NNO  * XDEPL . XMET . XnEPEP  * XDEPU . DX . DOPL  * DOPFP  # DOPU . 
AALFATD » ALFA I N * NDFPL » NDEPU 

DIMENSION  E(NND) *T(NND) 

DETERMINE  THE  INSIDE  DEPLETION  RFGION  GRIO  POINT  LIMITS 
NDEPL=XDEPL/DX+? 

NDFPUSXDEPU/DX+1 

INITIALIZE  GRID  pointer 
n=ndfpl-i 

INITIALIZE  the  ALFA  INTF6RAL 
ALFATNsO.O 

100  CONTINUE*' 

INTEGRATE  ALFA  FROM  N=NOEPL  TO  Nr NDEPU 
N=N+1 

X=(N-1»*DX 

IF ( X.GE . XMET)  GO  TO  14 

EVALUATE  ALFA  IN  LOWER  DEPLETION  REGION 
CALL  IONCOF  ( T ( N)  *E(N)  .DOPL.NPl  ,AI. FA . ALFATD) 

GO  TO  19 


14  CONTINUE 

IF (X.GT. XDEPEP)  GO  TO  IS 

EVALUATE  ALFA  IN  THE  EPITAXIAL  LAYER 
CALL  IONCOF(T(N) *E(N) .DOPFP.NPI .ALFA 
GO  TO  19 

15  CONTINUE 


ALFA. ALFATD) 


EVALUATE  ALFA  IN  UPPER  DEPLETION  REGION 
CALL  IONCOF(T(N) .E<nT.OOPU.NP1.AI  FA. ALFATD) 


19  CONTINUE 


SUBSTRACT  HALF  OF  FIRST  ALFA 
I F ( N . EG . NDEPL ) ALFA I N=ALF A I N-ALF A /2 . 0 

SUBSTRACT  HALF  OF  LAST  ALFA 

IF (N.EQ. NDEPU)  ALFA I NsALFA IN-ALFA/2.0 


ALFA IN-ALFA/2.0 


HAVE  ALL  ALFAS  BEEN  SUMMED? 
IF(N.LT. NDEPU)  GO  TO  100 


r 


24 

25 

26 

27 

28 

29 

30 

31 

32 

a 

35 

» 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

Si 

50 

51 

52 

53 

54 

55 

56 
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SUBROUTINE  IONCOF 


•TEMP ( 1 ) 
1 

l 

4 

5 

6 
7 
6 
9 

10 
11 
1 
1 

14 

15 

1$ 

18 

19 

20 
21 
2 


IONCOF 

C 
C 
C 
C 
C 
C 

r* 

G 
r 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


SUBROUT I NE  I ONCOF  L T , EE . OOP . NP1 , At F A ♦ ALFATD ) 

IONCOF  - EVALUATES  AVAL ANCHE  IONIZATION  COEFFICENT  AS  A FUNCTION 
OF  jFMpfKATURE  AND  ELECTRIC  FIELD.  THE  TEMPeRATURF  DEPFNDENCE  IS 

ALFATD.  THF  avalanche  COEFFICIENT  evaluated 
CORRESPONDS  TO  The  minority  CARRIER  conc  FOR  EPITAXIAL  layer 
MATERIAL* 


T - TEMPERATURE.  DFG  K 

it  - electric  field.  VOLTS/CM 

DOL  - OOP I NG  CONCENTRATION 


ALFA  - AVALANCHE  I0NI7k  ION  COErFICZc  NT 

Tfr,.  COEFFICIENT  TemPERATURF  DFPENDFNCF 


ALF.'TD  - AVALANCHE  I0NI7A 

PARAMETER.  1/DFG  K 
NP1  - DIODE  ORIENTATION. 


I - NP.  0 - PN 


C 

C 


C 

C 

C 

C 

C 

c 


c 

c 


initialize  parameters  r<*  ionization  coefficient  formula  .'ton 

DATA  EA.EB/3.8E6.1.7E6/' 

DATA  HA . HB/2 . 25E7 . 3 . 2oEft/ 

N-TYPE  OR  P-TYPEP 
IF (NP1.NE.1 ) GO  TO  10 

CAL.  N-TYPF  COEFF. 

ALFAro.O 

ESABS(EE) 

N-TYPE  COEFF.  = 0? 

IF (E.LT. 1 . 12E5)  RETURN 

~&L.  COEFF.  FOR  N-TYPE 

AlFA=EA*  { 1 . 0-Al.FATD*  ( T-300 . 0 ) ) *EXP L -EB/F ) 

ENFORCE  MINIMUM  VALUE  FOR  ALFA 
IFLALFA.LT. l.C)  ALFAsl.O 

RETURN 


10  CONTINUE 


EVALUATE  P-TYPF  COEFF 

ALFA=0.0 

E=ABS(EF) 

P-TYPE  COEFF.  = 0? 

IF (E.LT . 1 .93E5)  RETURN 

CAL  P-TYPE  COEFF. 

ALFA=HA* ( I . 0-ALFATD* ( T-300 . 0 ) ) *EXP ( -HB/E ) 

ENFORCE  MINIMUM  VALUE  FOR  ALFA 
IFLALFA.LT. 1.0)  ALFA=1.0 

RETURN 


END 


SUBROUTINE  LPLOT 
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ft 

i 


*TEMP(1> .LPLOT 
1 
2 

* s 

i i 

l cc 

9 C 

10  C 

11  c 

12  C 

1 3 C 
C 

is  c 

16  ( 

1/  c 


c 

c 


c 

c 


SUBROUTINE  LPLOT(NPTS*MPTS»NPLTS.MPLTS,MPLPT.NYCOR1»X, 
AXMIN, XMAX, Y » YMIN » YMAX * TITLE  * IGRID) 


SUB  LPLOT  GENERATES  A SINGLE  PAGF  LINE  PRINTER  PLOT  WHICH  MAY 
DISPLAY  A MAXIMUM  OF  26  DIFFERENT  CURVES.  THE  DIFFERENT  CURVES 

^ollTlols  COMMONHTOHMoSf  fhANA?fjF1CU^ERASEFcS6FDPWlVHEAN  ASTERISK 
DATA  POINTS  ARE  NOT  PLOTTED  ON,  OR  OUTSIDE  OF  THE  GRAPH  BOUNDARIES. 
THE  COORDINATE  BOUNDARIES  FOR  THF  PLOT  GENERATED  ARE  NOT 
ASSOCIATED  WITH  THE  MINIMUMS  AND  MAX? MUMS  FOR  THE  DATA  TO 


PLOTTED. 
A SINGLE 


HENCE#  LPLOT  HAS  A WINDOW  CAPABILITY. 
LINE  GRAPH  TITLE  OPTION. 


BE 


LPLOT  ALSO  HAS 


VARIABLE 

NPTS 

MPT? 


;cis  BMIh 


POINTS  OR  LFNGTH  OF 


NPLTS  - 
MPLTS  - 


MPLPT  - 
YCOR1  - 


■Stricted 

(Y.  (Y). 


BY  THE  DIMENSIONS 


NYC 


XHIN 


XMAX 


'MIN 


DEFINITIONS 

NUMBER  OF  ABSCISSA  DATA 
MAXIMUM  NUMBER  OF  ABSCl 
ABSCISSA  ARRAY. 

NUMBER  OF  CURVFS  TO  BE  PLOTTEI 
MAXIMUM  NUMBER  OF  CURVFS  AS  Rl 
OF  THE  DEPENDENT  VARIARIE  ARR 
MPTS*MPLTS 

ordinatf  data  format  sfntinel.  a zfro  vau 

STORED  COLUMN  WISE  AND  A NONZERO  VALUE  FOR 
STORED  ROW  WISF  IN  THE  DEPENDENT  VARIABLE 

- MINIMUM  ABSCISSA  VALUE  FOR  THE  PLOT,  BUT  NOT  NECESSARILY 

- WSt? *BUt  NOT 

- ffi  ^®paS,5CI|pYV^^0g8fNATVi  REVALUE,  

DIFFERENT  CURVFS  STORED  COLUMN  WISF  OR  ROW  WlSF. 

Y MUST  BE  DIMENSIONED  Y(MPTS*MPLTS)  FOR  THF  COLUMN 
“ RM£T  AND  Y(MPLTS,MPTS)  FOR  THE  ROW  WISE  FORMAT. 

VALUE  FOR  THE  PLOT  BUT  NOT  NECESSARILY 


VALUE  FOR  CURVES 

CURVES 

ARRAY  (Y). 
INDEPENDENT  VARIABLE, 


NFCESSARILY 


WITH  „ 
ARRAY 


WISE  FORMA 
MINIMUM  OR 


01 


YMAX 


AND 

NATE 

RDU 


TITLF  - 


_ m_3E  F°r  THE  PLOT  BUT 

THE  MINIMUM  OftDINATE  VALUE  FOR  THE  DATA. 

MAXIMUM  ORDINATE  VALUE  FOR  THE  PLOT  BUT  NOT  NECESSARILY 
THE  MAXIMUM  ORDINATE  VALUE  FOR  THE  DATA. 

A 22  ELEMENT  ONE  01 MENS! ON AL  ARRAY  FOR  THE  GRAPH  TITLE. 
THE  TITLE  SHOULD  BE  LEFT  ADJUSTIFIFD  WITH  6 CHARACTERS 
PER  WORD.  LPLOT  CENTERS  AND  PRINTS  THE  TITLE  TWO  LINES 
BELOW  THE  GRAPH.  SPECIFICATION  OF  THE  TITLE  IS  NOT 
NECESSARY,  BUT  THE  TITLF  ARRAY  MUST  BE  DIMENSIONED. 
IGRID  - GEID  SENTINEL.  GRID  LINES  ARE  NOT  PLOTTED  WHEN  I GRID 
IS  EQUAL  TO  ZERO. 


DIMENSION  CHAR( 2ft) ,XCMPTS) »Y (MPLPT) , ARRAY (51 , 101 ) 
A,XAXIS( 11) *YAXIS( 11) , TITLE ( 2? ) 


8 


ATA 

ATA 


IFLAG/O/ 

PLUS/ ♦ ♦ ♦ /ONF/ ♦ 1 • /AMIN/ • - * /Bl 


•PERIOD/* .»/ 
OATA  CHAR/* 


ANK/*  • /ZERO/OO /ASTR/ * A • / 


A*#*B*,*C*,»D*,»E*,»F*,*G***H*,*I»,*J*,*K*#*L» 


A,*M*.*N*,*0*,*P*,'Q*,*R*,*S*,*T*,*Ui,*V*,*W»,*X*,*Y*,*Z»/ 


GENERATE  GRAPH 


IfLAG.NE. 0 )G0*?0 


26 


IQ 


WRITE 

S^RR^Y 

ARRAY 
DO  10 
ARRAY 
ARRAY 


15 


CONTINUE 

CONTINUE 


vertical  AXES 
Jsl , 46, 5 
( J, 1)SPLUS 
( J, 101)=PLUS 
K=1 , 4 

(J+K,1)=0NF 
( J+K , 101 ) SONE 


IN  GRAPH  ARRAY 


JONTINUE 

ARRAY (51,1) =PLUS 
ARRAY (51,1 01 >=PLUS 


20 


WRITE  HORIZONTAL  AXES  IN  GRAPH  ARRAY 

DO  2S  J=1 *91  * 10 

ARRAy(1,J)=PLUS 

ARRAY (51 ,U) sPLUS 

DO  20  K=l,9 

ARRAY ( 1 , J+K ) =AMIN 

ARRAY ( 51  * J*K)=AMIN 

CONTINUE 


j 


SUBROUTINE  LPLOT  (Continued) 


25  CONTINUE 

26  CONTINUE 

BLANK  OUT  GRAPH  ARRAY 
DO  28  J=2» 100 


DO  27  K=2.50 
ARRA Y( K » J) rBLANK 

27  CONTINUE 

28  CONTINUE 

GENERATE  GRID 
IFUGRlO.EQ.OJGO  TO  29 
DO  31  J=6#46*5 
DO  31  K=2»100 
ARRAY (J'K)sPERIOD 

31  CONTINUE 

DO  32  J=11.9l»10 
DO  32  K=2»50 
ARRAY (K* JJrPERlOD 

32  CONTINUE 

29  CONTINUE 

EVALUATE  X AND  Y AXIS  SCALE 

DDX=(XMAX-XMIN)/10 

DDY=(YMAX-YMIN)/10 

DO  30  Jsl.ll 

XAXIS( J)SXM|N+DDX*( J«1 ) 

YAXIS(J)=YMIN  ♦ ODY  * { J-l ) 

30  CONTINUE 


BLANK 


VALUES 


(J-l  ) 


EVALUATE  X AND  Y INCREMENTS 
DX=  ( XMAX-XMIN) /l 00 
DY=  ( YMAX— YMIN) /50 

GENERATE  GRAPH  ARRAY 
IE(NYCORl.NE.O)  GO  TO  40 

ORDINATE  VALUES  STORED  COLUMN  WISE  IN  CALLING  PROGRAM 

DO  38  J=1.NPTS 

NX=  ( X ( J ) -XMI N > /DX4- 1 . 4999 

DO  37  K=J.MPLPT.MPTS 

KKS(K-J) /MPTS+1 

NY=( Y(K)-YMIN)/DY+1.4999 

IF  ( NX.lT .2. OR. NX. GT.  100 .OR.NY.LT •2*OR<1NY.GT.50 ) GO  TO  37 
IF  (ARRAY  (NY»NX) .EO.CHAR(KK) ) GO  TO  37 

IF  ( ARRAY(NY, NX) .EQ. BLANK. OR. ARRAY(NY»NX) .EG. PERIOD)  GO  TO  35 
ARRAY ( NY » NX ) r ASTR 
GO  TO  37 

35  ARRAY (NY 'NX)sCHAR(KK) 

37  CONTINUE 

38  CONTINUE 
GO  TO  49 


40  CONTINUE 


ORDINATE  values  storei 

DO  48  J= 1 » MPLPT  * MPLTS 
J Js ( J-l ) /MPLTS+ 1 
NXr(X( JJ)-XMIN)/DX  + 
KMAXs J+NPLTS-1 
DO  47  KsJtKMAX 
KK=K-J+1 


S STORED  ROW  WISE  IN  CALLING  PROGRAM 


1.4999 


ARRAY ( NY » NX ) sASTR 
GO  TO  47 

45  ARRAY(NY»NX)=CHAR(KK) 

47  CONTINUE 

48  CONTINUE 

49  CONTINUE 

WRITE  GRAPH  ARRAY 

WRITE (6»50)  (XAXTS(I).Isl.ll) 

50  FORMAT ( »1*/////.7X»11(1PE10.2>  #/) 
NNYS11 


.LT.^0R^NY.Gt.5°)  GO  TO  47 
RRAY(NY*NX) .EO. PERIOD)  GO  TO  45 


160 

161 

ill 

18 

166 

16? 

ill 

I?? 

B 

174 

IE 

177 

m 

iso 

1 82 

183 

184 

185 

186 

\ll 

it? 

191 


SUBROUTINE  LPLOT  (Continued) 


DO  70  JS51.6.-5 

WR1TF16.55)  YAXIS<NNY)t (ARRAY(J.I),lrl,101 1.YAXISCNNY) 
55  FORMAT (4X» 1PE8.P* IX »101A1»1X> 1PFR.2 ) 

NNY=NNY-1 
DO  65  K=l»4 
„ WRITF(6*60) (ARRAY 

60  FORMAT! 13X.101A1) 

65  CONTINUE 
70  CONTINUE 


B37 


<J-K.I).Isi,101> 


C 

C 


75 


80 


FORMAT! /*7X» 1 1 ( 1 PE10 . 2 ) ) 

WRITF  GRAPH  TITLE 

mm  ! J> .EQ. BLANK. OR. TITLE! J> ,EQ. ZERO. OR. J.GF.?2)G0  TO  85 
J— J+ 1 
GO  TO  80 
JUF 


85  CONTIN_ 

JJ=J-1 
JI=(?2-J>/2 
IF(JJ.LT.l)  GO  TO  95 

„ WRITF (6»90)  (BLANK#Is1.II).(TITLF(J).J=1.JJ» 
90  FORMAT ( //*  22A6  * / ) 

95  CONTINUE 

WRITF (6» 105) 

105  FORMAT! »1*> 

I FLAG: 1 

RETURN 

END 
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PROGRAM  TEMPTR10 


B38 


♦TEMP ( 1 ) 
1 


•TEMPTR1 J 


*****  TEMPTRIO  ***** 

DIMENSION  X< 101) .TEMPI  101) .TTEMP(lOl) »TEMPL< 101) .TEMPAR1 10. 101) 
DIMENSION  C< 101) ,E<101) .EEC  101) .PE(2.101> .BOPLOGI 101 ) .TITLE5122) 
DIMENSION  PLTIME(IO) • A ( 1 01 ,4) ,THKDPF<101 ) .THKHPF1101 ) 

DIMENSION  PTEMPX(SOO) ,PTIMEX(500) .TITLE1C22) ,TITLF2(?2) ,TITLE3(22> 
DIMENSION  PVOLT (4*500) .TITLE4<22> 

DIMENSION  TEMPH( 12.101). XXDH ( 14 ) 

I NTEGER  UPE » FLAGFB . PLOT . PDTP , PD IP . PDEP . PD VP . TFL AG . SCFLAG . SMAX . S 

NAMELIST  /DP ARMR/THKD » THKHX . THKHY » DDEN . HDEN . DSPFC . HSPEC . XDT . 

AXDh  . XLDEP . XMET . XEPDFP . XllDFP . XL . DOPL . DOPFP . DOPU » VEL » Al“ATD  . NP1 » 
ATHKDPF » THKHPF 

NAMEL 1ST  /DPARMtt/THKD » THKHX . THKHY . DDEN » HDEN . DSPFC . HSPEC . XDT . 


OPFP . SDOPU , VEL . ALF  A TO  . NP 1 


TIMEMX  . TMAX . TPMAX . I TSNMX . DTIME . 
GRID.TRMSMI , ITFRMX, ITLST. 


APLTImE.EMAXL.EMAXU.EMAXR.FINT  »UPF » AERMAX.EERMAX. ITCMAX. 
SlBB0,IDBl.IDB2,IDB3. IDB4.PDTP.PnTP. PDEP.PDVP 

DATA  T I TLE 1 /6HMA  X I Mli » 6HM  TEMP  * 6HER ATUR  . GHE  VFRS.6HUS  TIM. 

A6HE  # 16*1H  / 

DATA  TITLE2/6H0I0DE  .6HTEMPER.6HATURE  .GHVERSUS.6H  POSIT. 

S6HI0N  . 16*1H  / 

DATA  TITLE3/6HL0G  0F.6H  ELECT. 6HRIC  FI.6HELD  VE.GHRSUS  P. 
A6H0SITI0.6HN  .15*1H  / 

DATA  TITLE4/6HDI0DE  .6HV0L TA6.6HFS  VER.6HSUS  TI. 

&2HME. 17*1H  / 

DATA  TITLE5/6HL0G  OF.feH  DIODE. 6H  IMPUR.hHITY  CO.ftHNCENTR .6HATION 
A6HVERSUS.6H  POSIT. 6HI0N  . 13*lH  / 

ASSIGN  CONSTANTS 

Qsl.GE-19 

PERM=1.04E-12 

DEFINE  DIODE  AND  SIMULATION  PARAMETERS 

ASSIGN  DEFAULT  DIODE  PARAMETER  VALUFS 
THKDrl.O 

THKD  - WATTS/CM-K 


THKHX=0.46 

THKHY=0.46 

DDENr2 • 3 

HDEN=4*  0 

DSPEC-0.7 

HSPEC=0.79 

XDTsl .OE-4 

XDH=?50E-4 

XLDEPsO.O 

XMET=20.0E-4 

XEPDFPS40.0E-4 

XUDEP=0.0 

XL=4O.0E-4 

DOPL=1.0E17 

DOPEPS1.0E16 

DOPUrl . 0E16 


THKHX  - WATTS/CM-K 
THKHY  - WATTS/CM-K 
DDEN  - GM/CM3 
HDEN  - GM/CM3 
DSPEC  - J/GM-K 
HSPEC  - J/GM-K 
XDT  - CM 
XDH  - CM 
XLDEP  - CM 
XMET  - CM 
XEPDEP  - CM 
XUDEP  - CM 
XL  - CM 
DOPL  - CM- 3 
DOPEP  - CM- 3 


PROGRAM  TEMPTR10  (Continued) 


1 

I ft 
138 


C 

C 

c 

c 

c 

c 


c 

c 


c 

c 

c 


8 

c 

c 

c 

s 


VEL=1.0E7 
ALFATDs? • 5E-3 
NPl  = l 


DOPU  - CM— 3 

VEL  - CM/SEC 

ALFATD  - 

NPl  - 1 FOR  NP. 


0 FOR  PN 


INITIALIZE  THERMAL  CONDUCTIVITY  PERTURBATION  ARRAY 

DO  5 K— 1 » 101 

THKDPF (K)=1.0 

THKHPF ( K ) = 1 • 0 

CONTINUE 


ASSIGN  DEFAULT  SIMULATION  PARAMETERS. 
CUR=1.25E4 

CUR  - AMPS/CM2 

NND— 1 01 
NNHSO 
IBNDsO 
TIMEMXS1.0 


SPARM 


TMAXs  700.0 

TPMAX=800.0 

ITSNMXS200 

OTlMFSl.OE-9 

LTRNIN=2 

CtEMP=0 

LDHTFMrO 

PL0T=0 

PDTPsl 

PDIPrl 

PDEPsl 

PDVPrl 

PVOLTXS250.0 
I6RID=0 
TRMSMlslOOO.n 
ITERMXS10 
ITLST=10000 
IDHTM0S1 
MAXrlO 

THMAXSl.OE-4 
ITPRH=0 
EMAXL=1.0E5 
EMAXl  J=1 . 0E6 
EMAXR=1.0E5 
EINTslOO.O 
UPESO 

mum0.! 

=1.0E-4 


TIMFMX  - SEC 
TMAX  - DEG  K 
TPMAX  - DEG  K 

DTIME  - SEC 


eermaxs 

IDBOrO 


40 


IDB1=0 

m 

00  40  Jsl.10 

PLTIME(J)=100 

CONTINUE 

PLT I ME ( 1 ) S25  • OE-9 
PLTIME (2 >=50.06-9 
PLTlME(3)=75.0E-9 


' =150. OE-9 


PLTImE (5) = 
PLTIME (6) 


BEGIN  NEW  SIMULATION 
13  CONTINUE 


SKIP  TO  TOP 
WRITE (6» 9) 

9 FORMAT ( 1H1 ) 


OF  NEXT  PAGF  TO  BEGIN  NEW  SIMULATION 


read/write  simulation  parameters 

RE AD ( 5 » SPARM » END: 1000) 

WRITE<6, SPARM) 


PROGRAM  TEMPTR10  (Continued) 
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C 

C 

c 

c 

c 

c 

c 


c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


8 


READ  DIODE  PARAMETERS 
READ(5»0PARMR) 

EVALUATE  MOBILE  SPACE  CHARGE  ASSUMING  LIMITING  VELOCITY 

in  depletion  region 
SCONC=CUR/Q/VEL 

IF  EXCESSIVE  MOBILE  SPACE  CHARGE  SET  SCFLAG  TO  OMIT  THIS  SIMULATION 


SCFLAG=0 

IFtSCQNC.GE.DQF 

IF(SCONC.GE.DOF 


)PL ) SCFLAGsf 
_ >PEP>  SCFLAGsl 
IF(SCONC.GE.DOPU)  SCFALG=1 

EVALUATE  NET  SPACE  CHARGE 
SDOPL=DOPL-SCONC 
SDOPEP=DOPEP-SCONC 
SDOPU=DOPU-SCONC 

WRITE  DIODE  PARAMETERS 
WRITF(6.DPARMW) 

IF (SCFLAG. EQ.l)  WRITE<6,30> 

30  FORMAT ( /// » 5X » 45H*****  DRIVING  CURRENT  YIELDS  EXCESSIVE  MOBILE 
A39H  SPACE  CHARGE.  SIMULATION  OMITTED  ♦****) 

IF  SCFLAG. EQ.l  SKIP  TO  NEXT  SIMULATION 
IF(SCFLAG.EQ.l)  60  TO  13 

RESET  FLAG  FOR  HEADER  THERMAL  MODEL  TO  SIGNAL  NFW  DATA  SET 
IFLAG=0 

SET  FLAG  TO  WRITE  TRANSIENT  DATA  HEADER 
TFLAG=1 

INITIALIZE  PROGRAM  PARAMETERS 

NNDM1SNND-1 

NNDM?=NND-2 

NNHPj=NNH+l 

{?IS?5=0 

TIMEsO.O 


IPTINCsITSNMX/101+1 
PTEMPXt IPTIME >=0.0 
PTIMFX ( IPTIME ) =0 . 0 
DXSXL/NNDM1 
DXDX=DX*DX 
EFXDH=XOH 

IF (NNH.LT .3)  GO  TO  12 
EFXDH=XDH/NNH 
DO  1]  N=1.NNHP1 
XXDH(N) = (N-l ) *EFXOH 

11  CONTINUE 

12  CONTINUE 

GENERATE  POSITION  ARRAY  AND  INITIALIZE  TEMPFRATURF  ARRAY 
00  IS  1=1. NND 
X< I )=DX*< 1-1) 

TEMP ( I >=300,0 
TEMPL ( I ) =300 . 0 
TTEMP ( I ) =300 . 0 
TEMPH ( 1 » I ) =300 . 0 
15  CONTINUE 

™f  B!ni6um0sea68hb>Slue  |g»CF 

CALL 1 eK OF T NND? [>X?8lDEP  . XMET . XEPDEP . XUDEP . XL  . 

ASDOPL  » SDOPEP . SDOPU » FMAXL . E . I DBS • FLAGEB) 

£ALLpDOPLG  TO  GENERATE  LOG  OF  IMPURITY  CONCENTRATION  PROFILE 


CALL 


’LOTTING 

DOPLG(N 


LG ( NND . DX • XMET . XEPDEP . DOPL  * OOPU . DOPEP  » OOPLOG • I DBA ) 


BKDEP  TO  EVALUATE  DEPLETION  REGION  WIDTH  AND  F -FIELD 

PROFILE  FOR  THE  INITIAL  TEMPERATURE  PROFILE 
CALL  BKOEPL C NND . XLQEP » XMET . XEPDEP . XUDEP , XL . DX . 

ADOPL  * DOPEP . DOPU . ALFATD . EMAX . AERMA  X » EMAXL . 


PROGRAM  TEMPTR10  (Continued) 

AEMAXI I . E . TEMP , ALF  A IN  # VDEP . I TCM AX , 

A IDB1 . IDB2# IDB3.NP1 . SDOPL . SDOPEP. SDOPU) 

CALL  EFIELD  TO  EVALUATE  THE  BULK  REGION  E-FIFLD  PROFILES 

C ALLTE?  IELD  ( tEMpIlL^NoIl  TCMAX?Ff)5mAX  » XLDEP . XUDFP . CUR  »DOPL . DOPU. 
AOOPEP . XMET  . XEPOEP . NP 1 # VLBULK . VURl  ILK  # VO I OOE . FF . EMAXR » I0B4 # E I NT ) 

EVALUATE  TOTAL  OIODF  VOLTAGE 
VDIODESVLBULK+VOEP+VUBULK 

COMBINE  BULK  AND  DEPLETION  E-FIFI  D 
DO  1A  N=1»NND 
E (N)rAMAXl (E (N) »FE (N) ) 

STORF  INITIAL  E FIELD  FOR  PLOTTING 
PE(1.N)=E(N) 

18  CONTINUE 

CHECK  FOR  THERMAL  CONDUCTIVITY  PFRTURBATIONS 
DO  20  KKS1,NND 

IF((THKDP^(KK).NF.1.0).OR.(THKHPF(KK).NF.1.0))GO  TO  22 
20  CONTINUE 
GO  TO  100 

OUTPUT  THERMAL  CONDUCTIVITY  PERTURBATIONS 

22  CONTINUE 

23  rsiisitiB?  1 /// * 2 X • 4HN0DF  » 6X » 6HP0S I T I # 2H0N • 8X  » 6HD0P I NG  # 9X » 

A6HE  FlEL#2HD  » 7X , 6HTHKDPF » 9X  # 6HTHKHPF » // ) 

DETERMINE  DOPING  CONCENTRATIONS 

DO  24  KKK=1»NND 

IF ( X ( KKK ) ,LE.XMFT)DOP=DOPL 

IF( ( X(KKK) .GT.XMFT) .AND.(X(KKK) .LE.XEPDFP) IDOPxDOPEP 


I ?F ( 6» 2&?  KKK , X ( KKK ) ?|oP?£ < KKK  > , THKDPF ( KKK  ) , THKHPF ( KKK ) 

26  FORMAT(1X»I4.5X»F10.5»5X.F10.5#SX»E10.5.5X»F10.4#SX*F10.4) 
24  CONTINUE 

100  CONTINUE 

WRITE  TRANSIENT  DATA  HEAOFRF 
IF(TFLAG.NE.l)  GO  TO  110 

80  m&mm ! /// » T2 » 4HITSN » T9#  4HTIME » T18. 5HDTIME » T25 

A » 5HITERN. T33, 4HTRMS. T41 . 6HTEMPMX, T51 # 4HFMAX, T59, 6HVLBULK # 
AT69 . 4HVDEP . T77 » 6H VUBULK . TA6 * 6HVD I ODF # T95 . 5HXLDEP . 

ATI  04#  5HXUDEP#  T1 1 3*  6HALFA IN*  T121 » 1 HS» / ) 

RESET  TRANSIENT  DATA  HEADER  FLAG 
TFLAGsO 

TIMF  STEP  LOOP 


?ime=time+dtime 

ITSNsITSN+1 
I TERNS 0 

EVALUATE  POWER  DFNSITY 
DO  115  NS2.NNDM1 
C(N)sCUR*E(N) 

115  CONTINUE 


ITERATION  LOOP 
ITFRNSlTERN+1 

MiFT5ig}g^NE/!!t^TI0N 

BBS-THKHY/ ( DDEN*DSPEC*XDT*EFXDH ) 
DXDTRB=OXDX*DTlME*BB 


ME*AA 

ME*DXDX*BB 

TBB*300.0 


DO  130  NSI.NNDM? 


mmmmmrn 


PROGRAM  TEMPTR10  (Continued) 

A (N»  1 l=Al*THKDPF(N) 

A ( N . ? ) STHKHPF ( N 1 * A 2- ( DXDX+2 . 0*A ( N» 1 ) ) 

A (N, 31=A ( N, 1 1 

A ( N . U 1 s-DXDX* ( TFMPL  < N+ 1 1 +DTIME*C ( N+ l > / ( DDEN*DSPFC  1 1 
IFCNNH.GE.31  A(N»4)sA(N»4> +DxDTRR*TEMPH( 1 , N+ 1 ) 

IF(NNH.LT.31  A(N.4>=A(N,4>+A3 
10  CONTINUE 

evaluate  boundary  conditions 

I RND=0# CONSTANT  BOUNDARY  TEMPERATURES 

IFllRNO.NE.Ol  GO  TO  133 
A { 1 # 4 ) =A ( 1 . 4 ) -TEMP ( 1 ) *A C 1 . 1 ) 

A ( NNDM2 , 4 ) =A ( NN0M2  * 4 ) -TFMP ( NND ) * A ( NN0M2  * 3 ) 

133  CONTINUE 

IFC IRND.NF. l)GO  TO  135 

IBND=1 • BLOCKING  BOUNDARY  CONDITIONS 

A I NNDM2? £ I =1 ( NNOM? • £ 1 ♦ A ( NNOM? , 3 l 
135  CONTINUF 

SOLVF  SYSTEM  OF  LINFAR  FQUATIONS 

CALL  BANDA6CNNDM?,2.4,int ,101 ,4, A, TEMP, 1 ,PIVMIN) 

IF( IRND.NE.11  GO  TO  13B 

UPDATE  END  POINT  TEMPS  FOR  BLOCKING  BND 
TEMP(1)=TEMP(2) 

TEMP ( NNO ) STEMP { NNDM1 ) 

138  CONTINUE 

EVALUATE  MAXIMUM  TEMPERATURE  AND  RMS  TEMPERATURF  CHANGE 

TRMSrO.ODO 

TEMPMX=TEMP(1I 

9$m4!:?RMS*  ( VeMP ! Nl -TTEMP (N) ) **2 
TTEMP CN1=TEMP(N> 

TEMPmX=DMAX1(TEMPMX,TEMP<N) ) 

* WUZ'il \t  ( TRMS/NNDMP ) 

IF ( I TLS  T.E0.1)WRITE(6,150)ITSN,TTME»DTIME,I TFRN , TRMS » TEMPMX , 
AEMAX , VLBULK , VDEP , VURULK , VDl ODE , XI  DFP , XUDEP , AL FA IN , S 
>0  FORMAT ( I5.2E9.3, 15* 10E9.3, 13) 

IF(TRMS.LE.TRMSMl)  GO  TO  170 
IF ( jTERN.LT • ITERMX)  GO  TO  120 
WRITF (6,160 1 ITSN 

>0  FORMAT ( 29H  *****  CONVERGENCE  FAI|  URF  AT, 

17H  TIME  STEP 


. **  CONVERGENCE  FAIlURF  AT, 

A17H  TIME  STEP  NUMBER, 15, 7H  **•**) 

CONTINUE 

ITERATION  LOOP  COMPLETE,  ADVANCE  TO  NEXT  TIMF  STEP 

EVALUATE  SUBSTRATE  OR  HFADER  TEMPERATURE  PROFILFS 

QUASI-TwO  DIMENSIONAL  HFADER  THERMAL  MODFL 
IFCNNH.GT.3.ANO. IDHTMO.FQ.il  CALL  DHTEMPCIFLAG.NND, 


IFCNNH.GT.3.ANO. IDHTMO.FQ.il  CALL  DHTEMPl IFLAG.NND.NNH, 

« XOH , THKHY , HSPEC  # HDEN, DT I ME , TEMP , TEMPH 1 

two-dimensional  header  thermal  mooel.  headfr  tfmpfraturfs 

EVALUATED  COLUMN  WISE  (PROFILES  PERPENDICULAR  To  OIODE  AXIS! 


IF(NNH.GT.3.AND.IDhTM0.FQ. 
AXDH, THKHX , THKHY , HSPEC , HDEN 
AITPRH'IBND'S) 


ON  IS  r< 
2)  CALL 
»DTIMF ,' 


TEMP , TFMPHisMAX , PThMAX , 


TWO-DIMENSIONAL  HEADER  THERMAL  MODEL.  HEADFR  TEMPERATURES 

Bruise  Yc8b®TE-isfEA??ffiftifspfs?^(iwg  is  ?js??e  axis> 

ARE  evaluated. 

IF(NNH.GT.3.AND.IDHTM0.FQ.31  CALL  DHT2D1 ( IFLAG.NND.NNH, 

A XL . XDH , THKHX , THKHY , HSPEC . HDEN » DT IMF  * TEMP . TEMPH . SMAX . 


PROGRAM  TEMP TRIO  (Continued) 


B43 


&OTHMAX * I TPRH , IBND * S I 


two-dimensional  header  thermal  model,  headfr  tfmperaturfs 

EVALUATED  ROW  WISE  (PROFILES  PARALLEL  TO  DIODE  AXIS) 

WITH  NEWLY  COMPUTED  TEMPERATURES  INCORPERATfd  AS  THFY 


ARE  EVALUATED. 

IF ( NNH.GT . 3. AND. IDHTM0.FQ.4 ) CALL  DHT2D2 ( IFLAG.NND.NNH. 

A XL  * XDH . THKHX . THKHY . HSPEC . HDEN  # DT  T ME . TEMP . TEMPH . SMAX . 

SDTHMAX . I TPRH . IBND • S ) 

STORF  PRESENT  DIODE  TEMPERATURE  PROFILE  FOR  NEXT  SFT 
OF  ITERATIONS 
DO  175  N=2.NNDM1 
TEMPL(N)=TEMP(N) 

175  CONTINUE 

IF(UPE.LE.O)  GO  TO  185 

UPDATE  DEP  REG  BNDS  AND  E-FIELD  FOR  NEXT  TIMF  STEP 
CALL  BKDEPL ( NND . XLDEP . XMET » XEPDFP . 

AXUDEP » XL  t DX » DOPL . OOPEP » DOPU*  ALF ATD . EMAX . AERMAX . 
«EMAXL»EmAXU»E#TEmP*ALFATN»VDEP» ITCMAX. 

A IDB 1 . IDB2 , I DB3 . NP 1 , SDOPL , SDOPEP . SDOPU ) 

UPDATE  BULK  REG  E-FIELD  FOR  NEXT  TIME  STEP 

EVALUATE  TOTAL  DIODE  VOLTAGE 
vdiode=vlbulk+vdep+vubulk 

COMBINE  BULK  ANn  DEPLETION  RFGION  E-FIELDS 
DO  179  N=1.NND 
E(N)=AMAX1 (E(N) .FE(N) ) 

179  CONTINUE 

185  CONTINUE 

SAVE  MAXIMUM  TEMPERATURF  AND  DIODE  VOLTAGES  VERSUS  TIME 
IFTSI=0 

IF(ITSN-ITSN/IPTINC*IPTINC.NE.0)  GO  TO  17B 
IPTIME=IPTIME4-1 
PTEMPX (IPTIME ) STFMPMX 
PVOLT ( 1 , IPTIME ) rVLBULK 
PVOLT (2# IPTIME ) =VDEP 
PV0LT(3»!pTIME)=VUBULK 
PVOLT ( 4 . IPT I ME ) = VDIODE 
PTIMFXI IPTIME >=TIME 
IFTSisl 
178  CONTINUE 

OUTPUT  SECTION. 

LIST  TRANSIENT  DATA 

feffiSSJ.TSN/LTRN.N.LTRN.N  •NE.O ) GO  TO  180 
write (6# 150)  itsn.time.dtime.itfrn.trms.tempmx. 

AEMAX . VLBULK , VDEP , VUBULK . VDIODE . XI  DEP , XUOFP , ALFA IN , S 

m £ontinG£ 

TERMINATE  EXECUTION  IF  FITHER  max  TIME.  MAX  TIMF  STFP 
COUNT  OR  MAX  TEMP  IS  EXCEFDED* 

IF(TIME.GE.TIMEMX.OR.ITSN.GE.ITSNMX.OR.TEMPMX.GF.TMAX)  GO  TO  19( 
IF(TTME.LT.PLTIMFCKPLOT) ) GO  TO  230 


TMAX)  GO  TO  190 


STORF  TEMPERATURF  PROFILES 
190  CONTINUE 

KPLOTsKPLOT ♦ 1 

IF (KPLOT . LE . 11 ) GO  TO  210 

WRITF16.200) 

200  FORMAT (/* 40H  •**♦*  MORE  THAN  10  TEMPERATURE  PROFILFS. 
A16H  REQUESTED  ****♦. /.28H  *♦**♦  SIMULATION  TERMINATED. 
A27H  *****,/> 


A27H 

KPLOT=KPLOT-l 
GO  TO  240 

210  CONTINUE 


m 

S3  i 

540 

541 

is! 

545 

3S$ 

548 


* 


PROGRAM  TEMPTR10  (Continued) 

DEFINE  PLOT  TIME  FOR  EXTRA  PLOT 

IF(TIME.LT.PLTIMF(KPL0T-1) ) PLTIME(KPL0T-I)=TIMF 
WRITE  TIME  STEP  DATA  FOR  PLOT  TIME 

IF(ICONSN.NE.I)  WRITF(6#150)  ITSN, TIME, OTIME * ITERN, TRMS, TEMPMX , 
ftEMAX . VLBULK  » VDEP . VUBOLK • VD I ODE * XLDEP » XUDEP  # ALFA I N » S 

STORF  MAX  TEMP,  DIODE  VOLTAGES  AND  TIME  FOR  PLOTTING 
IF  IFTSI.EQ.1  VALUES  ALREADY  STORED  FOR  LAST  TIME  STEP 
IF(IFTSI.EQ.I)  GO  TO  215 
IPTIMErlPTIME+1 
PTEMPXf IPTIME )=TFMPMX 
PVOLT ( 1 , IPTIME ) sVLBULK 
PVOLT ( 2 , IPTIME ) rVDEP 
PVOLT (3, IPTIME )=VUBULK 
PVOLT (4, IPTIME)=VDI0DE 
PTIMEXt IPTIME)=TTME 
215  CONTINUE 

STORE  DIODE  TEMPERATURE  PROFILES 


DO  220  JS1.NND 
TEMPAR ( KPLOT-1 , J ) =TEMP ( J ) 

220  CONTINUE 

IF(NNH.LT,3.0R.LDHTEM.NF. J ) GO  TO  230 

WRITE  DIODE  AND  HEADER  TEMPERATURE  PROFILES 
WRITE (6,223)  ITSN. TIME » ( XXDH (N) ,N=1,NNHP1 ) 

223  FORMAT ( 9H1  ITSN  = , I5» 5X , 7HTIME  = ,Efl.3,//, 
AT12»5HXH  s » 14EA.3) 

WRITE (6, 224) 

224  FORMAT ( / , T7 , 1HN , T 1 2 , 2HXO , T18 , 5HDTEMP , T2G , 5HHTEMP , / ) 
DO  228  JS1.NND 

WRITE (6, 226)  J, X ( J) . TEMP ( J) , (TEMPH(K,J) ,K=1.NNH) 

226  FORMAT (I8,15E8.3) 

228  CONTINUE 

SET  TFLAG  FOR  NEW  TRANSIENT  DATA  HEADER 
TFLAG=1 


230  CONTINUE 


TERMINATE  EXECUTION  IF  FITHFR  MAX  TIME » M 
COUNT  OR  MAX  TEMP  IS  EXCEEDED 
IF(TTME.LT.TIMEMX.AND.ITSN.LT.ITSNMX.AND. 


MAX  TIMF  STFP 


TEMPMX.LT.TMAX)  GO  TO  100 


240  CONTINUE 


TERMINATE  EXECUTION  AND  OUTPUT  SIMULATION  RFSULTS 
PLOT  MAXIMUM  TEMPERATURE  VERSUS  TIME 

CALL  LPLOT (IPTIME. 500, 1,1, 500, 0,PTIMEX,0.n.PTIMFX( IPTIME ) ,PTEMPX 
8,300.0, TPMAX. TITLE 1 . IGRTD ) 

LIST  AND/OR  PLOT  TEMPERATURE  PROEILFS 


KPLOT=KPLOT-l 

IF(LTEMP.NE.l) 


GO  TO  280 


LIST  DIODE  TEMPERATURE  PROFILES 


WRITE(6.250)  (PLTIME(J) ,J=l.KPLOT) 

250  FORMAT (1H1,T8,1HN.T1 3. 1HX.T19.5HTIME1.T27.5HTIMF2 
8 , T35, 5HTIME3. T43, 5HTIME4, T51 , 5HTTMF5. T59. 5HTTME6 
« , T67 , 5HTIME7 , T75, 5HTIMtft, TgS. 5HTTMP 9, T91 . 6HTTME 1 0 
8,T99,6HTIME11,T107.6HTIME12,T115,6HTIME13,T1?3 
8.6HT IME14,//, 16X, 14E8.3) 

WRITf (6,251 ) 

251  FORMAT ( ) 

WRI?F(6?2i&VN8,X(K) . ( TEMPAR ( J.K) , Jrl.KPLOT) 

260  FORMAT ( IX, 17, 15EA.3) 

270  CONTINUE 

280  CONTINUE 

I F ( ( PLOT .NE.l).AND.(PDTP.NE.l)) AO  TO  283 

PLOT  DIODE  TEMPERATURE  PROFILES 

CALL  LPLOT (NNO. 101. KPLOT. 10, 1010, 1 ,X,0.0,XL. 


PROGRAM  TEMPTR10  (Continued) 


B45 


560 

561 

562 

563 

564 

565 


568 

569 

570 

571 

m 

575 

578 

579 

580 

m 

583 

584 

585 

586 

m 

589 

590 

591 


C 

C 

c 

c 


4 TEMP A R .300.0. TPMAX » TITLF2. IGRID ) 


‘(NND.101.1 .1. tOl *0.X.0,0.XL»DOPLOG.10.0.20.0.TlTLF5»lGRI 


C 

C 


m 

596 


C 

C 


C 

C 


C 

C 

C 


283  I F ( ( PLOT . NE . 1 ) . AND . ( PD I P . NE • 1 ) ) GO  TO  287 
PLOT  DjogEfIMPURlTY  PROFILE. 

AD) 

287  IF((PLOT.NE.l).AND.(PDEP.NE.l))GO  TO  350 
STORF  FINAL  E FIELD  PROFILE  FOR  PL0TTIN6 
DO  290  N=1.NND 
PE(2.N)=E(N) 

290  CONTINUE 

6ENERATE  LOG  OF  F FIELD  PLOTTING  ARRAY 
DO  310  M=1.2 
DO  320  N=!#NND 

IF(PF(M.N) .GE.1.0E-2)  TPE  =ALOG10(PE<M.N> ) 

IF(PF<M.N) .LT.1.0E-P)  TPE  =-2.0 
PE ( M#  N) =TPE 
320  CONTINUE 
310  CONTINUE 

PLOT  LOG  OF  FIRST  AND  LAST  E FlFl D PROFILES  VERSUS  POSITION 
CALL  LPLOT(NND. 101»2.2»202»l.X.n.O.XL»PE.-2.n»8.0. 
*TITLF3»IGRID) 

350  IF C ( PLOT .NE* 1 ) .AND • (PDVP.NE • 1 ) ) GO  TO  500 

PLOT  DIODE  VOLTAGES  VERSUS  TIME 

CALL  LPLOTt IPTIME » 500 »4. 4# 404 »l#PTIMEX»n.O» 

APTIMFX ( IPTIME ) . PVOLT . 0.0. PVOLTX . TITLE4 . IGRID ) 

500  CONTINUE 

GO  TO  13 

1000  STOP 
END 


